scSorter: assigning cells to known cell types according to marker genes

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.

A self-training interpretable cell type annotation framework using specific marker gene

MOTIVATION

Recent advances in sequencing technology provide opportunities to study biological processes at a higher resolution. Cell type annotation is an important step in scRNA-seq analysis, which often relies on established marker genes. However, most of the previous methods divide the identification of cell types into two stages, clustering and assignment, whose performances are susceptible to the clustering algorithm, and the marker information cannot effectively guide the clustering process. Furthermore, their linear heuristic-based cell assignment process is often insufficient to capture potential dependencies between cells and types.

RESULTS

Here, we present Interpretable Cell Type Annotation based on self-training (sICTA), a marker-based cell type annotation method that combines the self-training strategy with pseudo-labeling and the nonlinear association capturing capability of Transformer. In addition, we incorporate biological priori knowledge of genes and pathways into the classifier through an attention mechanism to enhance the transparency of the model. A benchmark analysis on 11 publicly available single-cell datasets demonstrates the superiority of sICTA compared to state-of-the-art methods. The robustness of our method is further validated by evaluating the prediction accuracy of the model on different cell types for each single-cell data. Moreover, ablation studies show that self-training and the ability to capture potential dependencies between cells and cell types, both of which are mutually reinforcing, work together to improve model performance. Finally, we apply sICTA to the pancreatic dataset, exemplifying the interpretable attention matrix captured by sICTA.

AVAILABILITY AND IMPLEMENTATION

The source code of sICTA is available in public at https://github.com/nbnbhwyy/sICTA. The processed datasets can be found at https://drive.google.com/drive/folders/1jbqSxacL_IDIZ4uPjq220C9Kv024m9eL. The final version of the model will be permanently available at https://doi.org/10.5281/zenodo.13474010.

Single Cell Expression Analysis of Ductal Carcinoma in Situ Identifies Complex Genotypic-Phenotypic Relationships Altering Epithelial Composition

To identify mechanisms underlying the growth of ductal carcinoma in situ (DCIS) and properties that lead to progression to invasive cancer, we performed single-cell RNA-sequencing (scRNA-seq) on DCIS lesions and matched synchronous normal breast tissue. Using inferred copy number variations (CNV), we identified neoplastic epithelial cells from the clinical specimens which contained a mixture of DCIS and normal ducts. Phylogenetic analysis based on the CNVs demonstrated intratumoral clonal heterogeneity was associated with significant gene expression differences. We also classified epithelial cells into mammary cell states and found that individual genetic clones contained a mixture of cell states suggesting an ongoing pattern of differentiation after neoplastic transformation. Cell state proportions were significantly different based on estrogen receptor (ER) expression with ER-DCIS more closely resembling the distribution in the normal breast, particularly with respect to cells with basal characteristics. Using deconvolution from bulk RNA-seq in archival DCIS specimens, we show that specific alterations in cell state proportions are associated with progression to invasive cancer. Loss of an intact basement membrane (BM) is the functional definition of invasive breast cancer (IBC) and scRNA-seq data demonstrated that ongoing transcription of key BM genes occurs in specific subsets of epithelial cell states. Examining BM in archival microinvasive breast cancers and an in vitro model of invasion, we found that passive loss of BM gene expression due to cell state proportion alterations is associated with loss of the structural integrity of the duct leading to an invasive phenotype. Our analyses provide detailed insight into DCIS biology.

SIGNIFICANCE

Single cell analysis reveals that preinvasive breast cancer is comprised of multiple genetic clones and there is substantial phenotypic diversity both within and between these clones. Ductal carcinoma in situ (DCIS) of the breast is a non-invasive condition commonly identified through mammographic screening. A primary diagnosis of DCIS carries little mortality risk on its own, but its presence is a risk factor for subsequent clonally related invasive breast cancer (IBC) (1-5).

Deconvolving organogenesis in space and time via spatial transcriptomics in thick tissues

Organ development is guided by a space-time landscape that constraints cell behavior. This landscape is challenging to characterize for the hair follicle - the most abundant mini organ - due to its complex microscopic structure and asynchronous development. We developed 3DEEP, a tissue clearing and spatial transcriptomic strategy for characterizing tissue blocks up to 400 µm in thickness. We captured 371 hair follicles at different stages of organogenesis in 1 mm 3 of skin of a 12-hour-old mouse with 6 million transcripts from 81 genes. From this single time point, we deconvoluted follicles by age based on whole-organ molecular pseudotimes to animate a stop-motion 3D atlas of follicle development along its trajectory. We defined molecular stages for hair follicle organogenesis and characterized the order of emergence for its structures, differential signaling dynamics at its top and bottom, morphogen shifts preceding and accompanying structural changes, and series of structural changes leading to the formation of its canal and opening. We further found that hair follicle stem cells and their niche are established and stratified early in organogenesis, before the formation of the hair bulb. Overall, this work demonstrates the power of increased depth of spatial transcriptomics to provide a four-dimensional analysis of organogenesis.

The use of ectopic volar fibroblasts to modify skin identity

Skin identity is controlled by intrinsic features of the epidermis and dermis and their interactions. Modifying skin identity has clinical potential, such as the conversion of residual limb and stump (nonvolar) skin of amputees to pressure-responsive palmoplantar (volar) skin to enhance prosthesis use and minimize skin breakdown. Greater keratin 9 (KRT9) expression, higher epidermal thickness, keratinocyte cytoplasmic size, collagen length, and elastin are markers of volar skin and likely contribute to volar skin resiliency. Given fibroblasts' capacity to modify keratinocyte differentiation, we hypothesized that volar fibroblasts influence these features. Bioprinted skin constructs confirmed the capacity of volar fibroblasts to induce volar keratinocyte features. A clinical trial of healthy volunteers demonstrated that injecting volar fibroblasts into nonvolar skin increased volar features that lasted up to 5 months, highlighting a potential cellular therapy.

Scywalker: scalable end-to-end data analysis workflow for long-read single-cell transcriptome sequencing

MOTIVATION

Existing nanopore single-cell data analysis tools showed severe limitations in handling current data sizes.

RESULTS

We introduce scywalker, an innovative and scalable package developed to comprehensively analyze long-read sequencing data of full-length single-cell or single-nuclei cDNA. We developed novel scalable methods for cell barcode demultiplexing and single-cell isoform calling and quantification and incorporated these in an easily deployable package. Scywalker streamlines the entire analysis process, from sequenced fragments in FASTQ format to demultiplexed pseudobulk isoform counts, into a single command suitable for execution on either server or cluster. Scywalker includes data quality control, cell type identification, and an interactive report. Assessment of datasets from the human brain, Arabidopsis leaves, and previously benchmarked data from mixed cell lines demonstrate excellent correlation with short-read analyses at both the cell-barcoding and gene quantification levels. At the isoform level, we show that scywalker facilitates the direct identification of cell-type-specific expression of novel isoforms.

AVAILABILITY AND IMPLEMENTATION

Scywalker is available on github.com/derijkp/scywalker under the GNU General Public License (GPL) and at https://zenodo.org/records/13359438/files/scywalker-0.108.0-Linux-x86_64.tar.gz.

A Cell Cycle-Aware Network for Data Integration and Label Transferring of Single-Cell RNA-Seq and ATAC-Seq

In recent years, the integration of single-cell multi-omics data has provided a more comprehensive understanding of cell functions and internal regulatory mechanisms from a non-single omics perspective, but it still suffers many challenges, such as omics-variance, sparsity, cell heterogeneity, and confounding factors. As it is known, the cell cycle is regarded as a confounder when analyzing other factors in single-cell RNA-seq data, but it is not clear how it will work on the integrated single-cell multi-omics data. Here, a cell cycle-aware network (CCAN) is developed to remove cell cycle effects from the integrated single-cell multi-omics data while keeping the cell type-specific variations. This is the first computational model to study the cell-cycle effects in the integration of single-cell multi-omics data. Validations on several benchmark datasets show the outstanding performance of CCAN in a variety of downstream analyses and applications, including removing cell cycle effects and batch effects of scRNA-seq datasets from different protocols, integrating paired and unpaired scRNA-seq and scATAC-seq data, accurately transferring cell type labels from scRNA-seq to scATAC-seq data, and characterizing the differentiation process from hematopoietic stem cells to different lineages in the integration of differentiation data.

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.

Inhibition of the eukaryotic initiation factor-2α kinase PERK decreases risk of autoimmune diabetes in mice

Preventing the onset of autoimmune type 1 diabetes (T1D) is feasible through pharmacological interventions that target molecular stress-responsive mechanisms. Cellular stresses, such as nutrient deficiency, viral infection, or unfolded proteins, trigger the integrated stress response (ISR), which curtails protein synthesis by phosphorylating eukaryotic translation initiation factor-2α (eIF2α). In T1D, maladaptive unfolded protein response (UPR) in insulin-producing β cells renders these cells susceptible to autoimmunity. We found that inhibition of the eIF2α kinase PKR-like ER kinase (PERK), a common component of the UPR and ISR, reversed the mRNA translation block in stressed human islets and delayed the onset of diabetes, reduced islet inflammation, and preserved β cell mass in T1D-susceptible mice. Single-cell RNA-Seq of islets from PERK-inhibited mice showed reductions in the UPR and PERK signaling pathways and alterations in antigen-processing and presentation pathways in β cells. Spatial proteomics of islets from these mice showed an increase in the immune checkpoint protein programmed death-ligand 1 (PD-L1) in β cells. Golgi membrane protein 1, whose levels increased following PERK inhibition in human islets and EndoC-βH1 human β cells, interacted with and stabilized PD-L1. Collectively, our studies show that PERK activity enhances β cell immunogenicity and that inhibition of PERK may offer a strategy for preventing or delaying the development of T1D.

Effective anti-tumor immune response against HCC is orchestrated by immune cell partnership network that functions through hepatic homeostatic pathways, not direct cytotoxicity

The liver harbors a diverse array of immune cells during both health and disease. The specific roles of these cells in nonalcoholic fatty liver disease (NAFLD) and hepatocellular carcinoma (HCC) remain unclear. Using a systems immunology approach, we demonstrate that reciprocal cell-cell communications function through dominant-subdominant pattern of ligand-receptor homeostatic pathways. In the healthy control, hepatocyte-dominated homeostatic pathways induce local immune responses to maintain liver homeostasis. Chronic intake of a Western diet (WD) alters hepatocytes and induces hepatic stellate cell (HSC), cancer cell and NKT cell-dominated interactions during NAFLD. During HCC, monocytes, hepatocytes, and myofibroblasts join the dominant cellular interactions network to restore liver homeostasis. Dietary correction during NAFLD results in nonlinear outcomes with various cellular rearrangements. When cancer cells and stromal cells dominate hepatic interactions network without inducing homeostatic immune responses, HCC progression occurs. Conversely, myofibroblast and fibroblast-dominated network orchestrates monocyte-dominated HCC-preventive immune responses. Tumor immune surveillance by 75% of immune cells successfully promoting liver homeostasis can create a tumor-inhibitory microenvironment, while only 5% of immune cells manifest apoptosis-inducing functions, primarily for facilitating homeostatic liver cell turnover rather than direct tumor killing. These data suggest that an effective immunotherapy should promote liver homeostasis rather than direct tumor killing.

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.

Inhibition of the Eukaryotic Initiation Factor-2-α Kinase PERK Decreases Risk of Autoimmune Diabetes in Mice

Preventing the onset of autoimmune type 1 diabetes (T1D) is feasible through pharmacological interventions that target molecular stress-responsive mechanisms. Cellular stresses, such as nutrient deficiency, viral infection, or unfolded proteins, trigger the integrated stress response (ISR), which curtails protein synthesis by phosphorylating eIF2α. In T1D, maladaptive unfolded protein response (UPR) in insulin-producing β cells renders these cells susceptible to autoimmunity. We show that inhibition of the eIF2α kinase PERK, a common component of the UPR and ISR, reverses the mRNA translation block in stressed human islets and delays the onset of diabetes, reduces islet inflammation, and preserves β cell mass in T1D-susceptible mice. Single-cell RNA sequencing of islets from PERK-inhibited mice shows reductions in the UPR and PERK signaling pathways and alterations in antigen processing and presentation pathways in β cells. Spatial proteomics of islets from these mice shows an increase in the immune checkpoint protein PD-L1 in β cells. Golgi membrane protein 1, whose levels increase following PERK inhibition in human islets and EndoC-βH1 human β cells, interacts with and stabilizes PD-L1. Collectively, our studies show that PERK activity enhances β cell immunogenicity, and inhibition of PERK may offer a strategy to prevent or delay the development of T1D.

WITHDRAWN: LonP1 Drives Proneural Mesenchymal Transition in IDH1-R132H Diffuse Glioma

The authors have withdrawn their manuscript owing to massive revision and data validation. Therefore, the authors do not wish this work to be cited as reference for the project. If you have any questions, please contact the corresponding author.

Modeling type 1 diabetes progression using machine learning and single-cell transcriptomic measurements in human islets

Type 1 diabetes (T1D) is a chronic condition in which beta cells are destroyed by immune cells. Despite progress in immunotherapies that could delay T1D onset, early detection of autoimmunity remains challenging. Here, we evaluate the utility of machine learning for early prediction of T1D using single-cell analysis of islets. Using gradient-boosting algorithms, we model changes in gene expression of single cells from pancreatic tissues in T1D and non-diabetic organ donors. We assess if mathematical modeling could predict the likelihood of T1D development in non-diabetic autoantibody-positive donors. While most autoantibody-positive donors are predicted to be non-diabetic, select donors with unique gene signatures are classified as T1D. Our strategy also reveals a shared gene signature in distinct T1D-associated models across cell types, suggesting a common effect of the disease on transcriptional outputs of these cells. Our study establishes a precedent for using machine learning in early detection of T1D.

SCIPAC: quantitative estimation of cell-phenotype associations

Numerous algorithms have been proposed to identify cell types in single-cell RNA sequencing data, yet a fundamental problem remains: determining associations between cells and phenotypes such as cancer. We develop SCIPAC, the first algorithm that quantitatively estimates the association between each cell in single-cell data and a phenotype. SCIPAC also provides a p-value for each association and applies to data with virtually any type of phenotype. We demonstrate SCIPAC's accuracy in simulated data. On four real cancerous or noncancerous datasets, insights from SCIPAC help interpret the data and generate new hypotheses. SCIPAC requires minimum tuning and is computationally very fast.

CD4+ T Cells From Individuals With Type 1 Diabetes Respond to a Novel Class of Deamidated Peptides Formed in Pancreatic Islets

The β-cell plays a crucial role in the pathogenesis of type 1 diabetes, in part through the posttranslational modification of self-proteins by biochemical processes such as deamidation. These neoantigens are potential triggers for breaking immune tolerance. We report the detection by LC-MS/MS of 16 novel Gln and 27 novel Asn deamidations in 14 disease-related proteins within inflammatory cytokine-stressed human islets of Langerhans. T-cell clones responsive against one Gln- and three Asn-deamidated peptides could be isolated from peripheral blood of individuals with type 1 diabetes. Ex vivo HLA class II tetramer staining detected higher T-cell frequencies in individuals with the disease compared with control individuals. Furthermore, there was a positive correlation between the frequencies of T cells specific for deamidated peptides, insulin antibody levels at diagnosis, and duration of disease. These results highlight that stressed human islets are prone to enzymatic and biochemical deamidation and suggest that both Gln- and Asn-deamidated peptides can promote the activation and expansion of autoreactive CD4+ T cells. These findings add to the growing evidence that posttranslational modifications undermine tolerance and may open the road for the development of new diagnostic and therapeutic applications for individuals living with type 1 diabetes.

ARTICLE HIGHLIGHTS

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.

A comparison of marker gene selection methods for single-cell RNA sequencing data

BACKGROUND

The development of single-cell RNA sequencing (scRNA-seq) has enabled scientists to catalog and probe the transcriptional heterogeneity of individual cells in unprecedented detail. A common step in the analysis of scRNA-seq data is the selection of so-called marker genes, most commonly to enable annotation of the biological cell types present in the sample. In this paper, we benchmark 59 computational methods for selecting marker genes in scRNA-seq data.

RESULTS

We compare the performance of the methods using 14 real scRNA-seq datasets and over 170 additional simulated datasets. Methods are compared on their ability to recover simulated and expert-annotated marker genes, the predictive performance and characteristics of the gene sets they select, their memory usage and speed, and their implementation quality. In addition, various case studies are used to scrutinize the most commonly used methods, highlighting issues and inconsistencies.

CONCLUSIONS

Overall, we present a comprehensive evaluation of methods for selecting marker genes in scRNA-seq data. Our results highlight the efficacy of simple methods, especially the Wilcoxon rank-sum test, Student's t-test, and logistic regression.

A Cell Cycle-aware Network for Data Integration and Label Transferring of Single-cell RNA-seq and ATAC-seq

In recent years, the integration of single-cell multi-omics data has provided a more comprehensive understanding of cell functions and internal regulatory mechanisms from a non-single omics perspective, but it still suffers many challenges, such as omics-variance, sparsity, cell heterogeneity and confounding factors. As we know, cell cycle is regarded as a confounder when analyzing other factors in single-cell RNA-seq data, but it's not clear how it will work on the integrated single-cell multi-omics data. Here, we developed a Cell Cycle-Aware Network (CCAN) to remove cell cycle effects from the integrated single-cell multi-omics data while keeping the cell type-specific variations. This is the first computational model to study the cell-cycle effects in the integration of single-cell multi-omics data. Validations on several benchmark datasets show the out-standing performance of CCAN in a variety of downstream analyses and applications, including removing cell cycle effects and batch effects of scRNA-seq datasets from different protocols, integrating paired and unpaired scRNA-seq and scATAC-seq data, accurately transferring cell type labels from scRNA-seq to scATAC-seq data, and characterizing the differentiation process from hematopoietic stem cells to different lineages in the integration of differentiation data.

SRT-Server: powering the analysis of spatial transcriptomic data

BACKGROUND

Spatial resolved transcriptomics (SRT) encompasses a rapidly developing set of technologies that enable the measurement of gene expression in tissue while retaining spatial localization information. SRT technologies and the enabled SRT studies have provided unprecedent insights into the structural and functional underpinnings of complex tissues. As SRT technologies have advanced and an increasing number of SRT studies have emerged, numerous sophisticated statistical and computational methods have been developed to facilitate the analysis and interpretation of SRT data. However, despite the growing popularity of SRT studies and the widespread availability of SRT analysis methods, analysis of large-scale and complex SRT datasets remains challenging and not easily accessible to researchers with limited statistical and computational backgrounds.

RESULTS

Here, we present SRT-Server, the first webserver designed to carry out comprehensive SRT analyses for a wide variety of SRT technologies while requiring minimal prior computational knowledge. Implemented with cutting-edge web development technologies, SRT-Server is user-friendly and features multiple analytic modules that can perform a range of SRT analyses. With a flowchart-style interface, these different analytic modules on the SRT-Server can be dragged into the main panel and connected to each other to create custom analytic pipelines. SRT-Server then automatically executes the desired analyses, generates corresponding figures, and outputs results-all without requiring prior programming knowledge. We demonstrate the advantages of SRT-Server through three case studies utilizing SRT data collected from two common platforms, highlighting its versatility and values to researchers with varying analytic expertise.

CONCLUSIONS

Overall, SRT-Server presents a user-friendly, efficient, effective, secure, and expandable solution for SRT data analysis, opening new doors for researchers in the field. SRT-Server is freely available at https://spatialtranscriptomicsanalysis.com/ .

Cell-type brain-region specific changes in prefrontal cortex of a mouse model of alcohol dependence

The prefrontal cortex is a crucial regulator of alcohol drinking, and dependence, and other behavioral phenotypes associated with AUD. Comprehensive identification of cell-type specific transcriptomic changes in alcohol dependence will improve our understanding of mechanisms underlying the excessive alcohol use associated with alcohol dependence and will refine targets for therapeutic development. We performed single nucleus RNA sequencing (snRNA-seq) and Visium spatial gene expression profiling on the medial prefrontal cortex (mPFC) obtained from C57BL/6 J mice exposed to the two-bottle choice-chronic intermittent ethanol (CIE) vapor exposure (2BC-CIE, defined as dependent group) paradigm which models phenotypes of alcohol dependence including escalation of alcohol drinking. Gene co-expression network analysis and differential expression analysis identified highly dysregulated co-expression networks in multiple cell types. Dysregulated modules and their hub genes suggest novel understudied targets for studying molecular mechanisms contributing to the alcohol dependence state. A subtype of inhibitory neurons was the most alcohol-sensitive cell type and contained a downregulated gene co-expression module; the hub gene for this module is Cpa6, a gene previously identified by GWAS to be associated with excessive alcohol consumption. We identified an astrocytic Gpc5 module significantly upregulated in the alcohol-dependent group. To our knowledge, there are no studies linking Cpa6 and Gpc5 to the alcohol-dependent phenotype. We also identified neuroinflammation related gene expression changes in multiple cell types, specifically enriched in microglia, further implicating neuroinflammation in the escalation of alcohol drinking. Here, we present a comprehensive atlas of cell-type specific alcohol dependence mediated gene expression changes in the mPFC and identify novel cell type-specific targets implicated in alcohol dependence.

Pathway centric analysis for single-cell RNA-seq and spatial transcriptomics data with GSDensity

Advances in single-cell technology have enabled molecular dissection of heterogeneous biospecimens at unprecedented scales and resolutions. Cluster-centric approaches are widely applied in analyzing single-cell data, however they have limited power in dissecting and interpreting highly heterogenous, dynamically evolving data. Here, we present GSDensity, a graph-modeling approach that allows users to obtain pathway-centric interpretation and dissection of single-cell and spatial transcriptomics (ST) data without performing clustering. Using pathway gene sets, we show that GSDensity can accurately detect biologically distinct cells and reveal novel cell-pathway associations ignored by existing methods. Moreover, GSDensity, combined with trajectory analysis can identify curated pathways that are active at various stages of mouse brain development. Finally, GSDensity can identify spatially relevant pathways in mouse brains and human tumors including those following high-order organizational patterns in the ST data. Particularly, we create a pan-cancer ST map revealing spatially relevant and recurrently active pathways across six different tumor types.

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.

Probabilistic cell/domain-type assignment of spatial transcriptomics data with SpatialAnno

In the analysis of both single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) data, classifying cells/spots into cell/domain types is an essential analytic step for many secondary analyses. Most of the existing annotation methods have been developed for scRNA-seq datasets without any consideration of spatial information. Here, we present SpatialAnno, an efficient and accurate annotation method for spatial transcriptomics datasets, with the capability to effectively leverage a large number of non-marker genes as well as 'qualitative' information about marker genes without using a reference dataset. Uniquely, SpatialAnno estimates low-dimensional embeddings for a large number of non-marker genes via a factor model while promoting spatial smoothness among neighboring spots via a Potts model. Using both simulated and four real spatial transcriptomics datasets from the 10x Visium, ST, Slide-seqV1/2, and seqFISH platforms, we showcase the method's improved spatial annotation accuracy, including its robustness to the inclusion of marker genes for irrelevant cell/domain types and to various degrees of marker gene misspecification. SpatialAnno is computationally scalable and applicable to SRT datasets from different platforms. Furthermore, the estimated embeddings for cellular biological effects facilitate many downstream analyses.

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.

Protection of β cells against pro-inflammatory cytokine stress by the GDF15-ERBB2 signaling

Aim/hypothesis

Growth/differentiation factor 15 (GDF15) is a therapeutic target for a variety of metabolic diseases, including type 1 diabetes (T1D). However, the nausea caused by GDF15 is a challenging point for therapeutic development. In addition, it is unknown why the endogenous GDF15 fails to protect from T1D development. Here, we investigate the GDF15 signaling in pancreatic islets towards opening possibilities for therapeutic targeting in β cells and to understand why this protection fails to occur naturally.

Methods

GDF15 signaling in islets was determined by proximity-ligation assay, untargeted proteomics, pathway analysis, and treatment of cells with specific inhibitors. To determine if GDF15 levels would increase prior to disease onset, plasma levels of GDF15 were measured in a longitudinal prospective study of children during T1D development (n=132 cases vs. n=40 controls) and in children with islet autoimmunity but normoglycemia (n=47 cases vs. n=40 controls) using targeted mass spectrometry. We also investigated the regulation of GDF15 production in islets by fluorescence microscopy and western blot analysis.

Results

The proximity-ligation assay identified ERBB2 as the GDF15 receptor in islets, which was confirmed using its specific antagonist, tucatinib. The untargeted proteomics analysis and caspase assay showed that ERBB2 activation by GDF15 reduces β cell apoptosis by downregulating caspase 8. In plasma, GDF15 levels were higher (p=0.0024) during T1D development compared to controls, but not in islet autoimmunity with normoglycemia. However, in the pancreatic islets GDF15 was depleted via sequestration of its mRNA into stress granules, resulting in translation halting.

Conclusions/interpretation

GDF15 protects against T1D via ERBB2-mediated decrease of caspase 8 expression in pancreatic islets. Circulating levels of GDF15 increases pre-T1D onset, which is insufficient to promote protection due to its localized depletion in the islets. These findings open opportunities for targeting GDF15 downstream signaling for pancreatic β cell protection in T1D and help to explain the lack of natural protection by the endogenous protein.

Annotation of cell types (ACT): a convenient web server for cell type annotation

BACKGROUND

The advancement of single-cell sequencing has progressed our ability to solve biological questions. Cell type annotation is of vital importance to this process, allowing for the analysis and interpretation of enormous single-cell datasets. At present, however, manual cell annotation which is the predominant approach remains limited by both speed and the requirement of expert knowledge.

METHODS

To address these challenges, we constructed a hierarchically organized marker map through manually curating over 26,000 cell marker entries from about 7000 publications. We then developed WISE, a weighted and integrated gene set enrichment method, to integrate the prevalence of canonical markers and ordered differentially expressed genes of specific cell types in the marker map. Benchmarking analysis suggested that our method outperformed state-of-the-art methods.

RESULTS

By integrating the marker map and WISE, we developed a user-friendly and convenient web server, ACT ( http://xteam.xbio.top/ACT/ or http://biocc.hrbmu.edu.cn/ACT/ ), which only takes a simple list of upregulated genes as input and provides interactive hierarchy maps, together with well-designed charts and statistical information, to accelerate the assignment of cell identities and made the results comparable to expert manual annotation. Besides, a pan-tissue marker map was constructed to assist in cell assignments in less-studied tissues. Applying ACT to three case studies showed that all cell clusters were quickly and accurately annotated, and multi-level and more refined cell types were identified.

CONCLUSIONS

We developed a knowledge-based resource and a corresponding method, together with an intuitive graphical web interface, for cell type annotation. We believe that ACT, emerging as a powerful tool for cell type annotation, would be widely used in single-cell research and considerably accelerate the process of cell type identification.

Widespread sex dimorphism across single-cell transcriptomes of adult African turquoise killifish tissues

The African turquoise killifish (Nothobranchius furzeri), the shortest-lived vertebrate that can be bred in captivity, is an emerging model organism for aging research. Here, we describe a multitissue, single-cell gene expression atlas of female and male blood, kidney, liver, and spleen. We annotate 22 cell types, define marker genes, and infer differentiation trajectories. We find pervasive sex-dimorphic gene expression across cell types. Sex-dimorphic genes tend to be linked to lipid metabolism, consistent with clear differences in lipid storage in female vs. male turquoise killifish livers. We use machine learning to predict sex using single-cell gene expression and identify potential markers for molecular sex identity. As a proof of principle, we show that our atlas can be used to deconvolute existing bulk RNA sequencing (RNA-seq) data to obtain accurate estimates of cell type proportions. This atlas can be a resource to the community that could be leveraged to develop cell-type-specific expression in transgenic animals.

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.

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.

The crosstalking immune cells network creates a collective function beyond the function of each cellular constituent during the progression of hepatocellular carcinoma

Abundance of data on the role of inflammatory immune responses in the progression or inhibition of hepatocellular carcinoma (HCC) has failed to offer a curative immunotherapy for HCC. This is largely because of focusing on detailed specific cell types and missing the collective function of the hepatic immune system. To discover the collective immune function, we take systems immunology approach by performing high-throughput analysis of snRNAseq data collected from the liver of DIAMOND mice during the progression of nonalcoholic fatty liver disease (NAFLD) to HCC. We report that mutual signaling interactions of the hepatic immune cells in a dominant-subdominant manner, as well as their interaction with structural cells shape the immunological pattern manifesting a collective function beyond the function of the cellular constituents. Such pattern discovery approach recognized direct role of the innate immune cells in the progression of NASH and HCC. These data suggest that discovery of the immune pattern not only detects the immunological mechanism of HCC in spite of dynamic changes in immune cells during the course of disease but also offers immune modulatory interventions for the treatment of NAFLD and HCC.

Identification of ulcerative colitis-specific immune cell signatures from public single-cell RNA-seq data

BACKGROUND

Single-cell RNA-seq enabled microscopic studies on tissue microenvironment of many diseases. Inflammatory bowel disease, an autoimmune disease, is involved with various dysfunction of immune cells, for which single-cell RNA-seq may provide us a deeper insight into the causes and mechanism of this complex disease.

OBJECTIVE

In this work, we used public single-cell RNA-seq data to study tissue microenvironment around ulcerative colitis, an inflammatory bowel disease causing chronic inflammation and ulcers in large intestine.

METHODS

Since not all the datasets provide cell-type annotations, we first identified cell identities to select cell populations of our interest. Differentially expressed genes and gene set enrichment analysis was then performed to infer the polarization/activation state of macrophages and T cells. Cell-to-cell interaction analysis was also performed to discover distinct interactions in ulcerative colitis.

RESULTS

Differentially expressed genes analysis of the two datasets confirmed the regulation of CTLA4, IL2RA, and CCL5 genes in the T cell subset and regulation of S100A8/A9, CLEC10A genes in macrophages. Cell-to-cell interaction analysis showed CD4+ T cells and macrophages interact actively to each other. We also identified IL-18 pathway activation in inflammatory macrophages, evidence that CD4+ T cells induce Th1 and Th2 differentiation, and also found that macrophages regulate T cell activation through different ligand-receptor pairs, viz. CD86-CTL4, LGALS9-CD47, SIRPA-CD47, and GRN-TNFRSF1B.

CONCLUSION

Analysis of these immune cell subsets may suggest novel strategies for the treatment of inflammatory bowel disease.

Consensus label propagation with graph convolutional networks for single-cell RNA sequencing cell type annotation

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) data, annotated by cell type, is useful in a variety of downstream biological applications, such as profiling gene expression at the single-cell level. However, manually assigning these annotations with known marker genes is both time-consuming and subjective.

RESULTS

We present a Graph Convolutional Network (GCN)-based approach to automate the annotation process. Our process builds upon existing labeling approaches, using state-of-the-art tools to find cells with highly confident label assignments through consensus and spreading these confident labels with a semi-supervised GCN. Using simulated data and two scRNA-seq datasets from different tissues, we show that our method improves accuracy over a simple consensus algorithm and the average of the underlying tools. We also compare our method to a nonparametric neighbor majority approach, showing comparable results. We then demonstrate that our GCN method allows for feature interpretation, identifying important genes for cell type classification. We present our completed pipeline, written in PyTorch, as an end-to-end tool for automating and interpreting the classification of scRNA-seq data.

AVAILABILITY AND IMPLEMENTATION

Our code for conducting the experiments in this paper and using our model is available at https://github.com/lewinsohndp/scSHARP.

Brain matters: unveiling the distinct contributions of region, age, and sex to glia diversity and CNS function

The myelinated white matter tracts of the central nervous system (CNS) are essential for fast transmission of electrical impulses and are often differentially affected in human neurodegenerative diseases across CNS region, age and sex. We hypothesize that this selective vulnerability is underpinned by physiological variation in white matter glia. Using single nucleus RNA sequencing of human post-mortem white matter samples from the brain, cerebellum and spinal cord and subsequent tissue-based validation we found substantial glial heterogeneity with tissue region: we identified region-specific oligodendrocyte precursor cells (OPCs) that retain developmental origin markers into adulthood, distinguishing them from mouse OPCs. Region-specific OPCs give rise to similar oligodendrocyte populations, however spinal cord oligodendrocytes exhibit markers such as SKAP2 which are associated with increased myelin production and we found a spinal cord selective population particularly equipped for producing long and thick myelin sheaths based on the expression of genes/proteins such as HCN2. Spinal cord microglia exhibit a more activated phenotype compared to brain microglia, suggesting that the spinal cord is a more pro-inflammatory environment, a difference that intensifies with age. Astrocyte gene expression correlates strongly with CNS region, however, astrocytes do not show a more activated state with region or age. Across all glia, sex differences are subtle but the consistent increased expression of protein-folding genes in male donors hints at pathways that may contribute to sex differences in disease susceptibility. These findings are essential to consider for understanding selective CNS pathologies and developing tailored therapeutic strategies.

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.

CIForm as a Transformer-based model for cell-type annotation of large-scale single-cell RNA-seq data

Single-cell omics technologies have made it possible to analyze the individual cells within a biological sample, providing a more detailed understanding of biological systems. Accurately determining the cell type of each cell is a crucial goal in single-cell RNA-seq (scRNA-seq) analysis. Apart from overcoming the batch effects arising from various factors, single-cell annotation methods also face the challenge of effectively processing large-scale datasets. With the availability of an increase in the scRNA-seq datasets, integrating multiple datasets and addressing batch effects originating from diverse sources are also challenges in cell-type annotation. In this work, to overcome the challenges, we developed a supervised method called CIForm based on the Transformer for cell-type annotation of large-scale scRNA-seq data. To assess the effectiveness and robustness of CIForm, we have compared it with some leading tools on benchmark datasets. Through the systematic comparisons under various cell-type annotation scenarios, we exhibit that the effectiveness of CIForm is particularly pronounced in cell-type annotation. The source code and data are available at https://github.com/zhanglab-wbgcas/CIForm.

Celloscope: a probabilistic model for marker-gene-driven cell type deconvolution in spatial transcriptomics data

Spatial transcriptomics maps gene expression across tissues, posing the challenge of determining the spatial arrangement of different cell types. However, spatial transcriptomics spots contain multiple cells. Therefore, the observed signal comes from mixtures of cells of different types. Here, we propose an innovative probabilistic model, Celloscope, that utilizes established prior knowledge on marker genes for cell type deconvolution from spatial transcriptomics data. Celloscope outperforms other methods on simulated data, successfully indicates known brain structures and spatially distinguishes between inhibitory and excitatory neuron types based in mouse brain tissue, and dissects large heterogeneity of immune infiltrate composition in prostate gland tissue.

Widespread sex-dimorphism across single-cell transcriptomes of adult African turquoise killifish tissues

The African turquoise killifish (Nothobranchius furzeri), the shortest-lived vertebrate that can be bred in captivity, is an emerging model organism to study vertebrate aging. Here we describe the first multi-tissue, single-cell gene expression atlas of female and male turquoise killifish tissues comprising immune and metabolic cells from the blood, kidney, liver, and spleen. We were able to annotate 22 distinct cell types, define associated marker genes, and infer differentiation trajectories. Using this dataset, we found pervasive sex-dimorphic gene expression across cell types, especially in the liver. Sex-dimorphic genes tended to be involved in processes related to lipid metabolism, and indeed, we observed clear differences in lipid storage in female vs. male turquoise killifish livers. Importantly, we use machine-learning to predict sex using single-cell gene expression in our atlas and identify potential transcriptional markers for molecular sex identity in this species. As proof-of-principle, we show that our atlas can be used to deconvolute existing liver bulk RNA-seq data in this species to obtain accurate estimates of cell type proportions across biological conditions. We believe that this single-cell atlas can be a resource to the community that could notably be leveraged to identify cell type-specific genes for cell type-specific expression in transgenic animals.

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.

Hierarchical cell-type identifier accurately distinguishes immune-cell subtypes enabling precise profiling of tissue microenvironment with single-cell RNA-sequencing

Single-cell RNA-seq enabled in-depth study on tissue micro-environment and immune-profiling, where a crucial step is to annotate cell identity. Immune cells play key roles in many diseases, whereas their activities are hard to track due to their diverse and highly variable nature. Existing cell-type identifiers had limited performance for this purpose. We present HiCAT, a hierarchical, marker-based cell-type identifier utilising gene set analysis for statistical scoring for given markers. It features successive identification of major-type, minor-type and subsets utilising subset markers structured in a three-level taxonomy tree. Comparison with manual annotation and pairwise match test showed HiCAT outperforms others in major- and minor-type identification. For subsets, we qualitatively evaluated the marker expression profile demonstrating that HiCAT provide the clearest immune-cell landscape. HiCAT was also used for immune-cell profiling in ulcerative colitis and discovered distinct features of the disease in macrophage and T-cell subsets that could not be identified previously.

Role of lymphoid lineage cells aberrantly expressing alarmins S100A8/A9 in determining the severity of COVID-19

BACKGROUND

Alarmins S100A8 and S100A9 are recognized as hallmarks of severe COVID-19 and are primarily produced in myeloid cells, such as monocytes and neutrophils. As single-cell RNA-sequencing (scRNA-seq) data from patients with COVID-19 revealed the expression of S100A8/A9 in lymphoid cells in patients with severe COVID-19.

OBJECTIVE

We investigated the characteristics of lymphoid cells expressing S100A8/A9 in COVID-19 patients.

METHODS

Publicly available scRNA-seq data from patients with mild (N = 12) or severe (N = 7) COVID-19 were reanalyzed. The data were further divided into the following two groups based on the time of sample collection (from infection-onset): within 6 days (early phase) and after 6 days (late phase). Differential expression and gene set enrichment analyses were performed between S100A8/A9^High and S100A8/A9^Low lymphoid cells. Finally, cell-cell interaction analysis was performed to investigate the role of lymphoid cells expressing high levels of S100A8/A9 in COVID-19.

RESULTS

S100A8/A9 overexpression was observed in lymphoid cells, including B cells, T cells, and NK cells, in patients with severe COVID-19 (compared to patients with mild COVID-19). Cells exhibiting strong interferon/cytokine responses were found to be associated with the severity of COVID-19. Furthermore, differences in S100A8/A9-TLR4/RAGE interactions were confirmed between patients with severe and mild disease.

CONCLUSIONS

Lymphoid cells overexpressing S100A8/A9 contribute to the dysregulation of the innate immune response in patients with severe COVID-19, specifically during the early phase of infection. This study fosters a better understanding of the hyper-induction of pro-inflammatory cytokine expression and the generation of a cytokine storm in response to COVID-19 infection.

Complex Analysis of Single-Cell RNA Sequencing Data

Single-cell RNA sequencing (scRNA-seq) is a revolutionary tool for studying the physiology of normal and pathologically altered tissues. This approach provides information about molecular features (gene expression, mutations, chromatin accessibility, etc.) of cells, opens up the possibility to analyze the trajectories/phylogeny of cell differentiation and cell-cell interactions, and helps in discovery of new cell types and previously unexplored processes. From a clinical point of view, scRNA-seq facilitates deeper and more detailed analysis of molecular mechanisms of diseases and serves as a basis for the development of new preventive, diagnostic, and therapeutic strategies. The review describes different approaches to the analysis of scRNA-seq data, discusses the advantages and disadvantages of bioinformatics tools, provides recommendations and examples of their successful use, and suggests potential directions for improvement. We also emphasize the need for creating new protocols, including multiomics ones, for the preparation of DNA/RNA libraries of single cells with the purpose of more complete understanding of individual cells.

ImmCluster: an ensemble resource for immunology cell type clustering and annotations in normal and cancerous tissues

Single-cell transcriptome has enabled the transcriptional profiling of thousands of immune cells in complex tissues and cancers. However, subtle transcriptomic differences in immune cell subpopulations and the high dimensionality of transcriptomic data make the clustering and annotation of immune cells challenging. Herein, we introduce ImmCluster (http://bio-bigdata.hrbmu.edu.cn/ImmCluster) for immunology cell type clustering and annotation. We manually curated 346 well-known marker genes from 1163 studies. ImmCluster integrates over 420 000 immune cells from nine healthy tissues and over 648 000 cells from different tumour samples of 17 cancer types to generate stable marker-gene sets and develop context-specific immunology references. In addition, ImmCluster provides cell clustering using seven reference-based and four marker gene-based computational methods, and the ensemble method was developed to provide consistent cell clustering than individual methods. Five major analytic modules were provided for interactively exploring the annotations of immune cells, including clustering and annotating immune cell clusters, gene expression of markers, functional assignment in cancer hallmarks, cell states and immune pathways, cell-cell communications and the corresponding ligand-receptor interactions, as well as online tools. ImmCluster generates diverse plots and tables, enabling users to identify significant associations in immune cell clusters simultaneously. ImmCluster is a valuable resource for analysing cellular heterogeneity in cancer microenvironments.

Computational workflow and interactive analysis of single-cell expression profiling of islets generated by the Human Pancreas Analysis Program

Type 1 and Type 2 diabetes are distinct genetic diseases of the pancreas which are defined by the abnormal level of blood glucose. Understanding the initial molecular perturbations that occur during the pathogenesis of diabetes is of critical importance in understanding these disorders. The inability to biopsy the human pancreas of living donors hampers insights into early detection, as the majority of diabetes studies have been performed on peripheral leukocytes from the blood, which is not the site of pathogenesis. Therefore, efforts have been made by various teams including the Human Pancreas Analysis Program (HPAP) to collect pancreatic tissues from deceased organ donors with different clinical phenotypes. HPAP is designed to define the molecular pathogenesis of islet dysfunction by generating detailed datasets of functional, cellular, and molecular information in pancreatic tissues of clinically well-defined organ donors with Type 1 and Type 2 diabetes. Moreover, data generated by HPAP continously become available through a centralized database, PANC-DB, thus enabling the diabetes research community to access these multi-dimensional data prepublication. Here, we present the computational workflow for single-cell RNA-seq data analysis of 258,379 high-quality cells from the pancreatic islets of 67 human donors generated by HPAP, the largest existing scRNA-seq dataset of human pancreatic tissues. We report various computational steps including preprocessing, doublet removal, clustering and cell type annotation across single-cell RNA-seq data from islets of four distintct classes of organ donors, i.e. non-diabetic control, autoantibody positive but normoglycemic, Type 1 diabetic, and Type 2 diabetic individuals. Moreover, we present an interactive tool, called CellxGene developed by the Chan Zuckerberg initiative, to navigate these high-dimensional datasets. Our data and interactive tools provide a reliable reference for singlecell pancreatic islet biology studies, especially diabetes-related conditions.

scWECTA: A weighted ensemble classification framework for cell type assignment based on single cell transcriptome

Rapid advances in single-cell transcriptome analysis provide deeper insights into the study of tissue heterogeneity at the cellular level. Unsupervised clustering can identify potential cell populations in single-cell RNA-sequencing (scRNA-seq) data, but fail to further determine the identity of each cell. Existing automatic annotation methods using scRNA-seq data based on machine learning mainly use single feature set and single classifier. In view of this, we propose a Weighted Ensemble classification framework for Cell Type Annotation, named scWECTA, which improves the accuracy of cell type identification. scWECTA uses five informative gene sets and integrates five classifiers based on soft weighted ensemble framework. And the ensemble weights are inferred through the constrained non-negative least squares. Validated on multiple pairs of scRNA-seq datasets, scWECTA is able to accurately annotate scRNA-seq data across platforms and across tissues, especially for imbalanced data containing rare cell types. Moreover, scWECTA outperforms other comparable methods in balancing the prediction accuracy of common cell types and the unassigned rate of non-common cell types at the same time. The source code of scWECTA is freely available at https://github.com/ttren-sc/scWECTA.

Data analysis guidelines for single-cell RNA-seq in biomedical studies and clinical applications

The application of single-cell RNA sequencing (scRNA-seq) in biomedical research has advanced our understanding of the pathogenesis of disease and provided valuable insights into new diagnostic and therapeutic strategies. With the expansion of capacity for high-throughput scRNA-seq, including clinical samples, the analysis of these huge volumes of data has become a daunting prospect for researchers entering this field. Here, we review the workflow for typical scRNA-seq data analysis, covering raw data processing and quality control, basic data analysis applicable for almost all scRNA-seq data sets, and advanced data analysis that should be tailored to specific scientific questions. While summarizing the current methods for each analysis step, we also provide an online repository of software and wrapped-up scripts to support the implementation. Recommendations and caveats are pointed out for some specific analysis tasks and approaches. We hope this resource will be helpful to researchers engaging with scRNA-seq, in particular for emerging clinical applications.

Biologic insights from single-cell studies of psoriasis and psoriatic arthritis

INTRODUCTION

Psoriasis (PSO) and psoriatic arthritis (PSA) represent a large burden of global inflammatory disease, but sustained treatment response and early diagnosis remain challenging. Both conditions arise from complex immune cell dysregulation. Single-cell techniques, including single-cell RNA sequencing (scRNA-seq), have revolutionized our understanding of pathogenesis by illuminating heterogeneous cell populations and their interactions.

AREAS COVERED

We discuss the transcriptional profiles and cellular interactions unique to PSO/PSA affecting T cells, myeloid cells, keratinocytes, innate lymphoid cells, and stromal cells. We also review advances, limitations, and future challenges associated with single-cell studies.

EXPERT OPINION

Following analyses of 22 single-cell studies, several themes emerged. A small subpopulation of cells can have a large impact on disease pathogenesis. Multiple cell types identified via scRNA-seq play supporting roles in PSO pathogenesis, contrary to the traditional paradigm focusing on IL-23/IL-17 signaling among dendritic cells and T cells. Immune cell states are dynamic, with psoriatic subpopulations aberrantly re-activating and differentiating into inflammatory phenotypes depending on surrounding signaling cues. Comparison of circulating immune cells with resident skin/joint cells has uncovered specific T cell clonotypes associated with the disease. Finally, machine learning models demonstrate great promise in identifying biomarkers to diagnose clinically ambiguous rashes and PSA at earlier stages.

Fully-automated and ultra-fast cell-type identification using specific marker combinations from single-cell transcriptomic data

Identification of cell populations often relies on manual annotation of cell clusters using established marker genes. However, the selection of marker genes is a time-consuming process that may lead to sub-optimal annotations as the markers must be informative of both the individual cell clusters and various cell types present in the sample. Here, we developed a computational platform, ScType, which enables a fully-automated and ultra-fast cell-type identification based solely on a given scRNA-seq data, along with a comprehensive cell marker database as background information. Using six scRNA-seq datasets from various human and mouse tissues, we show how ScType provides unbiased and accurate cell type annotations by guaranteeing the specificity of positive and negative marker genes across cell clusters and cell types. We also demonstrate how ScType distinguishes between healthy and malignant cell populations, based on single-cell calling of single-nucleotide variants, making it a versatile tool for anticancer applications. The widely applicable method is deployed both as an interactive web-tool (https://sctype.app), and as an open-source R-package.

Cell types are typically identified in single cell transcriptomic data by manual annotation of cell clusters using established marker genes. Here the authors present a fully-automated computational platform that can quickly and accurately distinguish between cell types.