SciBet as a portable and fast single cell type identifier

Role of PRC2 in the stochastic expression of Aire target genes and development of mimetic cells in the thymus

The transcriptional targets of Aire and the mechanisms controlling their expression in medullary thymic epithelial cells (mTECs) need to be clarified to understand Aire's tolerogenic function. By using a multi-omics single-cell approach coupled with deep scRNA-seq, we examined how Aire controls the transcription of a wide variety of genes in a small fraction of Aire-expressing cells. We found that chromatin repression by PRC2 is an important step for Aire to achieve stochastic gene expression. Aire unleashed the silenced chromatin configuration caused by PRC2, thereby increasing the expression of its functional targets. Besides this preconditioning for Aire's gene induction, we demonstrated that PRC2 also controls the composition of mTECs that mimic the developmental trait of peripheral tissues, i.e., mimetic cells. Of note, this action of PRC2 was independent of Aire and it was more apparent than Aire. Thus, our study uncovered the essential role of polycomb complex for Aire-mediated promiscuous gene expression and the development of mimetic cells.

Long-term engraftment of human stem and progenitor cells for large-scale production of functional immune cells in engineered pigs

Existing immunodeficient pig models have demonstrated limited success in supporting robust human haematopoietic engraftment and multilineage differentiation. Here we hypothesize that the targeted deletion of integrin-associated protein (Cd47) in severe combined immunodeficient pigs, with deletions in the X-linked interleukin-2 receptor gamma chain and recombination activating gene 1, would enable long-term haematopoietic engraftment following transplantation with human haematopoietic stem/progenitor cells. In Cd47-deficient pigs, we observed high levels of human haematopoietic chimerism in the blood and spleen, with functional T and B lymphocytes, natural killer and myeloid cells, as well as robust thymopoiesis. Our findings suggest that severe combined immunodeficient pigs with Cd47 deletion represent an improved preclinical model for studying human haematopoiesis, disease mechanisms and therapies, and offer potential as a bioreactor for large-scale production of human immune cells.

Integrating single-cell data with biological variables

Constructing single-cell atlases requires preserving differences attributable to biological variables, such as cell types, tissue origins, and disease states, while eliminating batch effects. However, existing methods are inadequate in explicitly modeling these biological variables. Here, we introduce SIGNAL, a general framework that leverages biological variables to disentangle biological and technical effects, thereby linking these metadata to data integration. SIGNAL employs a variant of principal component analysis to align multiple batches, enabling the integration of 1 million cells in approximately 2 min. SIGNAL, despite its computational simplicity, surpasses state-of-the-art methods across multiple integration scenarios: 1) heterogeneous datasets, 2) cross-species datasets, 3) simulated datasets, 4) integration on low-quality cell annotations, and 5) reference-based integration. Furthermore, we demonstrate that SIGNAL accurately transfers knowledge from reference to query datasets. Notably, we propose a self-adjustment strategy to restore annotated cell labels potentially distorted during integration. Finally, we apply SIGNAL to multiple large-scale atlases, including a human heart cell atlas containing 2.7 million cells, identifying tissue- and developmental stage-specific subtypes, as well as condition-specific cell states. This underscores SIGNAL's exceptional capability in multiscale analysis.

Single-cell transcriptomics reveals extracellular matrix remodeling and collagen dynamics during lactation in sheep mammary gland

The mammary gland is a dynamic organ with diverse cell populations that maintain glandular homeostasis, particularly during lactation. However, the cellular architecture and molecular mechanisms underlying lactational remodeling in the sheep mammary gland remain incompletely understood. Given similarities in mammary stromal structure, sheep serve as a valuable model for studying lactational changes relevant to the human breast, which experiences collagen loss and sagging during lactation. Utilizing single-cell transcriptomics (scRNA-seq), we mapped the sheep mammary gland's cellular landscape at postpartum days 60 and 150, identifying seven major cell types, including six distinct epithelial clusters. These clusters revealed differentiation among luminal progenitors, hormone-sensing, and myoepithelial cells across peak and late lactation stages. Transcriptomic analysis highlighted pivotal roles for epithelial integrity and ECM remodeling, with myoepithelial cells centrally involved in these processes. We observed significant collagen remodeling driven by fibroblast-epithelial crosstalk and ECM reorganization during late lactation. Comparative analysis with human mammary epithelial cells showed conserved basal and myoepithelial cell populations, while luminal cells diverged across species. This study provides insights into lactation biology and ECM remodeling, offering a framework to inform future studies on lactational adaptation and its implications for human health.

An overview of computational methods in single-cell transcriptomic cell type annotation

The rapid accumulation of single-cell RNA sequencing data has provided unprecedented computational resources for cell type annotation, significantly advancing our understanding of cellular heterogeneity. Leveraging gene expression profiles derived from transcriptomic data, researchers can accurately infer cell types, sparking the development of numerous innovative annotation methods. These methods utilize a range of strategies, including marker genes, correlation-based matching, and supervised learning, to classify cell types. In this review, we systematically examine these annotation approaches based on transcriptomics-specific gene expression profiles and provide a comprehensive comparison and categorization of these methods. Furthermore, we focus on the main challenges in the annotation process, especially the long-tail distribution problem arising from data imbalance in rare cell types. We discuss the potential of deep learning techniques to address these issues and enhance model capability in recognizing novel cell types within an open-world framework.

Deconvolution of cell types and states in spatial multiomics utilizing TACIT

Identifying cell types and states remains a time-consuming, error-prone challenge for spatial biology. While deep learning increasingly plays a role, it is difficult to generalize due to variability at the level of cells, neighborhoods, and niches in health and disease. To address this, we develop TACIT, an unsupervised algorithm for cell annotation using predefined signatures that operates without training data. TACIT uses unbiased thresholding to distinguish positive cells from background, focusing on relevant markers to identify ambiguous cells in multiomic assays. Using five datasets (5,000,000 cells; 51 cell types) from three niches (brain, intestine, gland), TACIT outperforms existing unsupervised methods in accuracy and scalability. Integrating TACIT-identified cell types reveals new phenotypes in two inflammatory gland diseases. Finally, using combined spatial transcriptomics and proteomics, we discover under- and overrepresented immune cell types and states in regions of interest, suggesting multimodality is essential for translating spatial biology to clinical applications.

Differentiation success of reprogrammed cells is heterogeneous in vivo and modulated by somatic cell identity memory

Nuclear reprogramming can change cellular fates. Yet, reprogramming efficiency is low, and the resulting cell types are often not functional. Here, we used nuclear transfer to eggs to follow single cells during reprogramming in vivo. We show that the differentiation success of reprogrammed cells varies across cell types and depends on the expression of genes specific to the previous cellular identity. We find subsets of reprogramming-resistant cells that fail to form functional cell types, undergo cell death, or disrupt normal body patterning. Reducing expression levels of genes specific to the cell type of origin leads to better reprogramming and improved differentiation trajectories. Thus, our work demonstrates that failing to reprogram in vivo is cell type specific and emphasizes the necessity of minimizing aberrant transcripts of the previous somatic identity for improving reprogramming.

Mapping Cell Identity from scRNA-seq: A primer on computational methods

Single cell (sc) technologies mark a conceptual and methodological breakthrough in our way to study cells, the base units of life. Thanks to these technological developments, large-scale initiatives are currently ongoing aimed at mapping of all the cell types in the human body, with the ambitious aim to gain a cell-level resolution of physiological development and disease. Since its broad applicability and ease of interpretation scRNA-seq is probably the most common sc-based application. This assay uses high throughput RNA sequencing to capture gene expression profiles at the sc-level. Subsequently, under the assumption that differences in transcriptional programs correspond to distinct cellular identities, ad-hoc computational methods are used to infer cell types from gene expression patterns. A wide array of computational methods were developed for this task. However, depending on the underlying algorithmic approach and associated computational requirements, each method might have a specific range of application, with implications that are not always clear to the end user. Here we will provide a concise overview on state-of-the-art computational methods for cell identity annotation in scRNA-seq, tailored for new users and non-computational scientists. To this end, we classify existing tools in five main categories, and discuss their key strengths, limitations and range of application.

stAI: a deep learning-based model for missing gene imputation and cell-type annotation of spatial transcriptomics

Spatial transcriptomics technology has revolutionized our understanding of cellular systems by capturing RNA transcript levels in their original spatial context. Single-cell spatial transcriptomics (scST) offers single-cell resolution expression level and precise spatial information of RNA transcripts, while it has a limited capacity for simultaneously detecting a wide range of RNA transcripts, hindering its broader applications. Characterizing the whole transcriptome level and comprehensively annotating cell types represent two significant challenges in scST applications. Despite several proposed methods for one or both tasks, their performance remains inadequate. In this work, we introduce stAI, a deep learning-based model designed to address both missing gene imputation and cell-type annotation for scST data. stAI leverages a joint embedding for the scST and the reference scRNA-seq data with two separate encoder-decoder modules. Both the imputation and annotation are performed within the latent space in a supervised manner, utilizing scRNA-seq data to guide the processes. Experiments for datasets generated from diverse platforms with varying numbers of measured genes were conducted and compared with the updated methods. The results demonstrate that stAI can predict the unmeasured genes, especially the marker genes, with much higher accuracy, and annotate the cell types, including those of small size, with high precision.

Consensus representation of multiple cell-cell graphs from gene signaling pathways for cell type annotation

BACKGROUND

Recent advancements in single-cell RNA sequencing have greatly expanded our knowledge of the heterogeneous nature of tissues. However, robust and accurate cell type annotation continues to be a major challenge, hindered by issues such as marker specificity, batch effects, and a lack of comprehensive spatial and interaction data. Traditional annotation methods often fail to adequately address the complexity of cellular interactions and gene regulatory networks.

RESULTS

We proposed scMCGraph, a comprehensive computational framework that integrates gene expression with pathway activity to accurately annotate cell types within diverse scRNA-seq datasets. Initially, our model constructs multiple pathway-specific views using various pathway databases, which reflect both gene expression and pathway activities. These pathway-specific views are then integrated into a consensus graph. The consensus graph is subsequently utilized to reconstruct the multiple pathway views. Our model demonstrated exceptional robustness and accuracy across various analyses, including cross-platform, cross-time, cross-sample, and clinical dataset evaluations.

CONCLUSIONS

scMCGraph represents a significant advance in cell type annotation. The experiments have demonstrated that introducing pathway information significantly improves the learning of cell-cell graphs, with their resulting consensus graph enhancing the predictive performance of cell type prediction. Different pathway databases provide complementary data, and an increase in the number of pathways can also boost model performance. Extensive testing shows that in various cross-dataset application scenarios, scMCGraph consistently exhibits both accuracy and robustness.

FOXR2 Targets LHX6+/DLX+ Neural Lineages to Drive Central Nervous System Neuroblastoma

Central nervous system neuroblastoma with forkhead box R2 (FOXR2) activation (NB-FOXR2) is a high-grade tumor of the brain hemispheres and a newly identified molecular entity. Tumors express dual neuronal and glial markers, leading to frequent misdiagnoses, and limited information exists on the role of FOXR2 in their genesis. To identify their cellular origins, we profiled the transcriptomes of NB-FOXR2 tumors at the bulk and single-cell levels and integrated these profiles with large single-cell references of the normal brain. NB-FOXR2 tumors mapped to LHX6+/DLX+ lineages derived from the medial ganglionic eminence, a progenitor domain in the ventral telencephalon. In vivo prenatal Foxr2 targeting to the ganglionic eminences in mice induced postnatal cortical tumors recapitulating human NB-FOXR2-specific molecular signatures. Profiling of FOXR2 binding on chromatin in murine models revealed an association with ETS transcriptional networks, as well as direct binding of FOXR2 at key transcription factors that coordinate initiation of gliogenesis. These data indicate that NB-FOXR2 tumors originate from LHX6+/DLX+ interneuron lineages, a lineage of origin distinct from that of other FOXR2-driven brain tumors, highlight the susceptibility of ventral telencephalon-derived interneurons to FOXR2-driven oncogenesis, and suggest that FOXR2-induced activation of glial programs may explain the mixed neuronal and oligodendroglial features in these tumors. More broadly, this work underscores systematic profiling of brain development as an efficient approach to orient oncogenic targeting for in vivo modeling, critical for the study of rare tumors and development of therapeutics. Significance: Profiling the developing brain enabled rationally guided modeling of FOXR2-activated CNS neuroblastoma, providing a strategy to overcome the heterogeneous origins of pediatric brain tumors that hamper tumor modeling and therapy development. See related commentary by Orr, p. 195.

Single-cell omics: experimental workflow, data analyses and applications

Cells are the fundamental units of biological systems and exhibit unique development trajectories and molecular features. Our exploration of how the genomes orchestrate the formation and maintenance of each cell, and control the cellular phenotypes of various organismsis, is both captivating and intricate. Since the inception of the first single-cell RNA technology, technologies related to single-cell sequencing have experienced rapid advancements in recent years. These technologies have expanded horizontally to include single-cell genome, epigenome, proteome, and metabolome, while vertically, they have progressed to integrate multiple omics data and incorporate additional information such as spatial scRNA-seq and CRISPR screening. Single-cell omics represent a groundbreaking advancement in the biomedical field, offering profound insights into the understanding of complex diseases, including cancers. Here, we comprehensively summarize recent advances in single-cell omics technologies, with a specific focus on the methodology section. This overview aims to guide researchers in selecting appropriate methods for single-cell sequencing and related data analysis.

Profiling cell identity and tissue architecture with single-cell and spatial transcriptomics

Single-cell transcriptomics has broadened our understanding of cellular diversity and gene expression dynamics in healthy and diseased tissues. Recently, spatial transcriptomics has emerged as a tool to contextualize single cells in multicellular neighbourhoods and to identify spatially recurrent phenotypes, or ecotypes. These technologies have generated vast datasets with targeted-transcriptome and whole-transcriptome profiles of hundreds to millions of cells. Such data have provided new insights into developmental hierarchies, cellular plasticity and diverse tissue microenvironments, and spurred a burst of innovation in computational methods for single-cell analysis. In this Review, we discuss recent advancements, ongoing challenges and prospects in identifying and characterizing cell states and multicellular neighbourhoods. We discuss recent progress in sample processing, data integration, identification of subtle cell states, trajectory modelling, deconvolution and spatial analysis. Furthermore, we discuss the increasing application of deep learning, including foundation models, in analysing single-cell and spatial transcriptomics data. Finally, we discuss recent applications of these tools in the fields of stem cell biology, immunology, and tumour biology, and the future of single-cell and spatial transcriptomics in biological research and its translation to the clinic.

scCTS: identifying the cell type-specific marker genes from population-level single-cell RNA-seq

Single-cell RNA-sequencing (scRNA-seq) provides gene expression profiles of individual cells from complex samples, facilitating the detection of cell type-specific marker genes. In scRNA-seq experiments with multiple donors, the population level variation brings an extra layer of complexity in cell type-specific gene detection, for example, they may not appear in all donors. Motivated by this observation, we develop a statistical model named scCTS to identify cell type-specific genes from population-level scRNA-seq data. Extensive data analyses demonstrate that the proposed method identifies more biologically meaningful cell type-specific genes compared to traditional methods.

Graph domain adaptation-based framework for gene expression enhancement and cell type identification in large-scale spatially resolved transcriptomics

Spatially resolved transcriptomics (SRT) technologies facilitate gene expression profiling with spatial resolution in a naïve state. Nevertheless, current SRT technologies exhibit limitations, manifesting as either low transcript detection sensitivity or restricted gene throughput. These constraints result in diminished precision and coverage in gene measurement. In response, we introduce SpaGDA, a sophisticated deep learning-based graph domain adaptation framework for both scenarios of gene expression imputation and cell type identification in spatially resolved transcriptomics data by impartially transferring knowledge from reference scRNA-seq data. Systematic benchmarking analyses across several SRT datasets generated from different technologies have demonstrated SpaGDA's superior effectiveness compared to state-of-the-art methods in both scenarios. Further applied to three SRT datasets of different biological contexts, SpaGDA not only better recovers the well-established knowledge sourced from public atlases and existing scientific literature but also yields a more informative spatial expression pattern of genes. Together, these results demonstrate that SpaGDA can be used to overcome the challenges of current SRT data and provide more accurate insights into biological processes or disease development. The SpaGDA is available in https://github.com/shenrb/SpaGDA.

Progress and opportunities of foundation models in bioinformatics

Bioinformatics has undergone a paradigm shift in artificial intelligence (AI), particularly through foundation models (FMs), which address longstanding challenges in bioinformatics such as limited annotated data and data noise. These AI techniques have demonstrated remarkable efficacy across various downstream validation tasks, effectively representing diverse biological entities and heralding a new era in computational biology. The primary goal of this survey is to conduct a general investigation and summary of FMs in bioinformatics, tracing their evolutionary trajectory, current research landscape, and methodological frameworks. Our primary focus is on elucidating the application of FMs to specific biological problems, offering insights to guide the research community in choosing appropriate FMs for tasks like sequence analysis, structure prediction, and function annotation. Each section delves into the intricacies of the targeted challenges, contrasting the architectures and advancements of FMs with conventional methods and showcasing their utility across different biological domains. Further, this review scrutinizes the hurdles and constraints encountered by FMs in biology, including issues of data noise, model interpretability, and potential biases. This analysis provides a theoretical groundwork for understanding the circumstances under which certain FMs may exhibit suboptimal performance. Lastly, we outline prospective pathways and methodologies for the future development of FMs in biological research, facilitating ongoing innovation in the field. This comprehensive examination not only serves as an academic reference but also as a roadmap for forthcoming explorations and applications of FMs in biology.

Immune landscape of oncohistone-mutant gliomas reveals diverse myeloid populations and tumor-promoting function

Histone H3-mutant gliomas are deadly brain tumors characterized by a dysregulated epigenome and stalled differentiation. In contrast to the extensive datasets available on tumor cells, limited information exists on their tumor microenvironment (TME), particularly the immune infiltrate. Here, we characterize the immune TME of H3.3K27M and G34R/V-mutant gliomas, and multiple H3.3K27M mouse models, using transcriptomic, proteomic and spatial single-cell approaches. Resolution of immune lineages indicates high infiltration of H3-mutant gliomas with diverse myeloid populations, high-level expression of immune checkpoint markers, and scarce lymphoid cells, findings uniformly reproduced in all H3.3K27M mouse models tested. We show these myeloid populations communicate with H3-mutant cells, mediating immunosuppression and sustaining tumor formation and maintenance. Dual inhibition of myeloid cells and immune checkpoint pathways show significant therapeutic benefits in pre-clinical syngeneic mouse models. Our findings provide a valuable characterization of the TME of oncohistone-mutant gliomas, and insight into the means for modulating the myeloid infiltrate for the benefit of patients.

Comparative dissection of transcriptional landscapes of human iPSC-NK differentiation and NK cell development

Clinical and preclinical research has demonstrated that iPSC-derived NK (iNK) cells have a high therapeutic potential, yet poor understanding of the detailed process of their differentiation in vitro and their counterpart cell development in vivo has hindered therapeutic iNK cell production and engineering. Here we dissect the crucial differentiation of both fetal liver NK cells and iNK cells to enable the rational design of advanced iNK production protocols. We use a comparative analysis of single-cell RNA-seq (scRNA-seq) to pinpoint key factors lacking in the induced setting which we hypothesized would hinder iNK differentiation and/ or functionality. By analyzing key transcription factor regulatory networks, we discovered the importance of TBX21, EOMES, and STAT5A in the differentiation timeline. This analysis provides a blueprint for further engineering new iPSC lines to obtain iNK cells with enhanced functions. We validated this approach by creating a new line of STAT5A-iPSCs which can be differentiated to STAT5A-expressing macrophages with both NK cell and macrophage features such as perforin production, phagocytosis, and anti-tumor functions.

Timing and location dictate monocyte fate and their transition to tumor-associated macrophages

Tumor-associated macrophages (TAMs) are a heterogeneous population of cells whose phenotypes and functions are shaped by factors that are incompletely understood. Herein, we asked when and where TAMs arise from blood monocytes and how they evolve during tumor development. We initiated pancreatic ductal adenocarcinoma (PDAC) in inducible monocyte fate-mapping mice and combined single-cell transcriptomics and high-dimensional flow cytometry to profile the monocyte-to-TAM transition. We revealed that monocytes differentiate first into a transient intermediate population of TAMs that generates two longer-lived lineages of terminally differentiated TAMs with distinct gene expression profiles, phenotypes, and intratumoral localization. Transcriptome datasets and tumor samples from patients with PDAC evidenced parallel TAM populations in humans and their prognostic associations. These insights will support the design of new therapeutic strategies targeting TAMs in PDAC.

Detecting novel cell type in single-cell chromatin accessibility data via open-set domain adaptation

Recent advances in single-cell technologies enable the rapid growth of multi-omics data. Cell type annotation is one common task in analyzing single-cell data. It is a challenge that some cell types in the testing set are not present in the training set (i.e. unknown cell types). Most scATAC-seq cell type annotation methods generally assign each cell in the testing set to one known type in the training set but neglect unknown cell types. Here, we present OVAAnno, an automatic cell types annotation method which utilizes open-set domain adaptation to detect unknown cell types in scATAC-seq data. Comprehensive experiments show that OVAAnno successfully identifies known and unknown cell types. Further experiments demonstrate that OVAAnno also performs well on scRNA-seq data. Our codes are available online at https://github.com/lisaber/OVAAnno/tree/master.

ReCIDE: robust estimation of cell type proportions by integrating single-reference-based deconvolutions

In this study, we introduce Robust estimation of Cell type proportions by Integrating single-reference-based DEconvolutions (ReCIDE), an innovative framework for robust estimation of cell type proportions by integrating single-reference-based deconvolutions. ReCIDE outperforms existing approaches in benchmark and real datasets, particularly excelling in estimating rare cell type proportions. Through exploratory analysis on public bulk data of triple-negative breast cancer (TNBC) patients using ReCIDE, we demonstrate a significant correlation between the prognosis of TNBC patients and the proportions of both T cell and perivascular-like cell subtypes. Built upon this discovery, we develop a prognostic assessment model for TNBC patients. Our contribution presents a novel framework for enhancing deconvolution accuracy, showcasing its effectiveness in medical research.

Single-cell discovery of m6A RNA modifications in the hippocampus

N 6-Methyladenosine (m6A) is a prevalent and highly regulated RNA modification essential for RNA metabolism and normal brain function. It is particularly important in the hippocampus, where m6A is implicated in neurogenesis and learning. Although extensively studied, its presence in specific cell types remains poorly understood. We investigated m6A in the hippocampus at a single-cell resolution, revealing a comprehensive landscape of m6A modifications within individual cells. Through our analysis, we uncovered transcripts exhibiting a dense m6A profile, notably linked to neurological disorders such as Alzheimer's disease. Our findings suggest a pivotal role of m6A-containing transcripts, particularly in the context of CAMK2A neurons. Overall, this work provides new insights into the molecular mechanisms underlying hippocampal physiology and lays the foundation for future studies investigating the dynamic nature of m6A RNA methylation in the healthy and diseased brain.

Neuromedin U receptor 1 deletion leads to impaired immunotherapy response and high malignancy in colorectal cancer

Colorectal cancer (CRC) exhibits significant heterogeneity, impacting immunotherapy efficacy, particularly in immune desert subtypes. Neuromedin U receptor 1 (NMUR1) has been reported to perform a vital function in immunity and inflammation. Through comprehensive multi-omics analyses, we have systematically characterized NMUR1 across various tumors, assessing expression patterns, genetic alterations, prognostic significance, immune infiltration, and pathway associations at both the bulk sequencing and single-cell scales. Our findings demonstrate a positive correlation between NMUR1 and CD8+ T cell infiltration, with elevated NMUR1 levels in CD8+ T cells linked to improved immunotherapy outcomes in patients with CRC. Further, we have validated the NMUR1 expression signature in CRC cell lines and patient-derived tissues, revealing its interaction with key immune checkpoints, including lymphocyte activation gene 3 and cytotoxic T-lymphocyte-associated protein 4. Additionally, NMUR1 suppression enhances CRC cell proliferation and invasiveness. Our integrated analyses and experiments open new avenues for personalized immunotherapy strategies in CRC treatment.

Single-cell transcriptomic-informed deconvolution of bulk data identifies immune checkpoint blockade resistance in urothelial cancer

Interactions within the tumor microenvironment (TME) significantly influence tumor progression and treatment responses. While single-cell RNA sequencing (scRNA-seq) and spatial genomics facilitate TME exploration, many clinical cohorts are assessed at the bulk tissue level. Integrating scRNA-seq and bulk tissue RNA-seq data through computational deconvolution is essential for obtaining clinically relevant insights. Our method, ProM, enables the examination of major and minor cell types. Through evaluation against existing methods using paired single-cell and bulk RNA sequencing of human urothelial cancer (UC) samples, ProM demonstrates superiority. Application to UC cohorts treated with immune checkpoint inhibitors reveals pre-treatment cellular features associated with poor outcomes, such as elevated SPP1 expression in macrophage/monocytes (MM). Our deconvolution method and paired single-cell and bulk tissue RNA-seq dataset contribute novel insights into TME heterogeneity and resistance to immune checkpoint blockade.

SIMS: A deep-learning label transfer tool for single-cell RNA sequencing analysis

Cell atlases serve as vital references for automating cell labeling in new samples, yet existing classification algorithms struggle with accuracy. Here we introduce SIMS (scalable, interpretable machine learning for single cell), a low-code data-efficient pipeline for single-cell RNA classification. We benchmark SIMS against datasets from different tissues and species. We demonstrate SIMS's efficacy in classifying cells in the brain, achieving high accuracy even with small training sets (<3,500 cells) and across different samples. SIMS accurately predicts neuronal subtypes in the developing brain, shedding light on genetic changes during neuronal differentiation and postmitotic fate refinement. Finally, we apply SIMS to single-cell RNA datasets of cortical organoids to predict cell identities and uncover genetic variations between cell lines. SIMS identifies cell-line differences and misannotated cell lineages in human cortical organoids derived from different pluripotent stem cell lines. Altogether, we show that SIMS is a versatile and robust tool for cell-type classification from single-cell datasets.

Patient-derived tumor-like cell clusters for personalized chemo- and immunotherapies in non-small cell lung cancer

Many patient-derived tumor models have emerged recently. However, their potential to guide personalized drug selection remains unclear. Here, we report patient-derived tumor-like cell clusters (PTCs) for non-small cell lung cancer (NSCLC), capable of conducting 100-5,000 drug tests within 10 days. We have established 283 PTC models with an 81% success rate. PTCs contain primary tumor epithelium self-assembled with endogenous stromal and immune cells and show a high degree of similarity to the original tumors in phenotypic and genotypic features. Utilizing standardized culture and drug-response assessment protocols, PTC drug-testing assays reveal 89% overall consistency in prospectively predicting clinical outcomes, with 98.1% accuracy distinguishing complete/partial response from progressive disease. Notably, PTCs enable accurate prediction of clinical outcomes for patients undergoing anti-PD1 therapy by combining cell viability and IFN-γ value assessments. These findings suggest that PTCs could serve as a valuable preclinical model for personalized medicine and basic research in NSCLC.

Single-cell RNA sequencing to explore cancer-associated fibroblasts heterogeneity: "Single" vision for "heterogeneous" environment

Cancer-associated fibroblasts (CAFs), a phenotypically and functionally heterogeneous stromal cell, are one of the most important components of the tumour microenvironment. Previous studies have consolidated it as a promising target against cancer. However, variable therapeutic efficacy-both protumor and antitumor effects have been observed not least owing to the strong heterogeneity of CAFs. Over the past 10 years, advances in single-cell RNA sequencing (scRNA-seq) technologies had a dramatic effect on biomedical research, enabling the analysis of single cell transcriptomes with unprecedented resolution and throughput. Specifically, scRNA-seq facilitates our understanding of the complexity and heterogeneity of diverse CAF subtypes. In this review, we discuss the up-to-date knowledge about CAF heterogeneity with a focus on scRNA-seq perspective to investigate the emerging strategies for integrating multimodal single-cell platforms. Furthermore, we summarized the clinical application of scRNA-seq on CAF research. We believe that the comprehensive understanding of the heterogeneity of CAFs form different visions will generate innovative solutions to cancer therapy and achieve clinical applications.

Single-cell and spatial RNA sequencing reveal the spatiotemporal trajectories of fruit senescence

The senescence of fruit is a complex physiological process, with various cell types within the pericarp, making it highly challenging to elucidate their individual roles in fruit senescence. In this study, a single-cell expression atlas of the pericarp of pitaya (Hylocereus undatus) is constructed, revealing exocarp and mesocarp cells undergoing the most significant changes during the fruit senescence process. Pseudotime analysis establishes cellular differentiation and gene expression trajectories during senescence. Early-stage oxidative stress imbalance is followed by the activation of resistance in exocarp cells, subsequently senescence-associated proteins accumulate in the mesocarp cells at late-stage senescence. The central role of the early response factor HuCMB1 is unveiled in the senescence regulatory network. This study provides a spatiotemporal perspective for a deeper understanding of the dynamic senescence process in plants.

Lipidomics reveals new lipid-based lung adenocarcinoma early diagnosis model

Lung adenocarcinoma (LUAD) continues to pose a significant mortality risk with a lack of dependable biomarkers for early noninvasive cancer detection. Here, we find that aberrant lipid metabolism is significantly enriched in lung cancer cells. Further, we identified four signature lipids highly associated with LUAD and developed a lipid signature-based scoring model (LSRscore). Evaluation of LSRscore in a discovery cohort reveals a robust predictive capability for LUAD (AUC: 0.972), a result further validated in an independent cohort (AUC: 0.92). We highlight one lipid signature biomarker, PE(18:0/18:1), consistently exhibiting altered levels both in cancer tissue and in plasma of LUAD patients, demonstrating significant predictive power for early-stage LUAD. Transcriptome analysis reveals an association between increased PE(18:0/18:1) levels and dysregulated glycerophospholipid metabolism, which consistently displays strong prognostic value across two LUAD cohorts. The combined utility of LSRscore and PE(18:0/18:1) holds promise for early-stage diagnosis and prognosis of LUAD.

scAnnoX: an R package integrating multiple public tools for single-cell annotation

Background

Single-cell annotation plays a crucial role in the analysis of single-cell genomics data. Despite the existence of numerous single-cell annotation algorithms, a comprehensive tool for integrating and comparing these algorithms is also lacking.

Methods

This study meticulously investigated a plethora of widely adopted single-cell annotation algorithms. Ten single-cell annotation algorithms were selected based on the classification of either reference dataset-dependent or marker gene-dependent approaches. These algorithms included SingleR, Seurat, sciBet, scmap, CHETAH, scSorter, sc.type, cellID, scCATCH, and SCINA. Building upon these algorithms, we developed an R package named scAnnoX for the integration and comparative analysis of single-cell annotation algorithms.

Results

The development of the scAnnoX software package provides a cohesive framework for annotating cells in scRNA-seq data, enabling researchers to more efficiently perform comparative analyses among the cell type annotations contained in scRNA-seq datasets. The integrated environment of scAnnoX streamlines the testing, evaluation, and comparison processes among various algorithms. Among the ten annotation tools evaluated, SingleR, Seurat, sciBet, and scSorter emerged as top-performing algorithms in terms of prediction accuracy, with SingleR and sciBet demonstrating particularly superior performance, offering guidance for users. Interested parties can access the scAnnoX package at https://github.com/XQ-hub/scAnnoX.

Childhood cancer mutagenesis caused by transposase-derived PGBD5

Genomic rearrangements are a hallmark of most childhood tumors, including medulloblastoma, one of the most common brain tumors in children, but their causes remain largely unknown. Here, we show that PiggyBac transposable element derived 5 (Pgbd5) promotes tumor development in multiple developmentally accurate mouse models of Sonic Hedgehog (SHH) medulloblastoma. Most Pgbd5-deficient mice do not develop tumors, while maintaining normal cerebellar development. Ectopic activation of SHH signaling is sufficient to enforce cerebellar granule cell progenitor-like cell states, which exhibit Pgbd5-dependent expression of distinct DNA repair and neurodevelopmental factors. Mouse medulloblastomas expressing Pgbd5 have increased numbers of somatic structural DNA rearrangements, some of which carry PGBD5-specific sequences at their breakpoints. Similar sequence breakpoints recurrently affect somatic DNA rearrangements of known tumor suppressors and oncogenes in medulloblastomas in 329 children. This identifies PGBD5 as a medulloblastoma mutator and provides a genetic mechanism for the generation of oncogenic DNA rearrangements in childhood cancer.

scCancer2: data-driven in-depth annotations of the tumor microenvironment at single-level resolution

SUMMARY

Single-cell RNA-seq (scRNA-seq) is a powerful technique for decoding the complex cellular compositions in the tumor microenvironment (TME). As previous studies have defined many meaningful cell subtypes in several tumor types, there is a great need to computationally transfer these labels to new datasets. Also, different studies used different approaches or criteria to define the cell subtypes for the same major cell lineages. The relationships between the cell subtypes defined in different studies should be carefully evaluated. In this updated package scCancer2, designed for integrative tumor scRNA-seq data analysis, we developed a supervised machine learning framework to annotate TME cells with annotated cell subtypes from 15 scRNA-seq datasets with 594 samples in total. Based on the trained classifiers, we quantitatively constructed the similarity maps between the cell subtypes defined in different references by testing on all the 15 datasets. Secondly, to improve the identification of malignant cells, we designed a classifier by integrating large-scale pan-cancer TCGA bulk gene expression datasets and scRNA-seq datasets (10 cancer types, 175 samples, 663 857 cells). This classifier shows robust performances when no internal confidential reference cells are available. Thirdly, scCancer2 integrated a module to process the spatial transcriptomic data and analyze the spatial features of TME.

AVAILABILITY AND IMPLEMENTATION

The package and user documentation are available at http://lifeome.net/software/sccancer2/ and https://doi.org/10.5281/zenodo.10477296.

scMMT: a multi-use deep learning approach for cell annotation, protein prediction and embedding in single-cell RNA-seq data

Accurate cell type annotation in single-cell RNA-sequencing data is essential for advancing biological and medical research, particularly in understanding disease progression and tumor microenvironments. However, existing methods are constrained by single feature extraction approaches, lack of adaptability to immune cell types with similar molecular profiles but distinct functions and a failure to account for the impact of cell label noise on model accuracy, all of which compromise the precision of annotation. To address these challenges, we developed a supervised approach called scMMT. We proposed a novel feature extraction technique to uncover more valuable information. Additionally, we constructed a multi-task learning framework based on the GradNorm method to enhance the recognition of challenging immune cells and reduce the impact of label noise by facilitating mutual reinforcement between cell type annotation and protein prediction tasks. Furthermore, we introduced logarithmic weighting and label smoothing mechanisms to enhance the recognition ability of rare cell types and prevent model overconfidence. Through comprehensive evaluations on multiple public datasets, scMMT has demonstrated state-of-the-art performance in various aspects including cell type annotation, rare cell identification, dropout and label noise resistance, protein expression prediction and low-dimensional embedding representation.

Presence of onco-fetal neighborhoods in hepatocellular carcinoma is associated with relapse and response to immunotherapy

Onco-fetal reprogramming of the tumor ecosystem induces fetal developmental signatures in the tumor microenvironment, leading to immunosuppressive features. Here, we employed single-cell RNA sequencing, spatial transcriptomics and bulk RNA sequencing to delineate specific cell subsets involved in hepatocellular carcinoma (HCC) relapse and response to immunotherapy. We identified POSTN+ extracellular matrix cancer-associated fibroblasts (EM CAFs) as a prominent onco-fetal interacting hub, promoting tumor progression. Cell-cell communication and spatial transcriptomics analysis revealed crosstalk and co-localization of onco-fetal cells, including POSTN+ CAFs, FOLR2+ macrophages and PLVAP+ endothelial cells. Further analyses suggest an association between onco-fetal reprogramming and epithelial-mesenchymal transition (EMT), tumor cell proliferation and recruitment of Treg cells, ultimately influencing early relapse and response to immunotherapy. In summary, our study identifies POSTN+ CAFs as part of the HCC onco-fetal niche and highlights its potential influence in EMT, relapse and immunotherapy response, paving the way for the use of onco-fetal signatures for therapeutic stratification.

scPML: pathway-based multi-view learning for cell type annotation from single-cell RNA-seq data

Recent developments in single-cell technology have enabled the exploration of cellular heterogeneity at an unprecedented level, providing invaluable insights into various fields, including medicine and disease research. Cell type annotation is an essential step in its omics research. The mainstream approach is to utilize well-annotated single-cell data to supervised learning for cell type annotation of new singlecell data. However, existing methods lack good generalization and robustness in cell annotation tasks, partially due to difficulties in dealing with technical differences between datasets, as well as not considering the heterogeneous associations of genes in regulatory mechanism levels. Here, we propose the scPML model, which utilizes various gene signaling pathway data to partition the genetic features of cells, thus characterizing different interaction maps between cells. Extensive experiments demonstrate that scPML performs better in cell type annotation and detection of unknown cell types from different species, platforms, and tissues.

Hierarchical and automated cell-type annotation and inference of cancer cell of origin with Census

MOTIVATION

Cell-type annotation is a time-consuming yet critical first step in the analysis of single-cell RNA-seq data, especially when multiple similar cell subtypes with overlapping marker genes are present. Existing automated annotation methods have a number of limitations, including requiring large reference datasets, high computation time, shallow annotation resolution, and difficulty in identifying cancer cells or their most likely cell of origin.

RESULTS

We developed Census, a biologically intuitive and fully automated cell-type identification method for single-cell RNA-seq data that can deeply annotate normal cells in mammalian tissues and identify malignant cells and their likely cell of origin. Motivated by the inherently stratified developmental programs of cellular differentiation, Census infers hierarchical cell-type relationships and uses gradient-boosted \decision trees that capitalize on nodal cell-type relationships to achieve high prediction speed and accuracy. When benchmarked on 44 atlas-scale normal and cancer, human and mouse tissues, Census significantly outperforms state-of-the-art methods across multiple metrics and naturally predicts the cell-of-origin of different cancers. Census is pretrained on the Tabula Sapiens to classify 175 cell-types from 24 organs; however, users can seamlessly train their own models for customized applications.

AVAILABILITY AND IMPLEMENTATION

Census is available at Zenodo https://zenodo.org/records/7017103 and on our Github https://github.com/sjdlabgroup/Census.

Single-Cell RNA-Sequencing Provides Insight into Skeletal Muscle Evolution during the Selection of Muscle Characteristics

Skeletal muscle comprises a large, heterogeneous assortment of cell populations that interact to maintain muscle homeostasis, but little is known about the mechanism that controls myogenic development in response to artificial selection. Different pig (Sus scrofa) breeds exhibit distinct muscle phenotypes resulting from domestication and selective breeding. Using unbiased single-cell transcriptomic sequencing analysis (scRNA-seq), the impact of artificial selection on cell profiles is investigated in neonatal skeletal muscle of pigs. This work provides panoramic muscle-resident cell profiles and identifies novel and breed-specific cells, mapping them on pseudotime trajectories. Artificial selection has elicited significant changes in muscle-resident cell profiles, while conserving signs of generational environmental challenges. These results suggest that fibro-adipogenic progenitors serve as a cellular interaction hub and that specific transcription factors identified here may serve as candidate target regulons for the pursuit of a specific muscle phenotype. Furthermore, a cross-species comparison of humans, mice, and pigs illustrates the conservation and divergence of mammalian muscle ontology. The findings of this study reveal shifts in cellular heterogeneity, novel cell subpopulations, and their interactions that may greatly facilitate the understanding of the mechanism underlying divergent muscle phenotypes arising from artificial selection.

Unraveling Neuronal Identities Using SIMS: A Deep Learning Label Transfer Tool for Single-Cell RNA Sequencing Analysis

Large single-cell RNA datasets have contributed to unprecedented biological insight. Often, these take the form of cell atlases and serve as a reference for automating cell labeling of newly sequenced samples. Yet, classification algorithms have lacked the capacity to accurately annotate cells, particularly in complex datasets. Here we present SIMS (Scalable, Interpretable Machine Learning for Single-Cell), an end-to-end data-efficient machine learning pipeline for discrete classification of single-cell data that can be applied to new datasets with minimal coding. We benchmarked SIMS against common single-cell label transfer tools and demonstrated that it performs as well or better than state of the art algorithms. We then use SIMS to classify cells in one of the most complex tissues: the brain. We show that SIMS classifies cells of the adult cerebral cortex and hippocampus at a remarkably high accuracy. This accuracy is maintained in trans-sample label transfers of the adult human cerebral cortex. We then apply SIMS to classify cells in the developing brain and demonstrate a high level of accuracy at predicting neuronal subtypes, even in periods of fate refinement, shedding light on genetic changes affecting specific cell types across development. Finally, we apply SIMS to single cell datasets of cortical organoids to predict cell identities and unveil genetic variations between cell lines. SIMS identifies cell-line differences and misannotated cell lineages in human cortical organoids derived from different pluripotent stem cell lines. When cell types are obscured by stress signals, label transfer from primary tissue improves the accuracy of cortical organoid annotations, serving as a reliable ground truth. Altogether, we show that SIMS is a versatile and robust tool for cell-type classification from single-cell datasets.

Decoding Human Biology and Disease Using Single-cell Omics Technologies

Over the past decade, advances in single-cell omics (SCO) technologies have enabled the investigation of cellular heterogeneity at an unprecedented resolution and scale, opening a new avenue for understanding human biology and disease. In this review, we summarize the developments of sequencing-based SCO technologies and computational methods, and focus on considerable insights acquired from SCO sequencing studies to understand normal and diseased properties, with a particular emphasis on cancer research. We also discuss the technological improvements of SCO and its possible contribution to fundamental research of the human, as well as its great potential in clinical diagnoses and personalized therapies of human disease.

Single-cell assignment using multiple-adversarial domain adaptation network with large-scale references

The rapid accumulation of single-cell RNA-seq data has provided rich resources to characterize various human cell populations. However, achieving accurate cell-type annotation using public references presents challenges due to inconsistent annotations, batch effects, and rare cell types. Here, we introduce SELINA (single-cell identity navigator), an integrative and automatic cell-type annotation framework based on a pre-curated reference atlas spanning various tissues. SELINA employs a multiple-adversarial domain adaptation network to remove batch effects within the reference dataset. Additionally, it enhances the annotation of less frequent cell types by synthetic minority oversampling and fits query data with the reference data using an autoencoder. SELINA culminates in the creation of a comprehensive and uniform reference atlas, encompassing 1.7 million cells covering 230 distinct human cell types. We substantiate its robustness and superiority across a multitude of human tissues. Notably, SELINA could accurately annotate cells within diverse disease contexts. SELINA provides a complete solution for human single-cell RNA-seq data annotation with both python and R packages.

scMHNN: a novel hypergraph neural network for integrative analysis of single-cell epigenomic, transcriptomic and proteomic data

Technological advances have now made it possible to simultaneously profile the changes of epigenomic, transcriptomic and proteomic at the single cell level, allowing a more unified view of cellular phenotypes and heterogeneities. However, current computational tools for single-cell multi-omics data integration are mainly tailored for bi-modality data, so new tools are urgently needed to integrate tri-modality data with complex associations. To this end, we develop scMHNN to integrate single-cell multi-omics data based on hypergraph neural network. After modeling the complex data associations among various modalities, scMHNN performs message passing process on the multi-omics hypergraph, which can capture the high-order data relationships and integrate the multiple heterogeneous features. Followingly, scMHNN learns discriminative cell representation via a dual-contrastive loss in self-supervised manner. Based on the pretrained hypergraph encoder, we further introduce the pre-training and fine-tuning paradigm, which allows more accurate cell-type annotation with only a small number of labeled cells as reference. Benchmarking results on real and simulated single-cell tri-modality datasets indicate that scMHNN outperforms other competing methods on both cell clustering and cell-type annotation tasks. In addition, we also demonstrate scMHNN facilitates various downstream tasks, such as cell marker detection and enrichment analysis.

GOWDL: gene ontology-driven wide and deep learning model for cell typing of scRNA-seq data

Single-cell RNA-sequencing (scRNA-seq) allows for obtaining genomic and transcriptomic profiles of individual cells. That data make it possible to characterize tissues at the cell level. In this context, one of the main analyses exploiting scRNA-seq data is identifying the cell types within tissue to estimate the quantitative composition of cell populations. Due to the massive amount of available scRNA-seq data, automatic classification approaches for cell typing, based on the most recent deep learning technology, are needed. Here, we present the gene ontology-driven wide and deep learning (GOWDL) model for classifying cell types in several tissues. GOWDL implements a hybrid architecture that considers the functional annotations found in Gene Ontology and the marker genes typical of specific cell types. We performed cross-validation and independent external testing, comparing our algorithm with 12 other state-of-the-art predictors. Classification scores demonstrated that GOWDL reached the best results over five different tissues, except for recall, where we got about 92% versus 97% of the best tool. Finally, we presented a case study on classifying immune cell populations in breast cancer using a hierarchical approach based on GOWDL.

CellAnn: a comprehensive, super-fast, and user-friendly single-cell annotation web server

MOTIVATION

Single-cell sequencing technology has become a routine in studying many biological problems. A core step of analyzing single-cell data is the assignment of cell clusters to specific cell types. Reference-based methods are proposed for predicting cell types for single-cell clusters. However, the scalability and lack of preprocessed reference datasets prevent them from being practical and easy to use.

RESULTS

Here, we introduce a reference-based cell annotation web server, CellAnn, which is super-fast and easy to use. CellAnn contains a comprehensive reference database with 204 human and 191 mouse single-cell datasets. These reference datasets cover 32 organs. Furthermore, we developed a cluster-to-cluster alignment method to transfer cell labels from the reference to the query datasets, which is superior to the existing methods with higher accuracy and higher scalability. Finally, CellAnn is an online tool that integrates all the procedures in cell annotation, including reference searching, transferring cell labels, visualizing results, and harmonizing cell annotation labels. Through the user-friendly interface, users can identify the best annotation by cross-validating with multiple reference datasets. We believe that CellAnn can greatly facilitate single-cell sequencing data analysis.

AVAILABILITY AND IMPLEMENTATION

The web server is available at www.cellann.io, and the source code is available at https://github.com/Pinlyu3/CellAnn_shinyapp.

Uncovering a mammalian neural-specific poly(A) binding protein with unique properties

The mRNA 3' poly(A) tail plays a critical role in regulating both mRNA translation and turnover. It is bound by the cytoplasmic poly(A) binding protein (PABPC), an evolutionarily conserved protein that can interact with translation factors and mRNA decay machineries to regulate gene expression. Mammalian PABPC1, the prototypical PABPC, is expressed in most tissues and interacts with eukaryotic translation initiation factor 4G (eIF4G) to stimulate translation in specific contexts. In this study, we uncovered a new mammalian PABPC, which we named neural PABP (neuPABP), as it is predominantly expressed in the brain. neuPABP maintains a unique architecture as compared with other PABPCs, containing only two RNA recognition motifs (RRMs) and maintaining a unique N-terminal domain of unknown function. neuPABP expression is activated in neurons as they mature during synaptogenesis, where neuPABP localizes to the soma and postsynaptic densities. neuPABP interacts with the noncoding RNA BC1, as well as mRNAs coding for ribosomal and mitochondrial proteins. However, in contrast to PABPC1, neuPABP does not associate with actively translating mRNAs in the brain. In keeping with this, we show that neuPABP has evolved such that it does not bind eIF4G and as a result fails to support protein synthesis in vitro. Taken together, these results indicate that mammals have expanded their PABPC repertoire in the brain and propose that neuPABP may support the translational repression of select mRNAs.

A Review of Single-Cell RNA-Seq Annotation, Integration, and Cell-Cell Communication

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular biology at an unprecedented resolution, enabling the characterization of cellular heterogeneity, identification of rare but significant cell types, and exploration of cell-cell communications and interactions. Its broad applications span both basic and clinical research domains. In this comprehensive review, we survey the current landscape of scRNA-seq analysis methods and tools, focusing on count modeling, cell-type annotation, data integration, including spatial transcriptomics, and the inference of cell-cell communication. We review the challenges encountered in scRNA-seq analysis, including issues of sparsity or low expression, reliability of cell annotation, and assumptions in data integration, and discuss the potential impact of suboptimal clustering and differential expression analysis tools on downstream analyses, particularly in identifying cell subpopulations. Finally, we discuss recent advancements and future directions for enhancing scRNA-seq analysis. Specifically, we highlight the development of novel tools for annotating single-cell data, integrating and interpreting multimodal datasets covering transcriptomics, epigenomics, and proteomics, and inferring cellular communication networks. By elucidating the latest progress and innovation, we provide a comprehensive overview of the rapidly advancing field of scRNA-seq analysis.

PreCanCell: An ensemble learning algorithm for predicting cancer and non-cancer cells from single-cell transcriptomes

We propose PreCanCell, a novel algorithm for predicting malignant and non-malignant cells from single-cell transcriptomes. PreCanCell first identifies the differentially expressed genes (DEGs) between malignant and non-malignant cells commonly in five common cancer types-associated single-cell transcriptome datasets. The five common cancer types include renal cell carcinoma (RCC), head and neck squamous cell carcinoma (HNSCC), melanoma, lung adenocarcinoma (LUAD), and breast cancer (BC). With each of the five datasets as the training set and the DEGs as the features, a single cell is classified as malignant or non-malignant by k-NN (k = 5). Finally, the single cell is determined as malignant or non-malignant by the majority vote of the five k-NN classification results. We tested the predictive performance of PreCanCell in 19 single-cell datasets, and reported classification accuracy, sensitivity, specificity, balanced accuracy (the average of sensitivity and specificity) and the area under the receiver operating characteristic curve (AUROC). In all these datasets, PreCanCell achieved above 0.8 accuracy, sensitivity, specificity, balanced accuracy and AUROC. Finally, we compared the predictive performance of PreCanCell with that of seven other algorithms, including CHETAH, SciBet, SCINA, scmap-cell, scmap-cluster, SingleR, and ikarus. Compared to these algorithms, PreCanCell displays the advantages of higher accuracy and simpler implementation. We have developed an R package for the PreCanCell algorithm, which is available at https://github.com/WangX-Lab/PreCanCell.

A Compendium of Syngeneic, Transplantable Pediatric High-Grade Glioma Models Reveals Subtype-Specific Therapeutic Vulnerabilities

Pediatric high-grade gliomas (pHGG) are lethal, incurable brain tumors frequently driven by clonal mutations in histone genes. They often harbor a range of additional genetic alterations that correlate with different ages, anatomic locations, and tumor subtypes. We developed models representing 16 pHGG subtypes driven by different combinations of alterations targeted to specific brain regions. Tumors developed with varying latencies and cell lines derived from these models engrafted in syngeneic, immunocompetent mice with high penetrance. Targeted drug screening revealed unexpected selective vulnerabilities-H3.3G34R/PDGFRAC235Y to FGFR inhibition, H3.3K27M/PDGFRAWT to PDGFRA inhibition, and H3.3K27M/PDGFRAWT and H3.3K27M/PPM1DΔC/PIK3CAE545K to combined inhibition of MEK and PIK3CA. Moreover, H3.3K27M tumors with PIK3CA, NF1, and FGFR1 mutations were more invasive and harbored distinct additional phenotypes, such as exophytic spread, cranial nerve invasion, and spinal dissemination. Collectively, these models reveal that different partner alterations produce distinct effects on pHGG cellular composition, latency, invasiveness, and treatment sensitivity.

SIGNIFICANCE

Histone-mutant pediatric gliomas are a highly heterogeneous tumor entity. Different histone mutations correlate with different ages of onset, survival outcomes, brain regions, and partner alterations. We have developed models of histone-mutant gliomas that reflect this anatomic and genetic heterogeneity and provide evidence of subtype-specific biology and therapeutic targeting. See related commentary by Lubanszky and Hawkins, p. 1516. This article is highlighted in the In This Issue feature, p. 1501.

scAnno: a deconvolution strategy-based automatic cell type annotation tool for single-cell RNA-sequencing data sets

Undoubtedly, single-cell RNA sequencing (scRNA-seq) has changed the research landscape by providing insights into heterogeneous, complex and rare cell populations. Given that more such data sets will become available in the near future, their accurate assessment with compatible and robust models for cell type annotation is a prerequisite. Considering this, herein, we developed scAnno (scRNA-seq data annotation), an automated annotation tool for scRNA-seq data sets primarily based on the single-cell cluster levels, using a joint deconvolution strategy and logistic regression. We explicitly constructed a reference profile for human (30 cell types and 50 human tissues) and a reference profile for mouse (26 cell types and 50 mouse tissues) to support this novel methodology (scAnno). scAnno offers a possibility to obtain genes with high expression and specificity in a given cell type as cell type-specific genes (marker genes) by combining co-expression genes with seed genes as a core. Of importance, scAnno can accurately identify cell type-specific genes based on cell type reference expression profiles without any prior information. Particularly, in the peripheral blood mononuclear cell data set, the marker genes identified by scAnno showed cell type-specific expression, and the majority of marker genes matched exactly with those included in the CellMarker database. Besides validating the flexibility and interpretability of scAnno in identifying marker genes, we also proved its superiority in cell type annotation over other cell type annotation tools (SingleR, scPred, CHETAH and scmap-cluster) through internal validation of data sets (average annotation accuracy: 99.05%) and cross-platform data sets (average annotation accuracy: 95.56%). Taken together, we established the first novel methodology that utilizes a deconvolution strategy for automated cell typing and is capable of being a significant application in broader scRNA-seq analysis. scAnno is available at https://github.com/liuhong-jia/scAnno.

SigPrimedNet: A Signaling-Informed Neural Network for scRNA-seq Annotation of Known and Unknown Cell Types

Single-cell RNA sequencing is increasing our understanding of the behavior of complex tissues or organs, by providing unprecedented details on the complex cell type landscape at the level of individual cells. Cell type definition and functional annotation are key steps to understanding the molecular processes behind the underlying cellular communication machinery. However, the exponential growth of scRNA-seq data has made the task of manually annotating cells unfeasible, due not only to an unparalleled resolution of the technology but to an ever-increasing heterogeneity of the data. Many supervised and unsupervised methods have been proposed to automatically annotate cells. Supervised approaches for cell-type annotation outperform unsupervised methods except when new (unknown) cell types are present. Here, we introduce SigPrimedNet an artificial neural network approach that leverages (i) efficient training by means of a sparsity-inducing signaling circuits-informed layer, (ii) feature representation learning through supervised training, and (iii) unknown cell-type identification by fitting an anomaly detection method on the learned representation. We show that SigPrimedNet can efficiently annotate known cell types while keeping a low false-positive rate for unseen cells across a set of publicly available datasets. In addition, the learned representation acts as a proxy for signaling circuit activity measurements, which provide useful estimations of the cell functionalities.

deCS: A Tool for Systematic Cell Type Annotations of Single-cell RNA Sequencing Data among Human Tissues

Single-cell RNA sequencing (scRNA-seq) is revolutionizing the study of complex and dynamic cellular mechanisms. However, cell type annotation remains a main challenge as it largely relies on a priori knowledge and manual curation, which is cumbersome and subjective. The increasing number of scRNA-seq datasets, as well as numerous published genetic studies, has motivated us to build a comprehensive human cell type reference atlas.Here, we present decoding Cell type Specificity (deCS), an automatic cell type annotation method augmented by a comprehensive collection of human cell type expression profiles and marker genes. We used deCS to annotate scRNA-seq data from various tissue types and systematically evaluated the annotation accuracy under different conditions, including reference panels, sequencing depth, and feature selection strategies. Our results demonstrate that expanding the references is critical for improving annotation accuracy. Compared to many existing state-of-the-art annotation tools, deCS significantly reduced computation time and increased accuracy. deCS can be integrated into the standard scRNA-seq analytical pipeline to enhance cell type annotation. Finally, we demonstrated the broad utility of deCS to identify trait-cell type associations in 51 human complex traits, providing deep insights into the cellular mechanisms underlying disease pathogenesis. All documents for deCS, including source code, user manual, demo data, and tutorials, are freely available at https://github.com/bsml320/deCS.