Comparative cellular analysis of motor cortex in human, marmoset and mouse

Predicting cell population-specific gene expression from genomic sequence

Most regulatory elements, especially enhancer sequences, are cell population-specific. One could even argue that a distinct set of regulatory elements is what defines a cell population. However, discovering which non-coding regions of the DNA are essential in which context, and as a result, which genes are expressed, is a difficult task. Some computational models tackle this problem by predicting gene expression directly from the genomic sequence. These models are currently limited to predicting bulk measurements and mainly make tissue-specific predictions. Here, we present a model that leverages single-cell RNA-sequencing data to predict gene expression. We show that cell population-specific models outperform tissue-specific models, especially when the expression profile of a cell population and the corresponding tissue are dissimilar. Further, we show that our model can prioritize GWAS variants and learn motifs of transcription factor binding sites. We envision that our model can be useful for delineating cell population-specific regulatory elements.

Single-cell multiplex chromatin and RNA interactions in aging human brain

The dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA (caRNA) 1-3. Chromatin complex compositions change during cellular differentiation and aging, and are expected to be highly heterogeneous among terminally differentiated single cells 4-7. Here we introduce the Multi-Nucleic Acid Interaction Mapping in Single Cell (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression, and RNA-chromatin associations within individual nuclei. Applied to 14 human frontal cortex samples from elderly donors, MUSIC delineates diverse cortical cell types and states. We observed the nuclei exhibiting fewer short-range chromatin interactions are correlated with an "older" transcriptomic signature and with Alzheimer's pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci (cis eQTLs) and a promoter tends to be the cell type where these cis eQTLs specifically affect their target gene's expression. Additionally, the female cortical cells exhibit highly heterogeneous interactions between the XIST non-coding RNA and Chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploring chromatin architecture and transcription at cellular resolution in complex tissues.

Genetic architecture of the structural connectome

Myelinated axons form long-range connections that enable rapid communication between distant brain regions, but how genetics governs the strength and organization of these connections remains unclear. We perform genome-wide association studies of 206 structural connectivity measures derived from diffusion magnetic resonance imaging tractography of 26,333 UK Biobank participants, each representing the density of myelinated connections within or between a pair of cortical networks, subcortical structures or cortical hemispheres. We identify 30 independent genome-wide significant variants after Bonferroni correction for the number of measures studied (126 variants at nominal genome-wide significance) implicating genes involved in myelination (SEMA3A), neurite elongation and guidance (NUAK1, STRN, DPYSL2, EPHA3, SEMA3A, HGF, SHTN1), neural cell proliferation and differentiation (GMNC, CELF4, HGF), neuronal migration (CCDC88C), cytoskeletal organization (CTTNBP2, MAPT, DAAM1, MYO16, PLEC), and brain metal transport (SLC39A8). These variants have four broad patterns of spatial association with structural connectivity: some have disproportionately strong associations with corticothalamic connectivity, interhemispheric connectivity, or both, while others are more spatially diffuse. Structural connectivity measures are highly polygenic, with a median of 9.1 percent of common variants estimated to have non-zero effects on each measure, and exhibited signatures of negative selection. Structural connectivity measures have significant genetic correlations with a variety of neuropsychiatric and cognitive traits, indicating that connectivity-altering variants tend to influence brain health and cognitive function. Heritability is enriched in regions with increased chromatin accessibility in adult oligodendrocytes (as well as microglia, inhibitory neurons and astrocytes) and multiple fetal cell types, suggesting that genetic control of structural connectivity is partially mediated by effects on myelination and early brain development. Our results indicate pervasive, pleiotropic, and spatially structured genetic control of white-matter structural connectivity via diverse neurodevelopmental pathways, and support the relevance of this genetic control to healthy brain function.

Drug targeting in psychiatric disorders - how to overcome the loss in translation?

In spite of major efforts and investment in development of psychiatric drugs, many clinical trials have failed in recent decades, and clinicians still prescribe drugs that were discovered many years ago. Although multiple reasons have been discussed for the drug development deadlock, we focus here on one of the major possible biological reasons: differences between the characteristics of drug targets in preclinical models and the corresponding targets in patients. Importantly, based on technological advances in single-cell analysis, we propose here a framework for the use of available and newly emerging knowledge from single-cell and spatial omics studies to evaluate and potentially improve the translational predictivity of preclinical models before commencing preclinical and, in particular, clinical studies. We believe that these recommendations will improve preclinical models and the ability to assess drugs in clinical trials, reducing failure rates in expensive late-stage trials and ultimately benefitting psychiatric drug discovery and development.

A concerted neuron-astrocyte program declines in ageing and schizophrenia

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA sequencing to analyse the prefrontal cortex of 191 human donors aged 22-97 years, including healthy individuals and people with schizophrenia. Latent-factor analysis of these data revealed that, in people whose cortical neurons more strongly expressed genes encoding synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the synaptic neuron and astrocyte program (SNAP). In schizophrenia and ageing-two conditions that involve declines in cognitive flexibility and plasticity1,2-cells divested from SNAP: astrocytes, glutamatergic (excitatory) neurons and GABAergic (inhibitory) neurons all showed reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy people of similar age, may underlie many aspects of normal human interindividual differences and may be an important point of convergence for multiple kinds of pathophysiology.

Noninvasive focal transgene delivery with viral neuronal tracers in the marmoset monkey

We establish a reliable method for selectively delivering adeno-associated viral vectors (AAVs) across the blood-brain barrier (BBB) in the marmoset without the need for neurosurgical injection. We focally perturbed the BBB (∼1 × 2 mm) in area 8aD of the frontal cortex in four adult marmoset monkeys using low-intensity transcranial focused ultrasound aided by microbubbles. Within an hour of opening the BBB, either AAV2 or AAV9 was delivered systemically via tail-vein injection. In all four marmosets, fluorescence-encoded neurons were observed at the site of BBB perturbation, with AAV2 showing a sparse distribution of transduced neurons when compared to AAV9. The results are compared to direct intracortical injections of anterograde tracers into area 8aD and similar (albeit sparser) long-range connectivity was observed. With evidence of transduced neurons specific to the region of BBB opening as well as long-distance tracing, we establish a framework for focal noninvasive transgene delivery to the marmoset brain.

Molecularly stratified hypothalamic astrocytes are cellular foci for obesity

Astrocytes safeguard the homeostasis of the central nervous system1,2. Despite their prominent morphological plasticity under conditions that challenge the brain's adaptive capacity3-5, the classification of astrocytes, and relating their molecular make-up to spatially devolved neuronal operations that specify behavior or metabolism, remained mostly futile6,7. Although it seems unexpected in the era of single-cell biology, the lack of a major advance in stratifying astrocytes under physiological conditions rests on the incompatibility of 'neurocentric' algorithms that rely on stable developmental endpoints, lifelong transcriptional, neurotransmitter, and neuropeptide signatures for classification6-8 with the dynamic functional states, anatomic allocation, and allostatic plasticity of astrocytes1. Simplistically, therefore, astrocytes are still grouped as 'resting' vs. 'reactive', the latter referring to pathological states marked by various inducible genes3,9,10. Here, we introduced a machine learning-based feature recognition algorithm that benefits from the cumulative power of published single-cell RNA-seq data on astrocytes as a reference map to stepwise eliminate pleiotropic and inducible cellular features. For the healthy hypothalamus, this walk-back approach revealed gene regulatory networks (GRNs) that specified subsets of astrocytes, and could be used as landmarking tools for their anatomical assignment. The core molecular censuses retained by astrocyte subsets were sufficient to stratify them by allostatic competence, chiefly their signaling and metabolic interplay with neurons. Particularly, we found differentially expressed mitochondrial genes in insulin-sensing astrocytes and demonstrated their reciprocal signaling with neurons that work antagonistically within the food intake circuitry. As a proof-of-concept, we showed that disrupting Mfn2 expression in astrocytes reduced their ability to support dynamic circuit reorganization, a time-locked feature of satiety in the hypothalamus, thus leading to obesity in mice. Overall, our results suggest that astrocytes in the healthy brain are fundamentally more heterogeneous than previously thought and topologically mirror the specificity of local neurocircuits.

CHRNA5 links chandelier cells to severity of amyloid pathology in aging and Alzheimer's disease

Changes in high-affinity nicotinic acetylcholine receptors are intricately connected to neuropathology in Alzheimer's Disease (AD). Protective and cognitive-enhancing roles for the nicotinic α5 subunit have been identified, but this gene has not been closely examined in the context of human aging and dementia. Therefore, we investigate the nicotinic α5 gene CHRNA5 and the impact of relevant single nucleotide polymorphisms (SNPs) in prefrontal cortex from 922 individuals with matched genotypic and post-mortem RNA sequencing in the Religious Orders Study and Memory and Aging Project (ROS/MAP). We find that a genotype robustly linked to increased expression of CHRNA5 (rs1979905A2) predicts significantly reduced cortical β-amyloid load. Intriguingly, co-expression analysis suggests CHRNA5 has a distinct cellular expression profile compared to other nicotinic receptor genes. Consistent with this prediction, single nucleus RNA sequencing from 22 individuals reveals CHRNA5 expression is disproportionately elevated in chandelier neurons, a distinct subtype of inhibitory neuron known for its role in excitatory/inhibitory (E/I) balance. We show that chandelier neurons are enriched in amyloid-binding proteins compared to basket cells, the other major subtype of PVALB-positive interneurons. Consistent with the hypothesis that nicotinic receptors in chandelier cells normally protect against β-amyloid, cell-type proportion analysis from 549 individuals reveals these neurons show amyloid-associated vulnerability only in individuals with impaired function/trafficking of nicotinic α5-containing receptors due to homozygosity of the missense CHRNA5 SNP (rs16969968A2). Taken together, these findings suggest that CHRNA5 and its nicotinic α5 subunit exert a neuroprotective role in aging and Alzheimer's disease centered on chandelier interneurons.

Calcium-Sensitive Subthreshold Oscillations and Electrical Coupling in Principal Cells of Mouse Dorsal Cochlear Nucleus

In higher sensory brain regions, slow oscillations (0.5-5 Hz) associated with quiet wakefulness and attention modulate multisensory integration, predictive coding, and perception. Although often assumed to originate via thalamocortical mechanisms, the extent to which subcortical sensory pathways are independently capable of slow oscillatory activity is unclear. We find that in the first station for auditory processing, the cochlear nucleus, fusiform cells from juvenile mice (of either sex) generate robust 1-2 Hz oscillations in membrane potential and exhibit electrical resonance. Such oscillations were absent prior to the onset of hearing, intrinsically generated by hyperpolarization-activated cyclic nucleotide-gated (HCN) and persistent Na+ conductances (NaP) interacting with passive membrane properties, and reflected the intrinsic resonance properties of fusiform cells. Cx36-containing gap junctions facilitated oscillation strength and promoted pairwise synchrony of oscillations between neighboring neurons. The strength of oscillations were strikingly sensitive to external Ca2+, disappearing at concentrations >1.7 mM, due in part to the shunting effect of small-conductance calcium-activated potassium (SK) channels. This effect explains their apparent absence in previous in vitro studies of cochlear nucleus which routinely employed high-Ca2+ extracellular solution. In contrast, oscillations were amplified in reduced Ca2+ solutions, due to relief of suppression by Ca2+ of Na+ channel gating. Our results thus reveal mechanisms for synchronous oscillatory activity in auditory brainstem, suggesting that slow oscillations, and by extension their perceptual effects, may originate at the earliest stages of sensory processing.

Transcriptional cartography integrates multiscale biology of the human cortex

The cerebral cortex underlies many of our unique strengths and vulnerabilities, but efforts to understand human cortical organization are challenged by reliance on incompatible measurement methods at different spatial scales. Macroscale features such as cortical folding and functional activation are accessed through spatially dense neuroimaging maps, whereas microscale cellular and molecular features are typically measured with sparse postmortem sampling. Here, we integrate these distinct windows on brain organization by building upon existing postmortem data to impute, validate, and analyze a library of spatially dense neuroimaging-like maps of human cortical gene expression. These maps allow spatially unbiased discovery of cortical zones with extreme transcriptional profiles or unusually rapid transcriptional change which index distinct microstructure and predict neuroimaging measures of cortical folding and functional activation. Modules of spatially coexpressed genes define a family of canonical expression maps that integrate diverse spatial scales and temporal epochs of human brain organization - ranging from protein-protein interactions to large-scale systems for cognitive processing. These module maps also parse neuropsychiatric risk genes into subsets which tag distinct cyto-laminar features and differentially predict the location of altered cortical anatomy and gene expression in patients. Taken together, the methods, resources, and findings described here advance our understanding of human cortical organization and offer flexible bridges to connect scientific fields operating at different spatial scales of human brain research.

Progress in single-cell multimodal sequencing and multi-omics data integration

With the rapid advance of single-cell sequencing technology, cell heterogeneity in various biological processes was dissected at different omics levels. However, single-cell mono-omics results in fragmentation of information and could not provide complete cell states. In the past several years, a variety of single-cell multimodal omics technologies have been developed to jointly profile multiple molecular modalities, including genome, transcriptome, epigenome, and proteome, from the same single cell. With the availability of single-cell multimodal omics data, we can simultaneously investigate the effects of genomic mutation or epigenetic modification on transcription and translation, and reveal the potential mechanisms underlying disease pathogenesis. Driven by the massive single-cell omics data, the integration method of single-cell multi-omics data has rapidly developed. Integration of the massive multi-omics single-cell data in public databases in the future will make it possible to construct a cell atlas of multi-omics, enabling us to comprehensively understand cell state and gene regulation at single-cell resolution. In this review, we summarized the experimental methods for single-cell multimodal omics data and computational methods for multi-omics data integration. We also discussed the future development of this field.

Dictionary learning for integrative, multimodal and scalable single-cell analysis

Mapping single-cell sequencing profiles to comprehensive reference datasets provides a powerful alternative to unsupervised analysis. However, most reference datasets are constructed from single-cell RNA-sequencing data and cannot be used to annotate datasets that do not measure gene expression. Here we introduce 'bridge integration', a method to integrate single-cell datasets across modalities using a multiomic dataset as a molecular bridge. Each cell in the multiomic dataset constitutes an element in a 'dictionary', which is used to reconstruct unimodal datasets and transform them into a shared space. Our procedure accurately integrates transcriptomic data with independent single-cell measurements of chromatin accessibility, histone modifications, DNA methylation and protein levels. Moreover, we demonstrate how dictionary learning can be combined with sketching techniques to improve computational scalability and harmonize 8.6 million human immune cell profiles from sequencing and mass cytometry experiments. Our approach, implemented in version 5 of our Seurat toolkit ( http://www.satijalab.org/seurat ), broadens the utility of single-cell reference datasets and facilitates comparisons across diverse molecular modalities.

Deciphering microglia phenotypes in health and disease

Microglia are the major immune cells of the central nervous system (CNS) that perform numerous adaptive functions required for normal CNS development and homeostasis but are also linked to neurodegenerative and behavioral diseases. Microglia development and function are strongly influenced by brain environmental signals that are integrated at the level of transcriptional enhancers to drive specific programs of gene expression. Here, we describe a conceptual framework for how lineage-determining and signal-dependent transcription factors interact to select and regulate the ensembles of enhancers that determine microglia development and function. We then highlight recent findings that advance these concepts and conclude with a consideration of open questions that represent some of the major hurdles to be addressed in the future.

Integrating single-cell transcriptomics with cellular phenotypes: cell morphology, Ca2+ imaging and electrophysiology

I review recent technological advancements in coupling single-cell transcriptomics with cellular phenotypes including morphology, calcium signaling, and electrophysiology. Single-cell RNA sequencing (scRNAseq) has revolutionized cell type classifications by capturing the transcriptional diversity of cells. A new wave of methods to integrate scRNAseq and biophysical measurements is facilitating the linkage of transcriptomic data to cellular function, which provides physiological insight into cellular states. I briefly discuss critical factors of these phenotypical characterizations such as timescales, information content, and analytical tools. Dedicated sections focus on the integration with cell morphology, calcium imaging, and electrophysiology (patch-seq), emphasizing their complementary roles. I discuss their application in elucidating cellular states, refining cell type classifications, and uncovering functional differences in cell subtypes. To illustrate the practical applications and benefits of these methods, I highlight their use in tissues with excitable cell-types such as the brain, pancreatic islets, and the retina. The potential of combining functional phenotyping with spatial transcriptomics for a detailed mapping of cell phenotypes in situ is explored. Finally, I discuss open questions and future perspectives, emphasizing the need for a shift towards broader accessibility through increased throughput.

Non-CG DNA methylation and MeCP2 stabilize repeated tuning of long genes that distinguish closely related neuron types

The extraordinary diversity of neuron types in the mammalian brain is delineated at the highest resolution by subtle gene expression differences that may require specialized molecular mechanisms to be maintained. Neurons uniquely express the longest genes in the genome and utilize neuron-enriched non-CG DNA methylation (mCA) together with the Rett syndrome protein, MeCP2, to control gene expression, but the function of these unique gene structures and machinery in regulating finely resolved neuron type-specific gene programs has not been explored. Here, we employ epigenomic and spatial transcriptomic analyses to discover a major role for mCA and MeCP2 in maintaining neuron type-specific gene programs at the finest scale of cellular resolution. We uncover differential susceptibility to MeCP2 loss in neuronal populations depending on global mCA levels and dissect methylation patterns and intragenic enhancer repression that drive overlapping and distinct gene regulation between neuron types. Strikingly, we show that mCA and MeCP2 regulate genes that are repeatedly tuned to differentiate neuron types at the highest cellular resolution, including spatially resolved, vision-dependent gene programs in the visual cortex. These repeatedly tuned genes display genomic characteristics, including long length, numerous intragenic enhancers, and enrichment for mCA, that predispose them to regulation by MeCP2. Thus, long gene regulation by the MeCP2 pathway maintains differential gene expression between closely-related neurons to facilitate the exceptional cellular diversity in the complex mammalian brain.

Cell type specificity of mosaic chromosome 1q gain resolved by snRNA-seq in a case of epilepsy with hyaline protoplasmic astrocytopathy

Introduction

Mosaic gain of chromosome 1q (chr1q) has been associated with malformation of cortical development (MCD) and epilepsy. Hyaline protoplasmic astrocytopathy (HPA) is a rare neuropathological finding seen in cases of epilepsy with MCD. The cell-type specificity of mosaic chr1q gain in the brain and the molecular signatures of HPA are unknown.

Methods

We present a child with pharmacoresistant epilepsy who underwent epileptic focus resections at age 3 and 5 years and was found to have mosaic chr1q gain and HPA. We performed single-nuclei RNA-sequencing (snRNA-seq) of brain tissue from the second resection.

Results

snRNA-seq showed increased expression of chr1q genes specifically in subsets of neurons and astrocytes. Differentially expressed genes associated with inferred chr1q gain included AKT3 and genes associated with cell adhesion or migration. A subpopulation of astrocytes demonstrated marked enrichment for synapse-associated transcripts, possibly linked to the astrocytic inclusions observed in HPA.

Discussion

snRNA-seq may be used to infer the cell type-specificity of mosaic chromosomal copy number changes and identify associated gene expression alterations, which in the case of chr1q gain may involve aberrations in cell migration. Future studies using spatial profiling could yield further insights on the molecular signatures of HPA.

Layer 1 neocortex: Gating and integrating multidimensional signals

Layer 1 (L1) of the neocortex acts as a nexus for the collection and processing of widespread information. By integrating ascending inputs with extensive top-down activity, this layer likely provides critical information regulating how the perception of sensory inputs is reconciled with expectation. This is accomplished by sorting, directing, and integrating the complex network of excitatory inputs that converge onto L1. These signals are combined with neuromodulatory afferents and gated by the wealth of inhibitory interneurons that either are embedded within L1 or send axons from other cortical layers. Together, these interactions dynamically calibrate information flow throughout the neocortex. This review will primarily focus on L1 within the primary sensory cortex and will use these insights to understand L1 in other cortical areas.

Concerted neuron-astrocyte gene expression declines in aging and schizophrenia

Human brains vary across people and over time; such variation is not yet understood in cellular terms. Here we describe a striking relationship between people's cortical neurons and cortical astrocytes. We used single-nucleus RNA-seq to analyze the prefrontal cortex of 191 human donors ages 22-97 years, including healthy individuals and persons with schizophrenia. Latent-factor analysis of these data revealed that in persons whose cortical neurons more strongly expressed genes for synaptic components, cortical astrocytes more strongly expressed distinct genes with synaptic functions and genes for synthesizing cholesterol, an astrocyte-supplied component of synaptic membranes. We call this relationship the Synaptic Neuron-and-Astrocyte Program (SNAP). In schizophrenia and aging - two conditions that involve declines in cognitive flexibility and plasticity 1,2 - cells had divested from SNAP: astrocytes, glutamatergic (excitatory) neurons, and GABAergic (inhibitory) neurons all reduced SNAP expression to corresponding degrees. The distinct astrocytic and neuronal components of SNAP both involved genes in which genetic risk factors for schizophrenia were strongly concentrated. SNAP, which varies quantitatively even among healthy persons of similar age, may underlie many aspects of normal human interindividual differences and be an important point of convergence for multiple kinds of pathophysiology.

STAB2: an updated spatio-temporal cell atlas of the human and mouse brain

The brain is constituted of heterogeneous types of neuronal and non-neuronal cells, which are organized into distinct anatomical regions, and show precise regulation of gene expression during development, aging and function. In the current database release, STAB2 provides a systematic cellular map of the human and mouse brain by integrating recently published large-scale single-cell and single-nucleus RNA-sequencing datasets from diverse regions and across lifespan. We applied a hierarchical strategy of unsupervised clustering on the integrated single-cell transcriptomic datasets to precisely annotate the cell types and subtypes in the human and mouse brain. Currently, STAB2 includes 71 and 61 different cell subtypes defined in the human and mouse brain, respectively. It covers 63 subregions and 15 developmental stages of human brain, and 38 subregions and 30 developmental stages of mouse brain, generating a comprehensive atlas for exploring spatiotemporal transcriptomic dynamics in the mammalian brain. We also augmented web interfaces for querying and visualizing the gene expression in specific cell types. STAB2 is freely available at https://mai.fudan.edu.cn/stab2.

Isolated loss of the AUTS2 long isoform, brain-wide or targeted to Calbindin-lineage cells, generates a specific suite of brain, behavioral, and molecular pathologies

Rearrangements within the AUTS2 region are associated with a rare syndromic disorder with intellectual disability, developmental delay, and behavioral abnormalities as core features. In addition, smaller regional variants are linked to wide range of neuropsychiatric disorders, underscoring the gene's essential role in brain development. Like many essential neurodevelopmental genes, AUTS2 is large and complex, generating distinct long (AUTS2-l) and short (AUTS2-s) protein isoforms from alternative promoters. Although evidence suggests unique isoform functions, the contributions of each isoform to specific AUTS2-linked phenotypes have not been clearly resolved. Furthermore, Auts2 is widely expressed across the developing brain, but cell populations most central to disease presentation have not been determined. In this study, we focused on the specific roles of AUTS2-l in brain development, behavior, and postnatal brain gene expression, showing that brain-wide AUTS2-l ablation leads to specific subsets of the recessive pathologies associated with mutations in 3' exons (exons 8-19) that disrupt both major isoforms. We identify downstream genes that could explain expressed phenotypes including hundreds of putative direct AUTS2-l target genes. Furthermore, in contrast to 3' Auts2 mutations which lead to dominant hypoactivity, AUTS2-l loss-of-function is associated with dominant hyperactivity and repetitive behaviors, phenotypes exhibited by many human patients. Finally, we show that AUTS2-l ablation in Calbindin 1-expressing cell lineages is sufficient to yield learning/memory deficits and hyperactivity with abnormal dentate gyrus granule cell maturation, but not other phenotypic effects. These data provide new clues to in vivo AUTS2-l functions and novel information relevant to genotype-phenotype correlations in the human AUTS2 region.

Uniform quantification of single-nucleus ATAC-seq data with Paired-Insertion Counting (PIC) and a model-based insertion rate estimator

Existing approaches to scoring single-nucleus assay for transposase-accessible chromatin with sequencing (snATAC-seq) feature matrices from sequencing reads are inconsistent, affecting downstream analyses and displaying artifacts. We show that, even with sparse single-cell data, quantitative counts are informative for estimating the regulatory state of a cell, which calls for a consistent treatment. We propose Paired-Insertion Counting as a uniform method for snATAC-seq feature characterization and provide a probability model for inferring latent insertion dynamics from snATAC-seq count matrices.

Cellular development and evolution of the mammalian cerebellum

The expansion of the neocortex, a hallmark of mammalian evolution1,2, was accompanied by an increase in cerebellar neuron numbers3. However, little is known about the evolution of the cellular programmes underlying the development of the cerebellum in mammals. In this study we generated single-nucleus RNA-sequencing data for around 400,000 cells to trace the development of the cerebellum from early neurogenesis to adulthood in human, mouse and the marsupial opossum. We established a consensus classification of the cellular diversity in the developing mammalian cerebellum and validated it by spatial mapping in the fetal human cerebellum. Our cross-species analyses revealed largely conserved developmental dynamics of cell-type generation, except for Purkinje cells, for which we observed an expansion of early-born subtypes in the human lineage. Global transcriptome profiles, conserved cell-state markers and gene-expression trajectories across neuronal differentiation show that cerebellar cell-type-defining programmes have been overall preserved for at least 160 million years. However, we also identified many orthologous genes that gained or lost expression in cerebellar neural cell types in one of the species or evolved new expression trajectories during neuronal differentiation, indicating widespread gene repurposing at the cell-type level. In sum, our study unveils shared and lineage-specific gene-expression programmes governing the development of cerebellar cells and expands our understanding of mammalian brain evolution.

Human neuronal maturation comes of age: cellular mechanisms and species differences

The delayed and prolonged postmitotic maturation of human neurons, compared with neurons from other species, may contribute to human-specific cognitive abilities and neurological disorders. Here we review the mechanisms of neuronal maturation, applying lessons from model systems to understand the specific features of protracted human cortical maturation and species differences. We cover cell-intrinsic features of neuronal maturation, including transcriptional, epigenetic and metabolic mechanisms, as well as cell-extrinsic features, including the roles of activity and synapses, the actions of glial cells and the contribution of the extracellular matrix. We discuss evidence for species differences in biochemical reaction rates, the proposed existence of an epigenetic maturation clock and the contributions of both general and modular mechanisms to species-specific maturation timing. Finally, we suggest approaches to measure, improve and accelerate the maturation of human neurons in culture, examine crosstalk and interactions among these different aspects of maturation and propose conceptual models to guide future studies.

Slide-tags enables single-nucleus barcoding for multimodal spatial genomics

Recent technological innovations have enabled the high-throughput quantification of gene expression and epigenetic regulation within individual cells, transforming our understanding of how complex tissues are constructed1-6. However, missing from these measurements is the ability to routinely and easily spatially localize these profiled cells. We developed a strategy, Slide-tags, in which single nuclei within an intact tissue section are tagged with spatial barcode oligonucleotides derived from DNA-barcoded beads with known positions. These tagged nuclei can then be used as an input into a wide variety of single-nucleus profiling assays. Application of Slide-tags to the mouse hippocampus positioned nuclei at less than 10 μm spatial resolution and delivered whole-transcriptome data that are indistinguishable in quality from ordinary single-nucleus RNA-sequencing data. To demonstrate that Slide-tags can be applied to a wide variety of human tissues, we performed the assay on brain, tonsil and melanoma. We revealed cell-type-specific spatially varying gene expression across cortical layers and spatially contextualized receptor-ligand interactions driving B cell maturation in lymphoid tissue. A major benefit of Slide-tags is that it is easily adaptable to almost any single-cell measurement technology. As a proof of principle, we performed multiomic measurements of open chromatin, RNA and T cell receptor (TCR) sequences in the same cells from metastatic melanoma, identifying transcription factor motifs driving cancer cell state transitions in spatially distinct microenvironments. Slide-tags offers a universal platform for importing the compendium of established single-cell measurements into the spatial genomics repertoire.

ALDH5A1-deficient iPSC-derived excitatory and inhibitory neurons display cell type specific alterations

Succinic semialdehyde dehydrogenase deficiency (SSADHD) is a neurometabolic disorder caused by ALDH5A1 mutations presenting with autism and epilepsy. SSADHD leads to impaired GABA metabolism and results in accumulation of GABA and γ-hydroxybutyrate (GHB), which alter neurotransmission and are thought to lead to neurobehavioral symptoms. However, why increased inhibitory neurotransmitters lead to seizures remains unclear. We used induced pluripotent stem cells from SSADHD patients (one female and two male) and differentiated them into GABAergic and glutamatergic neurons. SSADHD iGABA neurons show altered GABA metabolism and concomitant changes in expression of genes associated with inhibitory neurotransmission. In contrast, glutamatergic neurons display increased spontaneous activity and upregulation of mitochondrial genes. CRISPR correction of the pathogenic variants or SSADHD mRNA expression rescue various metabolic and functional abnormalities in human neurons. Our findings uncover a previously unknown role for SSADHD in excitatory human neurons and provide unique insights into the cellular and molecular basis of SSADHD and potential therapeutic interventions.

LRRC37B is a human modifier of voltage-gated sodium channels and axon excitability in cortical neurons

The enhanced cognitive abilities characterizing the human species result from specialized features of neurons and circuits. Here, we report that the hominid-specific gene LRRC37B encodes a receptor expressed in human cortical pyramidal neurons (CPNs) and selectively localized to the axon initial segment (AIS), the subcellular compartment triggering action potentials. Ectopic expression of LRRC37B in mouse CPNs in vivo leads to reduced intrinsic excitability, a distinctive feature of some classes of human CPNs. Molecularly, LRRC37B binds to the secreted ligand FGF13A and to the voltage-gated sodium channel (Nav) β-subunit SCN1B. LRRC37B concentrates inhibitory effects of FGF13A on Nav channel function, thereby reducing excitability, specifically at the AIS level. Electrophysiological recordings in adult human cortical slices reveal lower neuronal excitability in human CPNs expressing LRRC37B. LRRC37B thus acts as a species-specific modifier of human neuron excitability, linking human genome and cell evolution, with important implications for human brain function and diseases.

Patch-seq: Advances and Biological Applications

Multimodal analysis of gene-expression patterns, electrophysiological properties, and morphological phenotypes at the single-cell/single-nucleus level has been arduous because of the diversity and complexity of neurons. The emergence of Patch-sequencing (Patch-seq) directly links transcriptomics, morphology, and electrophysiology, taking neuroscience research to a multimodal era. In this review, we summarized the development of Patch-seq and recent applications in the cortex, hippocampus, and other nervous systems. Through generating multimodal cell type atlases, targeting specific cell populations, and correlating transcriptomic data with phenotypic information, Patch-seq has provided new insight into outstanding questions in neuroscience. We highlight the challenges and opportunities of Patch-seq in neuroscience and hope to shed new light on future neuroscience research.

Human disease-specific cell signatures in non-lesional tissue in Multiple Sclerosis detected by single-cell and spatial transcriptomics

Recent investigations of cell type changes in Multiple Sclerosis (MS) using single-cell profiling methods have focused on active lesional and peri-lesional brain tissue, and have implicated a number of peripheral and central nervous system cell types. However, an important question is the extent to which so-called "normal-appearing" non-lesional tissue in individuals with MS accumulate changes over the lifespan. Here, we compared post-mortem non-lesional brain tissue from donors with a pathological or clinical diagnosis of MS from the Religious Orders Study or Rush Memory and Aging Project (ROSMAP) cohorts to age- and sex-matched brains from persons without MS (controls). We profiled three brain regions using single-nucleus RNA-seq: dorsolateral prefrontal cortex (DLPFC), normal appearing white matter (NAWM) and the pulvinar in thalamus (PULV), from 15 control individuals, 8 individuals with MS, and 5 individuals with other detrimental pathologies accompanied by demyelination, resulting in a total of 78 samples. We identified region- and cell type-specific differences in non-lesional samples from individuals diagnosed with MS and/or exhibiting secondary demyelination with other neurological conditions, as compared to control donors. These differences included lower proportions of oligodendrocytes with expression of myelination related genes MOBP, MBP, PLP1, as well as higher proportions of CRYAB+ oligodendrocytes in all three brain regions. Among microglial signatures, we identified subgroups that were higher in both demyelination (TMEM163+/ERC2+), as well as those that were specifically higher in MS donors (HIF1A+/SPP1+) and specifically in donors with secondary demyelination (SOCS6+/MYO1E+), in both white and grey matter. To validate our findings, we generated Visium spatial transcriptomics data on matched tissue from 13 donors, and recapitulated our observations of gene expression differences in oligodendrocytes and microglia. Finally, we show that some of the differences observed between control and MS donors in NAWM mirror those previously reported between control WM and active lesions in MS donors. Overall, our investigation sheds additional light on cell type- and disease-specific differences present even in non-lesional white and grey matter tissue, highlighting widespread cellular signatures that may be associated with downstream pathological changes.

Cell type-specific enrichment of somatic aneuploidy in the mammalian brain

Somatic mutations alter the genomes of a subset of an individual's brain cells1-3, impacting gene regulation and contributing to disease processes4,5. Mosaic single nucleotide variants have been characterized with single-cell resolution in the brain2,3, but we have limited information about large-scale structural variation, including whole-chromosome duplication or loss1,6,7. We used a dataset of over 415,000 single-cell DNA methylation and chromatin conformation profiles across the adult mouse brain to identify aneuploid cells comprehensively. Whole-chromosome loss or duplication occurred in <1% of cells, with rates up to 1.8% in non-neuronal cell types, including oligodendrocyte precursors and pericytes. Among all aneuploidies, we observed a strong enrichment of trisomy on chromosome 16, which is syntenic with human chromosome 21 and constitutively trisomic in Down syndrome. Chromosome 16 trisomy occurred in multiple cell types and across brain regions, suggesting that nondisjunction is a recurrent feature of somatic variation in the brain.

Shaping the brain: The emergence of cortical structure and folding

The cerebral cortex-the brain's covering and largest region-has increased in size and complexity in humans and supports higher cognitive functions such as language and abstract thinking. There is a growing understanding of the human cerebral cortex, including the diversity and number of cell types that it contains, as well as of the developmental mechanisms that shape cortical structure and organization. In this review, we discuss recent progress in our understanding of molecular and cellular processes, as well as mechanical forces, that regulate the folding of the cerebral cortex. Advances in human genetics, coupled with experimental modeling in gyrencephalic species, have provided insights into the central role of cortical progenitors in the gyrification and evolutionary expansion of the cerebral cortex. These studies are essential for understanding the emergence of structural and functional organization during cortical development and the pathogenesis of neurodevelopmental disorders associated with cortical malformations.

Auditory brainstem responses in the nine-banded armadillo (Dasypus novemcinctus)

The auditory brainstem response (ABR) to tone burst stimuli of thirteen frequencies ranging from 0.5 to 48 kHz was recorded in the nine-banded armadillo (Dasypus novemcinctus), the only extant member of the placental mammal superorder Xenarthra in North America. The armadillo ABR consisted of five main peaks that were visible within the first 10 ms when stimuli were presented at high intensities. The latency of peak I of the armadillo ABR increased as stimulus intensity decreased by an average of 20 μs/dB. Estimated frequency-specific thresholds identified by the ABR were used to construct an estimate of the armadillo audiogram describing the mean thresholds of the eight animals tested. The majority of animals tested (six out of eight) exhibited clear responses to stimuli from 0.5 to 38 kHz, and two animals exhibited responses to stimuli of 48 kHz. Across all cases, the lowest thresholds were observed for frequencies from 8 to 12 kHz. Overall, we observed that the armadillo estimated audiogram bears a similar pattern as those observed using ABR in members of other mammalian clades, including marsupials and later-derived placental mammals.

Evolutionary conservation of hippocampal mossy fiber synapse properties

Various specialized structural/functional properties are considered essential for contextual memory encoding by hippocampal mossy fiber (MF) synapses. Although investigated to exquisite detail in model organisms, synapses, including MFs, have undergone minimal functional interrogation in humans. To determine the translational relevance of rodent findings, we evaluated MF properties within human tissue resected to treat epilepsy. Human MFs exhibit remarkably similar hallmark features to rodents, including AMPA receptor-dominated synapses with small contributions from NMDA and kainate receptors, large dynamic range with strong frequency facilitation, NMDA receptor-independent presynaptic long-term potentiation, and strong cyclic AMP (cAMP) sensitivity of release. Array tomography confirmed the evolutionary conservation of MF ultrastructure. The astonishing congruence of rodent and human MF core features argues that the basic MF properties delineated in animal models remain critical to human MF function. Finally, a selective deficit in GABAergic inhibitory tone onto human MF postsynaptic targets suggests that unrestrained detonator excitatory drive contributes to epileptic circuit hyperexcitability.

Enhanced production of mesencephalic dopaminergic neurons from lineage-restricted human undifferentiated stem cells

Current differentiation protocols for generating mesencephalic dopaminergic (mesDA) neurons from human pluripotent stem cells result in grafts containing only a small proportion of mesDA neurons when transplanted in vivo. In this study, we develop lineage-restricted undifferentiated stem cells (LR-USCs) from pluripotent stem cells, which enhances their potential for differentiating into caudal midbrain floor plate progenitors and mesDA neurons. Using a ventral midbrain protocol, 69% of LR-USCs become bona fide caudal midbrain floor plate progenitors, compared to only 25% of human embryonic stem cells (hESCs). Importantly, LR-USCs generate significantly more mesDA neurons under midbrain and hindbrain conditions in vitro and in vivo. We demonstrate that midbrain-patterned LR-USC progenitors transplanted into 6-hydroxydopamine-lesioned rats restore function in a clinically relevant non-pharmacological behavioral test, whereas midbrain-patterned hESC-derived progenitors do not. This strategy demonstrates how lineage restriction can prevent the development of undesirable lineages and enhance the conditions necessary for mesDA neuron generation.

Gene expression in the primate orbitofrontal cortex related to anxious temperament

Anxiety disorders are among the most prevalent psychiatric disorders, causing significant suffering and disability. Relative to other psychiatric disorders, anxiety disorders tend to emerge early in life, supporting the importance of developmental mechanisms in their emergence and maintenance. Behavioral inhibition (BI) is a temperament that emerges early in life and, when stable and extreme, is linked to an increased risk for the later development of anxiety disorders and other stress-related psychopathology. Understanding the neural systems and molecular mechanisms underlying this dispositional risk could provide insight into treatment targets for anxiety disorders. Nonhuman primates (NHPs) have an anxiety-related temperament, called anxious temperament (AT), that is remarkably similar to BI in humans, facilitating the design of highly translational models for studying the early risk for stress-related psychopathology. Because of the recent evolutionary divergence between humans and NHPs, many of the anxiety-related brain regions that contribute to psychopathology are highly similar in terms of their structure and function, particularly with respect to the prefrontal cortex. The orbitofrontal cortex plays a critical role in the flexible encoding and regulation of threat responses, in part through connections with subcortical structures like the amygdala. Here, we explore individual differences in the transcriptional profile of cells within the region, using laser capture microdissection and single nuclear sequencing, providing insight into the molecules underlying individual differences in AT-related function of the pOFC, with a particular focus on previously implicated cellular systems, including neurotrophins and glucocorticoid signaling.

Cell-type-specific 3D-genome organization and transcription regulation in the brain

3D organization of the genome plays a critical role in regulating gene expression. However, it remains unclear how chromatin organization differs among different cell types in the brain. Here we used genome-scale DNA and RNA imaging to investigate 3D-genome organization in transcriptionally distinct cell types in the primary motor cortex of the mouse brain. We uncovered a wide spectrum of differences in the nuclear architecture and 3D-genome organization among different cell types, ranging from the physical size of the cell nucleus to the active-inactive chromatin compartmentalization and radial positioning of chromatin loci within the nucleus. These cell-type-dependent variations in nuclear architecture and chromatin organization exhibited strong correlation with both total transcriptional activity of the cell and transcriptional regulation of cell-type-specific marker genes. Moreover, we found that the methylated-DNA-binding protein MeCP2 regulates transcription in a divergent manner, depending on the nuclear radial positions of chromatin loci, through modulating active-inactive chromatin compartmentalization.

Single-cell transcriptomics reveals the brain evolution of web-building spiders

Spiders are renowned for their efficient capture of flying insects using intricate aerial webs. How the spider nervous systems evolved to cope with this specialized hunting strategy and various environmental clues in an aerial space remains unknown. Here we report a brain-cell atlas of >30,000 single-cell transcriptomes from a web-building spider (Hylyphantes graminicola). Our analysis revealed the preservation of ancestral neuron types in spiders, including the potential coexistence of noradrenergic and octopaminergic neurons, and many peptidergic neuronal types that are lost in insects. By comparing the genome of two newly sequenced plesiomorphic burrowing spiders with three aerial web-building spiders, we found that the positively selected genes in the ancestral branch of web-building spiders were preferentially expressed (42%) in the brain, especially in the three mushroom body-like neuronal types. By gene enrichment analysis and RNAi experiments, these genes were suggested to be involved in the learning and memory pathway and may influence the spiders' web-building and hunting behaviour. Our results provide key sources for understanding the evolution of behaviour in spiders and reveal how molecular evolution drives neuron innovation and the diversification of associated complex behaviours.

Conserved and divergent gene regulatory programs of the mammalian neocortex

Divergence of cis-regulatory elements drives species-specific traits1, but how this manifests in the evolution of the neocortex at the molecular and cellular level remains unclear. Here we investigated the gene regulatory programs in the primary motor cortex of human, macaque, marmoset and mouse using single-cell multiomics assays, generating gene expression, chromatin accessibility, DNA methylome and chromosomal conformation profiles from a total of over 200,000 cells. From these data, we show evidence that divergence of transcription factor expression corresponds to species-specific epigenome landscapes. We find that conserved and divergent gene regulatory features are reflected in the evolution of the three-dimensional genome. Transposable elements contribute to nearly 80% of the human-specific candidate cis-regulatory elements in cortical cells. Through machine learning, we develop sequence-based predictors of candidate cis-regulatory elements in different species and demonstrate that the genomic regulatory syntax is highly preserved from rodents to primates. Finally, we show that epigenetic conservation combined with sequence similarity helps to uncover functional cis-regulatory elements and enhances our ability to interpret genetic variants contributing to neurological disease and traits.

Cell type specializations of the vocal-motor cortex in songbirds

Identifying molecular specializations in cortical circuitry supporting complex behaviors, like learned vocalizations, requires understanding of the neuroanatomical context from which these circuits arise. In songbirds, the robust arcopallial nucleus (RA) provides descending cortical projections for fine vocal-motor control. Using single-nuclei transcriptomics and spatial gene expression mapping in zebra finches, we have defined cell types and molecular specializations that distinguish RA from adjacent regions involved in non-vocal motor and sensory processing. We describe an RA-specific projection neuron, differential inhibitory subtypes, and glia specializations and have probed predicted GABAergic interneuron subtypes electrophysiologically within RA. Several cell-specific markers arise developmentally in a sex-dependent manner. Our interactive apps integrate cellular data with developmental and spatial distribution data from the gene expression brain atlas ZEBrA. Users can explore molecular specializations of vocal-motor neurons and support cells that likely reflect adaptations key to the physiology and evolution of vocal control circuits and refined motor skills.

Connecting single-cell transcriptomes to projectomes in mouse visual cortex

The mammalian brain is composed of diverse neuron types that play different functional roles. Recent single-cell RNA sequencing approaches have led to a whole brain taxonomy of transcriptomically-defined cell types, yet cell type definitions that include multiple cellular properties can offer additional insights into a neuron's role in brain circuits. While the Patch-seq method can investigate how transcriptomic properties relate to the local morphological and electrophysiological properties of cell types, linking transcriptomic identities to long-range projections is a major unresolved challenge. To address this, we collected coordinated Patch-seq and whole brain morphology data sets of excitatory neurons in mouse visual cortex. From the Patch-seq data, we defined 16 integrated morpho-electric-transcriptomic (MET)-types; in parallel, we reconstructed the complete morphologies of 300 neurons. We unified the two data sets with a multi-step classifier, to integrate cell type assignments and interrogate cross-modality relationships. We find that transcriptomic variations within and across MET-types correspond with morphological and electrophysiological phenotypes. In addition, this variation, along with the anatomical location of the cell, can be used to predict the projection targets of individual neurons. We also shed new light on infragranular cell types and circuits, including cell-type-specific, interhemispheric projections. With this approach, we establish a comprehensive, integrated taxonomy of excitatory neuron types in mouse visual cortex and create a system for integrated, high-dimensional cell type classification that can be extended to the whole brain and potentially across species.

Does decreased autophagy and dysregulation of LC3A in astrocytes play a role in major depressive disorder?

Astrocytic dysfunction contributes to the molecular pathogenesis of major depressive disorder (MDD). However, the astrocytic subtype that mainly contributes to MDD etiology and whether dysregulated autophagy in astrocytes is associated with MDD remain unknown. Using a single-nucleus RNA sequencing (snRNA-seq) atlas, three astrocyte subtypes were identified in MDD, while C2 State-1Q astrocytes showed aberrant changes in both cell proportion and most differentially expressed genes compared with other subtypes. Moreover, autophagy pathways were commonly inhibited in astrocytes in the prefrontal cortices (PFCs) of patients with MDD, especially in C2 State-1Q astrocytes. Furthermore, by integrating snRNA-seq and bulk transcriptomic data, we found significant reductions in LC3A expression levels in the PFC region of CUMS-induced depressed mice, as well as in postmortem PFC tissues and peripheral blood samples from patients with MDD. These results were further validated by qPCR using whole-blood samples from patients with MDD and healthy controls. Finally, LC3A expression in the whole blood of patients with MDD was negatively associated with the severity of depressive symptoms. Overall, our results underscore autophagy inhibition in PFC astrocytes as a common molecular characteristic in MDD and might reveal a novel potential diagnostic marker LC3A.

CellGO: a novel deep learning-based framework and webserver for cell-type-specific gene function interpretation

Interpreting the function of genes and gene sets identified from omics experiments remains a challenge, as current pathway analysis tools often fail to consider the critical biological context, such as tissue or cell-type specificity. To address this limitation, we introduced CellGO. CellGO tackles this challenge by leveraging the visible neural network (VNN) and single-cell gene expressions to mimic cell-type-specific signaling propagation along the Gene Ontology tree within a cell. This design enables a novel scoring system to calculate the cell-type-specific gene-pathway paired active scores, based on which, CellGO is able to identify cell-type-specific active pathways associated with single genes. In addition, by aggregating the activities of single genes, CellGO extends its capability to identify cell-type-specific active pathways for a given gene set. To enhance biological interpretation, CellGO offers additional features, including the identification of significantly active cell types and driver genes and community analysis of pathways. To validate its performance, CellGO was assessed using a gene set comprising mixed cell-type markers, confirming its ability to discern active pathways across distinct cell types. Subsequent benchmarking analyses demonstrated CellGO's superiority in effectively identifying cell types and their corresponding cell-type-specific pathways affected by gene knockouts, using either single genes or sets of genes differentially expressed between knockout and control samples. Moreover, CellGO demonstrated its ability to infer cell-type-specific pathogenesis for disease risk genes. Accessible as a Python package, CellGO also provides a user-friendly web interface, making it a versatile and accessible tool for researchers in the field.

Single-cell genomics reveals region-specific developmental trajectories underlying neuronal diversity in the human hypothalamus

The development and diversity of neuronal subtypes in the human hypothalamus has been insufficiently characterized. To address this, we integrated transcriptomic data from 241,096 cells (126,840 newly generated) in the prenatal and adult human hypothalamus to reveal a temporal trajectory from proliferative stem cell populations to mature hypothalamic cell types. Iterative clustering of the adult neurons identified 108 robust transcriptionally distinct neuronal subtypes representing 10 hypothalamic nuclei. Pseudotime trajectories provided insights into the genes driving formation of these nuclei. Comparisons to single-cell transcriptomic data from the mouse hypothalamus suggested extensive conservation of neuronal subtypes despite certain differences in species-enriched gene expression. The uniqueness of hypothalamic neuronal lineages was examined developmentally by comparing excitatory lineages present in cortex and inhibitory lineages in ganglionic eminence, revealing both distinct and shared drivers of neuronal maturation across the human forebrain. These results provide a comprehensive transcriptomic view of human hypothalamus development through gestation and adulthood at cellular resolution.

Cell lineage analysis with somatic mutations reveals late divergence of neuronal cell types and cortical areas in human cerebral cortex

The mammalian cerebral cortex shows functional specialization into regions with distinct neuronal compositions, most strikingly in the human brain, but little is known in about how cellular lineages shape cortical regional variation and neuronal cell types during development. Here, we use somatic single nucleotide variants (sSNVs) to map lineages of neuronal sub-types and cortical regions. Early-occurring sSNVs rarely respect Brodmann area (BA) borders, while late-occurring sSNVs mark neuron-generating clones with modest regional restriction, though descendants often dispersed into neighboring BAs. Nevertheless, in visual cortex, BA17 contains 30-70% more sSNVs compared to the neighboring BA18, with clones across the BA17/18 border distributed asymmetrically and thus displaying different cortex-wide dispersion patterns. Moreover, we find that excitatory neuron-generating clones with modest regional restriction consistently share low-mosaic sSNVs with some inhibitory neurons, suggesting significant co-generation of excitatory and some inhibitory neurons in the dorsal cortex. Our analysis reveals human-specific cortical cell lineage patterns, with both regional inhomogeneities in progenitor proliferation and late divergence of excitatory/inhibitory lineages.

Neurogenesis in primates versus rodents and the value of non-human primate models

Neurogenesis, the process of generating neurons from neural stem cells, occurs during both embryonic and adult stages, with each stage possessing distinct characteristics. Dysfunction in either stage can disrupt normal neural development, impair cognitive functions, and lead to various neurological disorders. Recent technological advancements in single-cell multiomics and gene-editing have facilitated investigations into primate neurogenesis. Here, we provide a comprehensive overview of neurogenesis across rodents, non-human primates, and humans, covering embryonic development to adulthood and focusing on the conservation and diversity among species. While non-human primates, especially monkeys, serve as valuable models with closer neural resemblance to humans, we highlight the potential impacts and limitations of non-human primate models on both physiological and pathological neurogenesis research.

Functional genomics and systems biology in human neuroscience

Neuroscience research has entered a phase of key discoveries in the realm of neurogenomics owing to strong financial and intellectual support for resource building and tool development. The previous challenge of tissue heterogeneity has been met with the application of techniques that can profile individual cells at scale. Moreover, the ability to perturb genes, gene regulatory elements and neuronal activity in a cell-type-specific manner has been integrated with gene expression studies to uncover the functional underpinnings of the genome at a systems level. Although these insights have necessarily been grounded in model systems, we now have the opportunity to apply these approaches in humans and in human tissue, thanks to advances in human genetics, brain imaging and tissue collection. We acknowledge that there will probably always be limits to the extent to which we can apply the genomic tools developed in model systems to human neuroscience; however, as we describe in this Perspective, the neuroscience field is now primed with an optimal foundation for tackling this ambitious challenge. The application of systems-level network analyses to these datasets will facilitate a deeper appreciation of human neurogenomics that cannot otherwise be achieved from directly observable phenomena.

Population-level integration of single-cell datasets enables multi-scale analysis across samples

The increasing generation of population-level single-cell atlases has the potential to link sample metadata with cellular data. Constructing such references requires integration of heterogeneous cohorts with varying metadata. Here we present single-cell population level integration (scPoli), an open-world learner that incorporates generative models to learn sample and cell representations for data integration, label transfer and reference mapping. We applied scPoli on population-level atlases of lung and peripheral blood mononuclear cells, the latter consisting of 7.8 million cells across 2,375 samples. We demonstrate that scPoli can explain sample-level biological and technical variations using sample embeddings revealing genes associated with batch effects and biological effects. scPoli is further applicable to single-cell sequencing assay for transposase-accessible chromatin and cross-species datasets, offering insights into chromatin accessibility and comparative genomics. We envision scPoli becoming an important tool for population-level single-cell data integration facilitating atlas use but also interpretation by means of multi-scale analyses.

Comparative single-cell transcriptomic analysis of primate brains highlights human-specific regulatory evolution

Enhanced cognitive function in humans is hypothesized to result from cortical expansion and increased cellular diversity. However, the mechanisms that drive these phenotypic innovations remain poorly understood, in part because of the lack of high-quality cellular resolution data in human and non-human primates. Here, we take advantage of single-cell expression data from the middle temporal gyrus of five primates (human, chimp, gorilla, macaque and marmoset) to identify 57 homologous cell types and generate cell type-specific gene co-expression networks for comparative analysis. Although orthologue expression patterns are generally well conserved, we find 24% of genes with extensive differences between human and non-human primates (3,383 out of 14,131), which are also associated with multiple brain disorders. To assess the functional significance of gene expression differences in an evolutionary context, we evaluate changes in network connectivity across meta-analytic co-expression networks from 19 animals. We find that a subset of these genes has deeply conserved co-expression across all non-human animals, and strongly divergent co-expression relationships in humans (139 out of 3,383, <1% of primate orthologues). Genes with human-specific cellular expression and co-expression profiles (such as NHEJ1, GTF2H2, C2 and BBS5) typically evolve under relaxed selective constraints and may drive rapid evolutionary change in brain function.

Cortical circuit dynamics underlying motor skill learning: from rodents to humans

Motor learning is crucial for the survival of many animals. Acquiring a new motor skill involves complex alterations in both local neural circuits in many brain regions and long-range connections between them. Such changes can be observed anatomically and functionally. The primary motor cortex (M1) integrates information from diverse brain regions and plays a pivotal role in the acquisition and refinement of new motor skills. In this review, we discuss how motor learning affects the M1 at synaptic, cellular, and circuit levels. Wherever applicable, we attempt to relate and compare findings in humans, non-human primates, and rodents. Understanding the underlying principles shared by different species will deepen our understanding of the neurobiological and computational basis of motor learning.

Cross-species imputation and comparison of single-cell transcriptomic profiles

Cross-species comparison and prediction of gene expression profiles are important to understand regulatory changes during evolution and to transfer knowledge learned from model organisms to humans. Single-cell RNA-seq (scRNA-seq) profiles enable us to capture gene expression profiles with respect to variations among individual cells; however, cross-species comparison of scRNA-seq profiles is challenging because of data sparsity, batch effects, and the lack of one-to-one cell matching across species. Moreover, single-cell profiles are challenging to obtain in certain biological contexts, limiting the scope of hypothesis generation. Here we developed Icebear, a neural network framework that decomposes single-cell measurements into factors representing cell identity, species, and batch factors. Icebear enables accurate prediction of single-cell gene expression profiles across species, thereby providing high-resolution cell type and disease profiles in under-characterized contexts. Icebear also facilitates direct cross-species comparison of single-cell expression profiles for conserved genes that are located on the X chromosome in eutherian mammals but on autosomes in chicken. This comparison, for the first time, revealed evolutionary and diverse adaptations of X-chromosome upregulation in mammals.

Discrete hippocampal projections are differentially regulated by parvalbumin and somatostatin interneurons

People with schizophrenia show hyperactivity in the ventral hippocampus (vHipp) and we have previously demonstrated distinct behavioral roles for vHipp cell populations. Here, we test the hypothesis that parvalbumin (PV) and somatostatin (SST) interneurons differentially innervate and regulate hippocampal pyramidal neurons based on their projection target. First, we use eGRASP to show that PV-positive interneurons form a similar number of synaptic connections with pyramidal cells regardless of their projection target while SST-positive interneurons preferentially target nucleus accumbens (NAc) projections. To determine if these anatomical differences result in functional changes, we used in vivo opto-electrophysiology to show that SST cells also preferentially regulate the activity of NAc-projecting cells. These results suggest vHipp interneurons differentially regulate that vHipp neurons that project to the medial prefrontal cortex (mPFC) and NAc. Characterization of these cell populations may provide potential molecular targets for the treatment schizophrenia and other psychiatric disorders associated with vHipp dysfunction.