Exploring the single-cell RNA-seq analysis landscape with the scRNA-tools database

Single-cell multiomics to advance cell therapy

Single-cell RNA-sequencing (scRNAseq) was first introduced in 2009 and has evolved with many technological advancements over the last decade. Not only are there several scRNAseq platforms differing in many aspects, but there are also a large number of computational pipelines available for downstream analyses which are being developed at an exponential rate. Such computational data appear in many scientific publications in virtually every field of study; thus, investigators should be able to understand and interpret data in this rapidly evolving field. Here, we discuss key differences in scRNAseq platforms, crucial steps in scRNAseq experiments, standard downstream analyses and introduce newly developed multimodal approaches. We then discuss how single-cell omics has been applied to advance the field of cell therapy.

Applications and emerging challenges of single-cell RNA sequencing technology in tumor drug discovery

Current therapeutic drugs are inadequate for curing tumors, highlighting the need for novel tumor drugs. The advancement of single-cell RNA sequencing (scRNA-seq) technology offers new opportunities for tumor drug discovery. This technology allows us to explore tumor heterogeneity and developmental mechanisms at the single-cell level. In this review, we outline the application of scRNA-seq in tumor drug discovery stages, including elucidating tumor mechanisms, identifying targets, screening drugs, and understanding drug action and resistance. We also discuss the challenges and future prospects of using scRNA-seq in drug development, providing a scientific foundation for advancing tumor therapies.

Application of bioinformatic tools in cell type classification for single-cell RNA-seq data

The advancements in single-cell RNA sequencing (scRNAseq) technology have significantly transformed genomics research, enabling the handling of thousands of cells in each experiment. As of now, 32,068 research studies have been cataloged in the Pubmed database. The primary aim of scRNAseq investigations is to identify cell types, understand the antitumor immune response, and identify new and uncommon cell types. Traditional techniques for identifying cell types include microscopy, histology, and pathological characteristics. However, the complexity of instruments and the need for precise experimental design make it difficult to fully capture the overall heterogeneity. Unsupervised clustering and supervised classification methods have been used to solve this task. Supervised cell type classification methods have gained popularity as large-scale, high-quality, well-annotated and more robust results compared to clustering methods. A recent study showed that support vector machine (SVM) gives a high-quality classification performance in different scenarios. In this article, we compare and evaluate the performance of four different kernels (sigmoid, linear, radial, polynomial) of SVM. The results of the experiments on three standard scRNA-seq datasets indicate that SVM with linear and SVM with sigmoid kernel classify the cells more accurately (approx. 99 %) where SVM linear kernel method has remarkably fast computation time and we also evaluate the results using some single cell specific evaluation matrices F-1 score, MCC, AUC value. Additionally, it sheds light on the potential use of kernels of SVM to give underlying information of single-cell RNA-Seq data more effectively.

Vitessce: integrative visualization of multimodal and spatially resolved single-cell data

Multiomics technologies with single-cell and spatial resolution make it possible to measure thousands of features across millions of cells. However, visual analysis of high-dimensional transcriptomic, proteomic, genome-mapped and imaging data types simultaneously remains a challenge. Here we describe Vitessce, an interactive web-based visualization framework for exploration of multimodal and spatially resolved single-cell data. We demonstrate integrative visualization of millions of data points, including cell-type annotations, gene expression quantities, spatially resolved transcripts and cell segmentations, across multiple coordinated views. The open-source software is available at http://vitessce.io .

Spatial transcriptomics in breast cancer: providing insight into tumor heterogeneity and promoting individualized therapy

A comprehensive understanding of tumor heterogeneity, tumor microenvironment and the mechanisms of drug resistance is fundamental to advancing breast cancer research. While single-cell RNA sequencing has resolved the issue of "temporal dynamic expression" of genes at the single-cell level, the lack of spatial information still prevents us from gaining a comprehensive understanding of breast cancer. The introduction and application of spatial transcriptomics addresses this limitation. As the annual technical method of 2020, spatial transcriptomics preserves the spatial location of tissues and resolves RNA-seq data to help localize and differentiate the active expression of functional genes within a specific tissue region, enabling the study of spatial location attributes of gene locations and cellular tissue environments. In the context of breast cancer, spatial transcriptomics can assist in the identification of novel breast cancer subtypes and spatially discriminative features that show promise for individualized precise treatment. This article summarized the key technical approaches, recent advances in spatial transcriptomics and its applications in breast cancer, and discusses the limitations of current spatial transcriptomics methods and the prospects for future development, with a view to advancing the application of this technology in clinical practice.

Unveiling Leydig cell heterogeneity and its role in male infertility: A single-cell transcriptomic study of human testicular tissue

Male infertility and impaired spermatogenesis are significant concerns in reproductive health, often linked to disruptions in the cellular and molecular processes within the testis. The cellular composition and transcriptional dynamics of human testicular tissue are crucial for understanding these issues. Previous studies have largely relied on bulk tissue analysis, which obscures the distinct roles and interactions of specific cell types. Here, through a comprehensive single-cell transcriptomic analysis of human testes across various developmental stages and pathological conditions, we reveal the intricate cellular heterogeneity and the molecular mechanisms underlying testicular function. Our study identifies significant disruptions in the differentiation trajectories of Germ cells in conditions such as Klinefelter syndrome (KS), AZFa deletion, and Sertoli-cell-only syndrome (SCOS). We further uncover key transcription factors and regulatory networks governing Leydig cell function, particularly those related to steroidogenesis and hormonal regulation. These findings highlight the organized yet complex cellular and molecular landscape of the testis and uncover critical pathways altered in male infertility. Collectively, our data suggest that targeted therapeutic strategies could be developed to address specific disruptions in testicular cell populations and their associated regulatory networks.

Dynamic immunoediting by macrophages in homologous recombination deficiency-stratified pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is a lethal disease, notably resistant to existing therapies. Current research indicates that PDAC patients deficient in homologous recombination (HR) benefit from platinum-based treatments and poly-ADP-ribose polymerase inhibitors (PARPi). However, the effectiveness of PARPi in HR-deficient (HRD) PDAC is suboptimal, and significant challenges remain in fully understanding the distinct characteristics and implications of HRD-associated PDAC. We analyzed 16 PDAC patient-derived tissues, categorized by their homologous recombination deficiency (HRD) scores, and performed high-plex immunofluorescence analysis to define 20 cell phenotypes, thereby generating an in-situ PDAC tumor-immune landscape. Spatial phenotypic-transcriptomic profiling guided by regions-of-interest (ROIs) identified a crucial regulatory mechanism through localized tumor-adjacent macrophages, potentially in an HRD-dependent manner. Cellular neighborhood (CN) analysis further demonstrated the existence of macrophage-associated high-ordered cellular functional units in spatial contexts. Using our multi-omics spatial profiling strategy, we uncovered a dynamic macrophage-mediated regulatory axis linking HRD status with SIGLEC10 and CD52. These findings demonstrate the potential of targeting CD52 in combination with PARPi as a therapeutic intervention for PDAC.

Single-Cell RNA-Sequencing: Opening New Horizons for Breast Cancer Research

Breast cancer is the most prevalent malignant tumor among women with high heterogeneity. Traditional techniques frequently struggle to comprehensively capture the intricacy and variety of cellular states and interactions within breast cancer. As global precision medicine rapidly advances, single-cell RNA sequencing (scRNA-seq) has become a highly effective technique, revolutionizing breast cancer research by offering unprecedented insights into the cellular heterogeneity and complexity of breast cancer. This cutting-edge technology facilitates the analysis of gene expression profiles at the single-cell level, uncovering diverse cell types and states within the tumor microenvironment. By dissecting the cellular composition and transcriptional signatures of breast cancer cells, scRNA-seq provides new perspectives for understanding the mechanisms behind tumor therapy, drug resistance and metastasis in breast cancer. In this review, we summarized the working principle and workflow of scRNA-seq and emphasized the major applications and discoveries of scRNA-seq in breast cancer research, highlighting its impact on our comprehension of breast cancer biology and its potential for guiding personalized treatment strategies.

Single-Cell RNA Sequencing and Its Applications in Pituitary Research

BACKGROUND

Mounting evidence underscores the significance of cellular diversity within the endocrine system and the intricate interplay between different cell types and tissues, essential for preserving physiological balance and influencing disease trajectories. The pituitary gland, a central player in the endocrine orchestra, exemplifies this complexity with its assortment of hormone-secreting and nonsecreting cells.

SUMMARY

The pituitary gland houses several types of cells responsible for hormone production, alongside nonsecretory cells like fibroblasts and endothelial cells, each playing a crucial role in the gland's function and regulatory mechanisms. Despite the acknowledged importance of these cellular interactions, the detailed mechanisms by which they contribute to pituitary gland physiology and pathology remain largely uncharted. The last decade has seen the emergence of groundbreaking technologies such as single-cell RNA sequencing, offering unprecedented insights into cellular heterogeneity and interactions. However, the application of this advanced tool in exploring the pituitary gland's complexities has been scant. This review provides an overview of this methodology, highlighting its strengths and limitations, and discusses future possibilities for employing it to deepen our understanding of the pituitary gland and its dysfunction in disease states.

KEY MESSAGE

Single-cell RNA sequencing technology offers an unprecedented means to study the heterogeneity and interactions of pituitary cells, though its application has been limited thus far. Further utilization of this tool will help uncover the complex physiological and pathological mechanisms of the pituitary, advancing research and treatment of pituitary diseases.

Inflammatory Alterations to Renal Lymphatic Endothelial Cell Gene Expression in Mouse Models of Hypertension

INTRODUCTION

Hypertension (HTN) is a major cardiovascular disease that can cause and be worsened by renal damage and inflammation. We previously reported that renal lymphatic endothelial cells (LECs) increase in response to HTN and that augmenting lymphangiogenesis in the kidneys reduces blood pressure and renal pro-inflammatory immune cells in mice with various forms of HTN. Our aim was to evaluate the specific changes that renal LECs undergo in HTN.

METHODS

We performed single-cell RNA sequencing. Using the angiotensin II-induced and salt-sensitive mouse models of HTN, we isolated renal CD31+ and podoplanin+ cells.

RESULTS

Sequencing of these cells revealed three distinct cell types with unique expression profiles, including LECs. The number and transcriptional diversity of LECs increased in samples from mice with HTN, as demonstrated by 597 differentially expressed genes (p < 0.01), 274 significantly enriched pathways (p < 0.01), and 331 regulons with specific enrichment in HTN LECs. These changes demonstrate a profound inflammatory response in renal LECs in HTN, leading to an increase in genes and pathways associated with inflammation-driven growth and immune checkpoint activity in LECs.

CONCLUSION

These results reinforce and help to further explain the benefits of renal LECs and lymphangiogenesis in HTN.

SingleCellGGM enables gene expression program identification from single-cell transcriptomes and facilitates universal cell label transfer

Gene co-expression analysis of single-cell transcriptomes, aiming to define functional relationships between genes, is challenging due to excessive dropout values. Here, we developed a single-cell graphical Gaussian model (SingleCellGGM) algorithm to conduct single-cell gene co-expression network analysis. When applied to mouse single-cell datasets, SingleCellGGM constructed networks from which gene co-expression modules with highly significant functional enrichment were identified. We considered the modules as gene expression programs (GEPs). These GEPs enable direct cell-type annotation of individual cells without cell clustering, and they are enriched with genes required for the functions of the corresponding cells, sometimes at levels greater than 10-fold. The GEPs are conserved across datasets and enable universal cell-type label transfer across different studies. We also proposed a dimension-reduction method through averaging by GEPs for single-cell analysis, enhancing the interpretability of results. Thus, SingleCellGGM offers a unique GEP-based perspective to analyze single-cell transcriptomes and reveals biological insights shared by different single-cell datasets.

A systematic overview of single-cell transcriptomics databases, their use cases, and limitations

Rapid advancements in high-throughput single-cell RNA-seq (scRNA-seq) technologies and experimental protocols have led to the generation of vast amounts of transcriptomic data that populates several online databases and repositories. Here, we systematically examined large-scale scRNA-seq databases, categorizing them based on their scope and purpose such as general, tissue-specific databases, disease-specific databases, cancer-focused databases, and cell type-focused databases. Next, we discuss the technical and methodological challenges associated with curating large-scale scRNA-seq databases, along with current computational solutions. We argue that understanding scRNA-seq databases, including their limitations and assumptions, is crucial for effectively utilizing this data to make robust discoveries and identify novel biological insights. Such platforms can help bridge the gap between computational and wet lab scientists through user-friendly web-based interfaces needed for democratizing access to single-cell data. These platforms would facilitate interdisciplinary research, enabling researchers from various disciplines to collaborate effectively. This review underscores the importance of leveraging computational approaches to unravel the complexities of single-cell data and offers a promising direction for future research in the field.

Systematic evaluation with practical guidelines for single-cell and spatially resolved transcriptomics data simulation under multiple scenarios

BACKGROUND

Single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) have led to groundbreaking advancements in life sciences. To develop bioinformatics tools for scRNA-seq and SRT data and perform unbiased benchmarks, data simulation has been widely adopted by providing explicit ground truth and generating customized datasets. However, the performance of simulation methods under multiple scenarios has not been comprehensively assessed, making it challenging to choose suitable methods without practical guidelines.

RESULTS

We systematically evaluated 49 simulation methods developed for scRNA-seq and/or SRT data in terms of accuracy, functionality, scalability, and usability using 152 reference datasets derived from 24 platforms. SRTsim, scDesign3, ZINB-WaVE, and scDesign2 have the best accuracy performance across various platforms. Unexpectedly, some methods tailored to scRNA-seq data have potential compatibility for simulating SRT data. Lun, SPARSim, and scDesign3-tree outperform other methods under corresponding simulation scenarios. Phenopath, Lun, Simple, and MFA yield high scalability scores but they cannot generate realistic simulated data. Users should consider the trade-offs between method accuracy and scalability (or functionality) when making decisions. Additionally, execution errors are mainly caused by failed parameter estimations and appearance of missing or infinite values in calculations. We provide practical guidelines for method selection, a standard pipeline Simpipe ( https://github.com/duohongrui/simpipe ; https://doi.org/10.5281/zenodo.11178409 ), and an online tool Simsite ( https://www.ciblab.net/software/simshiny/ ) for data simulation.

CONCLUSIONS

No method performs best on all criteria, thus a good-yet-not-the-best method is recommended if it solves problems effectively and reasonably. Our comprehensive work provides crucial insights for developers on modeling gene expression data and fosters the simulation process for users.

Senescent-like microglia limit remyelination through the senescence associated secretory phenotype

The capacity to regenerate myelin in the central nervous system (CNS) diminishes with age. This decline is particularly evident in multiple sclerosis (MS), which has been suggested to exhibit features of accelerated biological aging. Whether cellular senescence, a hallmark of aging, contributes to remyelination impairment remains unknown. Here, we show that senescent cells (SCs) accumulate within demyelinated lesions after injury, and their elimination enhances remyelination in young mice but not in aged mice. In young mice, we observed the upregulation of senescence-associated transcripts primarily in microglia after demyelination, followed by their reduction during remyelination. However, in aged mice, senescence-associated factors persisted within lesions, correlating with inefficient remyelination. We found that SC elimination enhanced remyelination in young mice but was ineffective in aged mice. Proteomic analysis of senescence-associated secretory phenotype (SASP) revealed elevated levels of CCL11/Eotaxin-1 in lesions, which was found to inhibit efficient oligodendrocyte maturation. These results suggest therapeutic targeting of SASP components, such as CCL11, may improve remyelination in aging and MS.

Establishing a conceptual framework for holistic cell states and state transitions

Cell states were traditionally defined by how they looked, where they were located, and what functions they performed. In this post-genomic era, the field is largely focused on a molecular view of cell state. Moving forward, we anticipate that the observables used to define cell states will evolve again as single-cell imaging and analytics are advancing at a breakneck pace via the collection of large-scale, systematic cell image datasets and the application of quantitative image-based data science methods. This is, therefore, a key moment in the arc of cell biological research to develop approaches that integrate the spatiotemporal observables of the physical structure and organization of the cell with molecular observables toward the concept of a holistic cell state. In this perspective, we propose a conceptual framework for holistic cell states and state transitions that is data-driven, practical, and useful to enable integrative analyses and modeling across many data types.

Data normalization for addressing the challenges in the analysis of single-cell transcriptomic datasets

BACKGROUND

Normalization is a critical step in the analysis of single-cell RNA-sequencing (scRNA-seq) datasets. Its main goal is to make gene counts comparable within and between cells. To do so, normalization methods must account for technical and biological variability. Numerous normalization methods have been developed addressing different sources of dispersion and making specific assumptions about the count data.

MAIN BODY

The selection of a normalization method has a direct impact on downstream analysis, for example differential gene expression and cluster identification. Thus, the objective of this review is to guide the reader in making an informed decision on the most appropriate normalization method to use. To this aim, we first give an overview of the different single cell sequencing platforms and methods commonly used including isolation and library preparation protocols. Next, we discuss the inherent sources of variability of scRNA-seq datasets. We describe the categories of normalization methods and include examples of each. We also delineate imputation and batch-effect correction methods. Furthermore, we describe data-driven metrics commonly used to evaluate the performance of normalization methods. We also discuss common scRNA-seq methods and toolkits used for integrated data analysis.

CONCLUSIONS

According to the correction performed, normalization methods can be broadly classified as within and between-sample algorithms. Moreover, with respect to the mathematical model used, normalization methods can further be classified into: global scaling methods, generalized linear models, mixed methods, and machine learning-based methods. Each of these methods depict pros and cons and make different statistical assumptions. However, there is no better performing normalization method. Instead, metrics such as silhouette width, K-nearest neighbor batch-effect test, or Highly Variable Genes are recommended to assess the performance of normalization methods.

Towards early diagnosis of Alzheimer's disease: advances in immune-related blood biomarkers and computational approaches

Alzheimer's disease has an increasing prevalence in the population world-wide, yet current diagnostic methods based on recommended biomarkers are only available in specialized clinics. Due to these circumstances, Alzheimer's disease is usually diagnosed late, which contrasts with the currently available treatment options that are only effective for patients at an early stage. Blood-based biomarkers could fill in the gap of easily accessible and low-cost methods for early diagnosis of the disease. In particular, immune-based blood-biomarkers might be a promising option, given the recently discovered cross-talk of immune cells of the central nervous system with those in the peripheral immune system. Here, we give a background on recent advances in research on brain-immune system cross-talk in Alzheimer's disease and review machine learning approaches, which can combine multiple biomarkers with further information (e.g. age, sex, APOE genotype) into predictive models supporting an earlier diagnosis. In addition, mechanistic modeling approaches, such as agent-based modeling open the possibility to model and analyze cell dynamics over time. This review aims to provide an overview of the current state of immune-system related blood-based biomarkers and their potential for the early diagnosis of Alzheimer's disease.

Library size confounds biology in spatial transcriptomics data

Spatial molecular data has transformed the study of disease microenvironments, though, larger datasets pose an analytics challenge prompting the direct adoption of single-cell RNA-sequencing tools including normalization methods. Here, we demonstrate that library size is associated with tissue structure and that normalizing these effects out using commonly applied scRNA-seq normalization methods will negatively affect spatial domain identification. Spatial data should not be specifically corrected for library size prior to analysis, and algorithms designed for scRNA-seq data should be adopted with caution.

A Systematic Overview of Single-Cell Transcriptomics Databases, their Use cases, and Limitations

Rapid advancements in high-throughput single-cell RNA-seq (scRNA-seq) technologies and experimental protocols have led to the generation of vast amounts of genomic data that populates several online databases and repositories. Here, we systematically examined large-scale scRNA-seq databases, categorizing them based on their scope and purpose such as general, tissue-specific databases, disease-specific databases, cancer-focused databases, and cell type-focused databases. Next, we discuss the technical and methodological challenges associated with curating large-scale scRNA-seq databases, along with current computational solutions. We argue that understanding scRNA-seq databases, including their limitations and assumptions, is crucial for effectively utilizing this data to make robust discoveries and identify novel biological insights. Furthermore, we propose that bridging the gap between computational and wet lab scientists through user-friendly web-based platforms is needed for democratizing access to single-cell data. These platforms would facilitate interdisciplinary research, enabling researchers from various disciplines to collaborate effectively. This review underscores the importance of leveraging computational approaches to unravel the complexities of single-cell data and offers a promising direction for future research in the field.

SuperCellCyto: enabling efficient analysis of large scale cytometry datasets

Advancements in cytometry technologies have enabled quantification of up to 50 proteins across millions of cells at single cell resolution. Analysis of cytometry data routinely involves tasks such as data integration, clustering, and dimensionality reduction. While numerous tools exist, many require extensive run times when processing large cytometry data containing millions of cells. Existing solutions, such as random subsampling, are inadequate as they risk excluding rare cell subsets. To address this, we propose SuperCellCyto, an R package that builds on the SuperCell tool which groups highly similar cells into supercells. SuperCellCyto is available on GitHub ( https://github.com/phipsonlab/SuperCellCyto ) and Zenodo ( https://doi.org/10.5281/zenodo.10521294 ).

Defining and benchmarking open problems in single-cell analysis

With the growing number of single-cell analysis tools, benchmarks are increasingly important to guide analysis and method development. However, a lack of standardisation and extensibility in current benchmarks limits their usability, longevity, and relevance to the community. We present Open Problems, a living, extensible, community-guided benchmarking platform including 10 current single-cell tasks that we envision will raise standards for the selection, evaluation, and development of methods in single-cell analysis.

Single-cell and spatial transcriptomics: Bridging current technologies with long-read sequencing

Single-cell technologies have transformed biomedical research over the last decade, opening up new possibilities for understanding cellular heterogeneity, both at the genomic and transcriptomic level. In addition, more recent developments of spatial transcriptomics technologies have made it possible to profile cells in their tissue context. In parallel, there have been substantial advances in sequencing technologies, and the third generation of methods are able to produce reads that are tens of kilobases long, with error rates matching the second generation short reads. Long reads technologies make it possible to better map large genome rearrangements and quantify isoform specific abundances. This further improves our ability to characterize functionally relevant heterogeneity. Here, we show how researchers have begun to combine single-cell, spatial transcriptomics, and long-read technologies, and how this is resulting in powerful new approaches to profiling both the genome and the transcriptome. We discuss the achievements so far, and we highlight remaining challenges and opportunities.

A machine-readable specification for genomics assays

MOTIVATION

Understanding the structure of sequenced fragments from genomics libraries is essential for accurate read preprocessing. Currently, different assays and sequencing technologies require custom scripts and programs that do not leverage the common structure of sequence elements present in genomics libraries.

RESULTS

We present seqspec, a machine-readable specification for libraries produced by genomics assays that facilitates standardization of preprocessing and enables tracking and comparison of genomics assays.

AVAILABILITY AND IMPLEMENTATION

The specification and associated seqspec command line tool is available at https://www.doi.org/10.5281/zenodo.10213865.

Recent advances in differential expression analysis for single-cell RNA-seq and spatially resolved transcriptomic studies

Differential expression (DE) analysis is a necessary step in the analysis of single-cell RNA sequencing (scRNA-seq) and spatially resolved transcriptomics (SRT) data. Unlike traditional bulk RNA-seq, DE analysis for scRNA-seq or SRT data has unique characteristics that may contribute to the difficulty of detecting DE genes. However, the plethora of DE tools that work with various assumptions makes it difficult to choose an appropriate one. Furthermore, a comprehensive review on detecting DE genes for scRNA-seq data or SRT data from multi-condition, multi-sample experimental designs is lacking. To bridge such a gap, here, we first focus on the challenges of DE detection, then highlight potential opportunities that facilitate further progress in scRNA-seq or SRT analysis, and finally provide insights and guidance in selecting appropriate DE tools or developing new computational DE methods.

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.

The impacts of active and self-supervised learning on efficient annotation of single-cell expression data

A crucial step in the analysis of single-cell data is annotating cells to cell types and states. While a myriad of approaches has been proposed, manual labeling of cells to create training datasets remains tedious and time-consuming. In the field of machine learning, active and self-supervised learning methods have been proposed to improve the performance of a classifier while reducing both annotation time and label budget. However, the benefits of such strategies for single-cell annotation have yet to be evaluated in realistic settings. Here, we perform a comprehensive benchmarking of active and self-supervised labeling strategies across a range of single-cell technologies and cell type annotation algorithms. We quantify the benefits of active learning and self-supervised strategies in the presence of cell type imbalance and variable similarity. We introduce adaptive reweighting, a heuristic procedure tailored to single-cell data-including a marker-aware version-that shows competitive performance with existing approaches. In addition, we demonstrate that having prior knowledge of cell type markers improves annotation accuracy. Finally, we summarize our findings into a set of recommendations for those implementing cell type annotation procedures or platforms. An R package implementing the heuristic approaches introduced in this work may be found at https://github.com/camlab-bioml/leader .

Deciphering the immune modulation through deep transcriptomic profiling and therapeutic implications of DNA damage repair pattern in hepatocellular carcinoma

AIMS

DNA damage repair (DDR) plays a pivotal role in hepatocellular carcinoma (HCC), driving oncogenesis, progression, and therapeutic response. However, the mechanisms of DDR mediated immune cells and immuno-modulatory pathways in HCC are yet ill-defined.

METHODS

Our study introduces an innovative deep machine learning framework for precise DDR assessment, utilizing single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data. Single-cell RNA sequencing data were obtained and in total 85,628 cells of primary or post-immunotherapy cases were analyzed. Large-scale HCC datasets, including 1027 patients in house together with public datasets, were used for 101 machine-learning models and a novel DDR feature was derived at single-cell resolution (DDRscore). Druggable targets were predicted using the reverse phase protein array (RPPA) proteomic profiling of 169 HCC patients and RNA-seq data from 22 liver cancer cell lines.

RESULTS

Our investigation reveals a dynamic interplay of DDR with natural killer cells and B cells in the primary HCC microenvironment, shaping a tumor-promoting immune milieu through metabolic programming. Analysis of HCC post-immunotherapy demonstrates elevated DDR levels that induces epithelial-mesenchymal transition and fibroblast-like transformation, reshaping the fibrotic tumor microenvironment. Conversely, attenuated DDR promotes antigen cross-presentation by dendritic cells and CD8+ T cells, modulating the inflammatory tumor microenvironment. Regulatory network analysis identifies the CXCL10-CXCR3 axis as a key determinant of immunotherapeutic response in low DDR HCC, potentially regulated by transcription factors GATA3, REL, and TBX21. Using machine learning techniques by combining bulk RNA-seq data in house together with public datasets, we introduce DDRscore, a robust consensus DDR scoring system to predict overall survival and resistance to PD-1 therapy in HCC patients. Finally, we identify BRAF as a potential therapeutic target for high DDRscore patients.

CONCLUSION

Our comprehensive findings advance our understanding of DDR and the tumor microenvironment in HCC, providing insights into immune regulatory mechanisms mediated via DDR pathways.

A review on gene regulatory network reconstruction algorithms based on single cell RNA sequencing

BACKGROUND

Understanding gene regulatory networks (GRNs) is essential for unraveling the molecular mechanisms governing cellular behavior. With the advent of high-throughput transcriptome measurement technology, researchers have aimed to reverse engineer the biological systems, extracting gene regulatory rules from their outputs, which represented by gene expression data. Bulk RNA sequencing, a widely used method for measuring gene expression, has been employed for GRN reconstruction. However, it falls short in capturing dynamic changes in gene expression at the level of individual cells since it averages gene expression across mixed cell populations.

OBJECTIVE

In this review, we provide an overview of 15 GRN reconstruction tools and discuss their respective strengths and limitations, particularly in the context of single cell RNA sequencing (scRNA-seq).

METHODS

Recent advancements in scRNA-seq break new ground of GRN reconstruction. They offer snapshots of the individual cell transcriptomes and capturing dynamic changes. We emphasize how these technological breakthroughs have enhanced GRN reconstruction.

CONCLUSION

GRN reconstructors can be classified based on their requirement for cellular trajectory, which represents a dynamical cellular process including differentiation, aging, or disease progression. Benchmarking studies support the superiority of GRN reconstructors that do not require trajectory analysis in identifying regulator-target relationships. However, methods equipped with trajectory analysis demonstrate better performance in identifying key regulatory factors. In conclusion, researchers should select a suitable GRN reconstructor based on their specific research objectives.

Fishing Innate Immune System Properties through the Transcriptomic Single-Cell Data of Teleostei

The innate immune system is the first line of defense in multicellular organisms. Danio rerio is widely considered a promising model for IIS-related research, with the most amount of scRNAseq data available among Teleostei. We summarized the scRNAseq and spatial transcriptomics experiments related to the IIS for zebrafish and other Teleostei from the GEO NCBI and the Single-Cell Expression Atlas. We found a considerable number of scRNAseq experiments at different stages of zebrafish development in organs such as the kidney, liver, stomach, heart, and brain. These datasets could be further used to conduct large-scale meta-analyses and to compare the IIS of zebrafish with the mammalian one. However, only a small number of scRNAseq datasets are available for other fish (turbot, salmon, cavefish, and dark sleeper). Since fish biology is very diverse, it would be a major mistake to use zebrafish alone in fish immunology studies. In particular, there is a special need for new scRNAseq experiments involving nonmodel Teleostei, e.g., long-lived species, cancer-resistant fish, and various fish ecotypes.

SingleScan: a comprehensive resource for single-cell sequencing data processing and mining

Single-cell sequencing has shed light on previously inaccessible biological questions from different fields of research, including organism development, immune function, and disease progression. The number of single-cell-based studies increased dramatically over the past decade. Several new methods and tools have been continuously developed, making it extremely tricky to navigate this research landscape and develop an up-to-date workflow to analyze single-cell sequencing data, particularly for researchers seeking to enter this field without computational experience. Moreover, choosing appropriate tools and optimal parameters to meet the demands of researchers represents a major challenge in processing single-cell sequencing data. However, a specific resource for easy access to detailed information on single-cell sequencing methods and data processing pipelines is still lacking. In the present study, an online resource called SingleScan was developed to curate all up-to-date single-cell transcriptome/genome analyzing tools and pipelines. All the available tools were categorized according to their main tasks, and several typical workflows for single-cell data analysis were summarized. In addition, spatial transcriptomics, which is a breakthrough molecular analysis method that enables researchers to measure all gene activity in tissue samples and map the site of activity, was included along with a portion of single-cell and spatial analysis solutions. For each processing step, the available tools and specific parameters used in published articles are provided and how these parameters affect the results is shown in the resource. All information used in the resource was manually extracted from related literature. An interactive website was designed for data retrieval, visualization, and download. By analyzing the included tools and literature, users can gain insights into the trends of single-cell studies and easily grasp the specific usage of a specific tool. SingleScan will facilitate the analysis of single-cell sequencing data and promote the development of new tools to meet the growing and diverse needs of the research community. The SingleScan database is publicly accessible via the website at http://cailab.labshare.cn/SingleScan .

Single-Cell RNA-Seq Analysis Identifies Angiotensinogen and Galanin as Unique Molecular Markers of Acinar Cells in Murine Salivary Glands

The submandibular gland (SMG) and sublingual gland (SLG) are two of three major salivary glands in mammals and comprise serous and mucous acinar cells. The two glands share some functional properties, which are largely dependent on the types of acinar cells. In recent years, while ScRNA-seq (single-cell sequencing) with a 10 × platform has been used to explore molecular markers in salivary glands, few studies have examined the acinar heterogeneity and unique molecular markers between SMG and SLG. This study aimed to identify the molecular markers of acinar cells in the SLG and SMG. We performed ScRNA-seq analyses in 4-week-old mice and verified the screened molecular markers using reverse transcription-quantitative real-time PCR, immunohistochemistry, and immunofluorescence. Our results showed prominently heterogeneous acinar cells, although there was great similarity in the cluster composition between the two glands at 4 weeks. Furthermore, we demonstrated that Agt is a specific marker of SMG serous acinar cells, whereas Gal is a specific marker of SLG mucous acinar cells. Trajectory inference revealed that Agt and Gal represent two types of differential acinar cell clusters during late development in adults. Thus, we reveal previously unknown specific markers for salivary acinar cell diversity, which has extensive implications for their further functional research.

Novel insights into molecular signatures and pathogenic cell populations shared by systemic lupus erythematosus and vascular dementia

Although vascular dementia (VD) and systemic lupus erythematosus (SLE) may share immune-mediated pathophysiologic processes, the underlying mechanisms are unclear. This study investigated shared gene signatures in SLE versus VD, as well as their potential molecular mechanisms. Bulk RNA sequencing (RNAseq) and single-cell or single-nucleus RNAseq (sc/snRNAseq) datasets from SLE blood samples and VD brain samples were obtained from Gene Expression Omnibus. The identification of genes associated with both SLE and VD was performed using the weighted gene co-expression network analysis (WGCNA) and machine learning algorithms. For the sc/snRNAseq data, an unbiased clustering pipeline based on Seurat and CellChat was used to determine the cellular landscape profile and examine intracellular communication, respectively. The results were subsequently validated using a mice model of SLE with cognitive dysfunction (female MRL/lpr mice). WGCNA and machine learning identified C1QA, LY96, CD163, and MS4A4A as key genes for SLE and VD. sc/snRNAseq analyses revealed that CD163 and MS4A4A were upregulated in mononuclear phagocytes (MPs) from SLE and VD samples and were associated with monocyte-macrophage differentiation. Intriguingly, LGALS9-associated molecular pathway, as the only signaling pathway common between SLE and VD via CellChat analysis, exhibited significant upregulation in cortical microglia of MRL/lpr mice. Our analyses identified C1QA, LY96, CD163, and MS4A4A as potential biomarkers for SLE and VD. Moreover, the upregulation of CD163/MS4A4A and activation of LGALS9 signaling in MPs may contribute to the pathogenesis of VD with SLE. These findings offer novel insight into the mechanisms underlying VD in SLE patients.

Single-cell transcriptomics of the ocular anterior segment: a comprehensive review

Elucidating the cellular and genetic composition of ocular tissues is essential for uncovering the pathophysiology of ocular diseases. Since the introduction of single-cell RNA sequencing (scRNA-seq) in 2009, vision researchers have performed extensive single-cell analyses to better understand transcriptome complexity and heterogeneity of ocular structures. This technology has revolutionized our ability to identify rare cell populations and to make cross-species comparisons of gene expression in both steady state and disease conditions. Importantly, single-cell transcriptomic analyses have enabled the identification of cell-type specific gene markers and signalling pathways between ocular cell populations. While most scRNA-seq studies have been conducted on retinal tissues, large-scale transcriptomic atlases pertaining to the ocular anterior segment have also been constructed in the past three years. This timely review provides vision researchers with an overview of scRNA-seq experimental design, technical limitations, and clinical applications in a variety of anterior segment-related ocular pathologies. We review open-access anterior segment-related scRNA-seq datasets and illustrate how scRNA-seq can be an indispensable tool for the development of targeted therapeutics.

SCALA: A complete solution for multimodal analysis of single-cell Next Generation Sequencing data

Analysis and interpretation of high-throughput transcriptional and chromatin accessibility data at single-cell (sc) resolution are still open challenges in the biomedical field. The existence of countless bioinformatics tools, for the different analytical steps, increases the complexity of data interpretation and the difficulty to derive biological insights. In this article, we present SCALA, a bioinformatics tool for analysis and visualization of single-cell RNA sequencing (scRNA-seq) and Assay for Transposase-Accessible Chromatin using sequencing (scATAC-seq) datasets, enabling either independent or integrative analysis of the two modalities. SCALA combines standard types of analysis by integrating multiple software packages varying from quality control to the identification of distinct cell populations and cell states. Additional analysis options enable functional enrichment, cellular trajectory inference, ligand-receptor analysis, and regulatory network reconstruction. SCALA is fully parameterizable, presenting data in tabular format and producing publication-ready visualizations. The different available analysis modules can aid biomedical researchers in exploring, analyzing, and visualizing their data without any prior experience in coding. We demonstrate the functionality of SCALA through two use-cases related to TNF-driven arthritic mice, handling both scRNA-seq and scATAC-seq datasets. SCALA is developed in R, Shiny and JavaScript and is mainly available as a standalone version, while an online service of more limited capacity can be found at http://scala.pavlopouloslab.info or https://scala.fleming.gr.

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.

Artificial Intelligence Models for Cell Type and Subtype Identification Based on Single-Cell RNA Sequencing Data in Vision Science

Single-cell RNA sequencing (scRNA-seq) provides a high throughput, quantitative and unbiased framework for scientists in many research fields to identify and characterize cell types within heterogeneous cell populations from various tissues. However, scRNA-seq based identification of discrete cell-types is still labor intensive and depends on prior molecular knowledge. Artificial intelligence has provided faster, more accurate, and user-friendly approaches for cell-type identification. In this review, we discuss recent advances in cell-type identification methods using artificial intelligence techniques based on single-cell and single-nucleus RNA sequencing data in vision science. The main purpose of this review paper is to assist vision scientists not only to select suitable datasets for their problems, but also to be aware of the appropriate computational tools to perform their analysis. Developing novel methods for scRNA-seq data analysis remains to be addressed in future studies.

A Roadmap for the Human Gut Cell Atlas

The number of studies investigating the human gastrointestinal tract using various single-cell profiling methods has increased substantially in the past few years. Although this increase provides a unique opportunity for the generation of the first comprehensive Human Gut Cell Atlas (HGCA), there remains a range of major challenges ahead. Above all, the ultimate success will largely depend on a structured and coordinated approach that aligns global efforts undertaken by a large number of research groups. In this Roadmap, we discuss a comprehensive forward-thinking direction for the generation of the HGCA on behalf of the Gut Biological Network of the Human Cell Atlas. Based on the consensus opinion of experts from across the globe, we outline the main requirements for the first complete HGCA by summarizing existing data sets and highlighting anatomical regions and/or tissues with limited coverage. We provide recommendations for future studies and discuss key methodologies and the importance of integrating the healthy gut atlas with related diseases and gut organoids. Importantly, we critically overview the computational tools available and provide recommendations to overcome key challenges.

Single-Cell RNA Sequencing for Studying Human Cancers

Since the first publication a decade ago describing the use of single-cell RNA sequencing (scRNA-seq) in the context of cancer, over 200 datasets and thousands of scRNA-seq studies have been published in cancer biology. scRNA-seq technologies have been applied across dozens of cancer types and a diverse array of study designs to improve our understanding of tumor biology, the tumor microenvironment, and therapeutic responses, and scRNA-seq is on the verge of being used to improve decision-making in the clinic. Computational methodologies and analytical pipelines are key in facilitating scRNA-seq research. Numerous computational methods utilizing the most advanced tools in data science have been developed to extract meaningful insights. Here, we review the advancements in cancer biology gained by scRNA-seq and discuss the computational challenges of the technology that are specific to cancer research.

Exploring the optimization of autoencoder design for imputing single-cell RNA sequencing data

Autoencoders are the backbones of many imputation methods that aim to relieve the sparsity issue in single-cell RNA sequencing (scRNA-seq) data. The imputation performance of an autoencoder relies on both the neural network architecture and the hyperparameter choice. So far, literature in the single-cell field lacks a formal discussion on how to design the neural network and choose the hyperparameters. Here, we conducted an empirical study to answer this question. Our study used many real and simulated scRNA-seq datasets to examine the impacts of the neural network architecture, the activation function, and the regularization strategy on imputation accuracy and downstream analyses. Our results show that (i) deeper and narrower autoencoders generally lead to better imputation performance; (ii) the sigmoid and tanh activation functions consistently outperform other commonly used functions including ReLU; (iii) regularization improves the accuracy of imputation and downstream cell clustering and DE gene analyses. Notably, our results differ from common practices in the computer vision field regarding the activation function and the regularization strategy. Overall, our study offers practical guidance on how to optimize the autoencoder design for scRNA-seq data imputation.

A machine-readable specification for genomics assays

Understanding the structure of sequenced fragments from genomics libraries is essential for accurate read preprocessing. Currently, different assays and sequencing technologies require custom scripts and programs that do not leverage the common structure of sequence elements present in genomics libraries. We present seqspec, a machine-readable specification for libraries produced by genomics assays that facilitates standardization of preprocessing and enables tracking and comparison of genomics assays. The specification and associated seqspec command line tool is available at https://github.com/IGVF/seqspec.

Respiratory epithelial cell types, states and fates in the era of single-cell RNA-sequencing

Standalone and consortia-led single-cell atlases of healthy and diseased human airways generated with single-cell RNA-sequencing (scRNA-seq) have ushered in a new era in respiratory research. Numerous discoveries, including the pulmonary ionocyte, potentially novel cell fates, and a diversity of cell states among common and rare epithelial cell types have highlighted the extent of cellular heterogeneity and plasticity in the respiratory tract. scRNA-seq has also played a pivotal role in our understanding of host-virus interactions in coronavirus disease 2019 (COVID-19). However, as our ability to generate large quantities of scRNA-seq data increases, along with a growing number of scRNA-seq protocols and data analysis methods, new challenges related to the contextualisation and downstream applications of insights are arising. Here, we review the fundamental concept of cellular identity from the perspective of single-cell transcriptomics in the respiratory context, drawing attention to the need to generate reference annotations and to standardise the terminology used in literature. Findings about airway epithelial cell types, states and fates obtained from scRNA-seq experiments are compared and contrasted with information accumulated through the use of conventional methods. This review attempts to discuss major opportunities and to outline some of the key limitations of the modern-day scRNA-seq that need to be addressed to enable efficient and meaningful integration of scRNA-seq data from different platforms and studies, with each other as well as with data from other high-throughput sequencing-based genomic, transcriptomic and epigenetic analyses.

Biomarkers as Biomedical Bioindicators: Approaches and Techniques for the Detection, Analysis, and Validation of Novel Biomarkers of Diseases

A biomarker is any measurable biological moiety that can be assessed and measured as a potential index of either normal or abnormal pathophysiology or pharmacological responses to some treatment regimen. Every tissue in the body has a distinct biomolecular make-up, which is known as its biomarkers, which possess particular features, viz., the levels or activities (the ability of a gene or protein to carry out a particular body function) of a gene, protein, or other biomolecules. A biomarker refers to some feature that can be objectively quantified by various biochemical samples and evaluates the exposure of an organism to normal or pathological procedures or their response to some drug interventions. An in-depth and comprehensive realization of the significance of these biomarkers becomes quite important for the efficient diagnosis of diseases and for providing the appropriate directions in case of multiple drug choices being presently available, which can benefit any patient. Presently, advancements in omics technologies have opened up new possibilities to obtain novel biomarkers of different types, employing genomic strategies, epigenetics, metabolomics, transcriptomics, lipid-based analysis, protein studies, etc. Particular biomarkers for specific diseases, their prognostic capabilities, and responses to therapeutic paradigms have been applied for screening of various normal healthy, as well as diseased, tissue or serum samples, and act as appreciable tools in pharmacology and therapeutics, etc. In this review, we have summarized various biomarker types, their classification, and monitoring and detection methods and strategies. Various analytical techniques and approaches of biomarkers have also been described along with various clinically applicable biomarker sensing techniques which have been developed in the recent past. A section has also been dedicated to the latest trends in the formulation and designing of nanotechnology-based biomarker sensing and detection developments in this field.

Meta-analysis of (single-cell method) benchmarks reveals the need for extensibility and interoperability

Computational methods represent the lifeblood of modern molecular biology. Benchmarking is important for all methods, but with a focus here on computational methods, benchmarking is critical to dissect important steps of analysis pipelines, formally assess performance across common situations as well as edge cases, and ultimately guide users on what tools to use. Benchmarking can also be important for community building and advancing methods in a principled way. We conducted a meta-analysis of recent single-cell benchmarks to summarize the scope, extensibility, and neutrality, as well as technical features and whether best practices in open data and reproducible research were followed. The results highlight that while benchmarks often make code available and are in principle reproducible, they remain difficult to extend, for example, as new methods and new ways to assess methods emerge. In addition, embracing containerization and workflow systems would enhance reusability of intermediate benchmarking results, thus also driving wider adoption.

Multiomics Empowers Predictive Pancreatic Cancer Immunotherapy

Advances in cancer immunotherapy, particularly immune checkpoint inhibitors, have dramatically improved the prognosis for patients with metastatic melanoma and other previously incurable cancers. However, patients with pancreatic ductal adenocarcinoma (PDAC) generally do not respond to these therapies. PDAC is exceptionally difficult to treat because of its often late stage at diagnosis, modest mutation burden, and notoriously complex and immunosuppressive tumor microenvironment. Simultaneously interrogating features of cancer, immune, and other cellular components of the PDAC tumor microenvironment is therefore crucial for identifying biomarkers of immunotherapeutic resistance and response. Notably, single-cell and multiomics technologies, along with the analytical tools for interpreting corresponding data, are facilitating discoveries of the systems-level cellular and molecular interactions contributing to the overall resistance of PDAC to immunotherapy. Thus, in this review, we will explore how multiomics and single-cell analyses provide the unprecedented opportunity to identify biomarkers of resistance and response to successfully sensitize PDAC to immunotherapy.

The shaky foundations of simulating single-cell RNA sequencing data

BACKGROUND

With the emergence of hundreds of single-cell RNA-sequencing (scRNA-seq) datasets, the number of computational tools to analyze aspects of the generated data has grown rapidly. As a result, there is a recurring need to demonstrate whether newly developed methods are truly performant-on their own as well as in comparison to existing tools. Benchmark studies aim to consolidate the space of available methods for a given task and often use simulated data that provide a ground truth for evaluations, thus demanding a high quality standard results credible and transferable to real data.

RESULTS

Here, we evaluated methods for synthetic scRNA-seq data generation in their ability to mimic experimental data. Besides comparing gene- and cell-level quality control summaries in both one- and two-dimensional settings, we further quantified these at the batch- and cluster-level. Secondly, we investigate the effect of simulators on clustering and batch correction method comparisons, and, thirdly, which and to what extent quality control summaries can capture reference-simulation similarity.

CONCLUSIONS

Our results suggest that most simulators are unable to accommodate complex designs without introducing artificial effects, they yield over-optimistic performance of integration and potentially unreliable ranking of clustering methods, and it is generally unknown which summaries are important to ensure effective simulation-based method comparisons.

Single-cell RNA sequencing to identify cellular heterogeneity and targets in cardiovascular diseases: from bench to bedside

The mechanisms of cardiovascular diseases (CVDs) remain incompletely elucidated. Single-cell RNA sequencing (scRNA-seq) has enabled the profiling of single-cell transcriptomes at unprecedented resolution and throughput, which is critical for deciphering cardiovascular cellular heterogeneity and underlying disease mechanisms, thereby facilitating the development of therapeutic strategies. In this review, we summarize cellular heterogeneity in cardiovascular homeostasis and diseases as well as the discovery of potential disease targets based on scRNA-seq, and yield new insights into the promise of scRNA-seq technology in precision medicine and clinical application.

Computational Methods for Single-cell DNA Methylome Analysis

Dissecting intercellular epigenetic differences is key to understanding tissue heterogeneity. Recent advances in single-cell DNA methylome profiling have presented opportunities to resolve this heterogeneity at the maximum resolution. While these advances enable us to explore frontiers of chromatin biology and better understand cell lineage relationships, they pose new challenges in data processing and interpretation. This review surveys the current state of computational tools developed for single-cell DNA methylome data analysis. We discuss critical components of single-cell DNA methylome data analysis, including data preprocessing, quality control, imputation, dimensionality reduction, cell clustering, supervised cell annotation, cell lineage reconstruction, gene activity scoring, and integration with transcriptome data. We also highlight unique aspects of single-cell DNA methylome data analysis and discuss how techniques common to other single-cell omics data analyses can be adapted to analyze DNA methylomes. Finally, we discuss existing challenges and opportunities for future development.

A novel Bayesian framework for harmonizing information across tissues and studies to increase cell type deconvolution accuracy

Computational cell type deconvolution on bulk transcriptomics data can reveal cell type proportion heterogeneity across samples. One critical factor for accurate deconvolution is the reference signature matrix for different cell types. Compared with inferring reference signature matrices from cell lines, rapidly accumulating single-cell RNA-sequencing (scRNA-seq) data provide a richer and less biased resource. However, deriving cell type signature from scRNA-seq data is challenging due to high biological and technical noises. In this article, we introduce a novel Bayesian framework, tranSig, to improve signature matrix inference from scRNA-seq by leveraging shared cell type-specific expression patterns across different tissues and studies. Our simulations show that tranSig is robust to the number of signature genes and tissues specified in the model. Applications of tranSig to bulk RNA sequencing data from peripheral blood, bronchoalveolar lavage and aorta demonstrate its accuracy and power to characterize biological heterogeneity across groups. In summary, tranSig offers an accurate and robust approach to defining gene expression signatures of different cell types, facilitating improved in silico cell type deconvolutions.

Non-negative low-rank representation based on dictionary learning for single-cell RNA-sequencing data analysis

In the analysis of single-cell RNA-sequencing (scRNA-seq) data, how to effectively and accurately identify cell clusters from a large number of cell mixtures is still a challenge. Low-rank representation (LRR) method has achieved excellent results in subspace clustering. But in previous studies, most LRR-based methods usually choose the original data matrix as the dictionary. In addition, the methods based on LRR usually use spectral clustering algorithm to complete cell clustering. Therefore, there is a matching problem between the spectral clustering method and the affinity matrix, which is difficult to ensure the optimal effect of clustering. Considering the above two points, we propose the DLNLRR method to better identify the cell type. First, DLNLRR can update the dictionary during the optimization process instead of using the predefined fixed dictionary, so it can realize dictionary learning and LRR learning at the same time. Second, DLNLRR can realize subspace clustering without relying on spectral clustering algorithm, that is, we can perform clustering directly based on the low-rank matrix. Finally, we carry out a large number of experiments on real single-cell datasets and experimental results show that DLNLRR is superior to other scRNA-seq data analysis algorithms in cell type identification.