Cell type hierarchy reconstruction via reconciliation of multi-resolution cluster tree

Transient expression of the neuropeptide galanin modulates peripheral‑to‑central connectivity in the somatosensory thalamus during whisker development in mice

The significance of transient neuropeptide expression during postnatal brain development is unknown. Here, we show that galanin expression in the ventrobasal thalamus of infant mice coincides with whisker map development and modulates subcortical circuit wiring. Time-resolved neuroanatomy and single-nucleus RNA-seq identified complementary galanin (Gal) and galanin receptor 1 (Galr1) expression in the ventrobasal thalamus and the principal sensory nucleus of the trigeminal nerve (Pr5), respectively. Somatodendritic galanin release from the ventrobasal thalamus was time-locked to the first postnatal week, when Gal1R+ Pr5 afferents form glutamatergic (Slc17a6+) synapses for the topographical whisker map to emerge. RNAi-mediated silencing of galanin expression disrupted glutamatergic synaptogenesis, which manifested as impaired whisker-dependent exploratory behaviors in infant mice, with behavioral abnormalities enduring into adulthood. Pharmacological probing of receptor selectivity in vivo corroborated that target recognition and synaptogenesis in the thalamus, at least in part, are reliant on agonist-induced Gal1R activation in inbound excitatory axons. Overall, we suggest a neuropeptide-dependent developmental mechanism to contribute to the topographical specification of a fundamental sensory neurocircuit in mice.

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.

CHOIR improves significance-based detection of cell types and states from single-cell data

Clustering is a critical step in the analysis of single-cell data, as it enables the discovery and characterization of putative cell types and states. However, most popular clustering tools do not subject clustering results to statistical inference testing, leading to risks of overclustering or underclustering data and often resulting in ineffective identification of cell types with widely differing prevalence. To address these challenges, we present CHOIR (clustering hierarchy optimization by iterative random forests), which applies a framework of random forest classifiers and permutation tests across a hierarchical clustering tree to statistically determine which clusters represent distinct populations. We demonstrate the enhanced performance of CHOIR through extensive benchmarking against 14 existing clustering methods across 100 simulated and 4 real single-cell RNA-seq, ATAC-seq, spatial transcriptomic, and multi-omic datasets. CHOIR can be applied to any single-cell data type and provides a flexible, scalable, and robust solution to the important challenge of identifying biologically relevant cell groupings within heterogeneous single-cell data.

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.

Self-supervised deep clustering of single-cell RNA-seq data to hierarchically detect rare cell populations

Single-cell RNA sequencing (scRNA-seq) is a widely used technique for characterizing individual cells and studying gene expression at the single-cell level. Clustering plays a vital role in grouping similar cells together for various downstream analyses. However, the high sparsity and dimensionality of large scRNA-seq data pose challenges to clustering performance. Although several deep learning-based clustering algorithms have been proposed, most existing clustering methods have limitations in capturing the precise distribution types of the data or fully utilizing the relationships between cells, leaving a considerable scope for improving the clustering performance, particularly in detecting rare cell populations from large scRNA-seq data. We introduce DeepScena, a novel single-cell hierarchical clustering tool that fully incorporates nonlinear dimension reduction, negative binomial-based convolutional autoencoder for data fitting, and a self-supervision model for cell similarity enhancement. In comprehensive evaluation using multiple large-scale scRNA-seq datasets, DeepScena consistently outperformed seven popular clustering tools in terms of accuracy. Notably, DeepScena exhibits high proficiency in identifying rare cell populations within large datasets that contain large numbers of clusters. When applied to scRNA-seq data of multiple myeloma cells, DeepScena successfully identified not only previously labeled large cell types but also subpopulations in CD14 monocytes, T cells and natural killer cells, respectively.

scAce: an adaptive embedding and clustering method for single-cell gene expression data

MOTIVATION

Since the development of single-cell RNA sequencing (scRNA-seq) technologies, clustering analysis of single-cell gene expression data has been an essential tool for distinguishing cell types and identifying novel cell types. Even though many methods have been available for scRNA-seq clustering analysis, the majority of them are constrained by the requirement on predetermined cluster numbers or the dependence on selected initial cluster assignment.

RESULTS

In this article, we propose an adaptive embedding and clustering method named scAce, which constructs a variational autoencoder to simultaneously learn cell embeddings and cluster assignments. In the scAce method, we develop an adaptive cluster merging approach which achieves improved clustering results without the need to estimate the number of clusters in advance. In addition, scAce provides an option to perform clustering enhancement, which can update and enhance cluster assignments based on previous clustering results from other methods. Based on computational analysis of both simulated and real datasets, we demonstrate that scAce outperforms state-of-the-art clustering methods for scRNA-seq data, and achieves better clustering accuracy and robustness.

AVAILABILITY AND IMPLEMENTATION

The scAce package is implemented in python 3.8 and is freely available from https://github.com/sldyns/scAce.

Multiscale geometric and topological analyses for characterizing and predicting immune responses from single cell data

Single cell genomics has revolutionized our ability to map immune heterogeneity and responses. With the influx of large-scale data sets from diverse modalities, the resolution achieved has supported the long-held notion that immune cells are naturally organized into hierarchical relationships, characterized at multiple levels. Such a multigranular structure corresponds to key geometric and topological features. Given that differences between an effective and ineffective immunological response may not be found at one level, there is vested interest in characterizing and predicting outcomes from such features. In this review, we highlight single cell methods and principles for learning geometric and topological properties of data at multiple scales, discussing their contributions to immunology. Ultimately, multiscale approaches go beyond classical clustering, revealing a more comprehensive picture of cellular heterogeneity.

Accurate estimation of rare cell type fractions from tissue omics data via hierarchical deconvolution

Bulk transcriptomics in tissue samples reflects the average expression levels across different cell types and is highly influenced by cellular fractions. As such, it is critical to estimate cellular fractions to both deconfound differential expression analyses and infer cell type-specific differential expression. Since experimentally counting cells is infeasible in most tissues and studies, in silico cellular deconvolution methods have been developed as an alternative. However, existing methods are designed for tissues consisting of clearly distinguishable cell types and have difficulties estimating highly correlated or rare cell types. To address this challenge, we propose Hierarchical Deconvolution (HiDecon) that uses single-cell RNA sequencing references and a hierarchical cell type tree, which models the similarities among cell types and cell differentiation relationships, to estimate cellular fractions in bulk data. By coordinating cell fractions across layers of the hierarchical tree, cellular fraction information is passed up and down the tree, which helps correct estimation biases by pooling information across related cell types. The flexible hierarchical tree structure also enables estimating rare cell fractions by splitting the tree to higher resolutions. Through simulations and real data applications with the ground truth of measured cellular fractions, we demonstrate that HiDecon significantly outperforms existing methods and accurately estimates cellular fractions.

Identification of a unique subset of tissue-resident memory CD4(+) T cells in Crohn's disease

T cells differentiate into highly diverse subsets and display plasticity depending on the environment. Although lymphocytes are key mediators of inflammation, functional specialization of T cells in inflammatory bowel disease (IBD) has not been effectively described. Here, we performed deep profiling of T cells in the intestinal mucosa of IBD and identified a CD4⁺ tissue-resident memory T cell (Trm) subset that is increased in Crohn's disease (CD) showing unique inflammatory properties. Functionally and transcriptionally distinct CD4⁺ Trm subsets are observed in the inflamed gut mucosa, among which a CD-specific CD4⁺ Trm subset, expressing CD161 and CCR5 along with CD103, displays previously unrecognized pleiotropic signatures of innate and effector activities. These inflammatory features are further enhanced by their spatial proximity to gut epithelial cells. Furthermore, the CD-specific CD4⁺ Trm subset is the most predominant producer of type 1 inflammatory cytokines upon various stimulations among all CD4⁺ T cells, suggesting that the accumulation of this T cell subset is a pathological hallmark of CD. Our results provide comprehensive insights into the pathogenesis of IBD, paving the way for decoding of the molecular mechanisms underlying this disease.

Rare coding variation provides insight into the genetic architecture and phenotypic context of autism

Some individuals with autism spectrum disorder (ASD) carry functional mutations rarely observed in the general population. We explored the genes disrupted by these variants from joint analysis of protein-truncating variants (PTVs), missense variants and copy number variants (CNVs) in a cohort of 63,237 individuals. We discovered 72 genes associated with ASD at false discovery rate (FDR) ≤ 0.001 (185 at FDR ≤ 0.05). De novo PTVs, damaging missense variants and CNVs represented 57.5%, 21.1% and 8.44% of association evidence, while CNVs conferred greatest relative risk. Meta-analysis with cohorts ascertained for developmental delay (DD) (n = 91,605) yielded 373 genes associated with ASD/DD at FDR ≤ 0.001 (664 at FDR ≤ 0.05), some of which differed in relative frequency of mutation between ASD and DD cohorts. The DD-associated genes were enriched in transcriptomes of progenitor and immature neuronal cells, whereas genes showing stronger evidence in ASD were more enriched in maturing neurons and overlapped with schizophrenia-associated genes, emphasizing that these neuropsychiatric disorders may share common pathways to risk.