How to Get Started with Single Cell RNA Sequencing Data Analysis

scIALM: A method for sparse scRNA-seq expression matrix imputation using the Inexact Augmented Lagrange Multiplier with low error

Single-cell RNA sequencing (scRNA-seq) is a high-throughput sequencing technology that quantifies gene expression profiles of specific cell populations at the single-cell level, providing a foundation for studying cellular heterogeneity and patient pathological characteristics. It is effective for developmental, fertility, and disease studies. However, the cell-gene expression matrix of single-cell sequencing data is often sparse and contains numerous zero values. Some of the zero values derive from noise, where dropout noise has a large impact on downstream analysis. In this paper, we propose a method named scIALM for imputation recovery of sparse single-cell RNA data expression matrices, which employs the Inexact Augmented Lagrange Multiplier method to use sparse but clean (accurate) data to recover unknown entries in the matrix. We perform experimental analysis on four datasets, calling the expression matrix after Quality Control (QC) as the original matrix, and comparing the performance of scIALM with six other methods using mean squared error (MSE), mean absolute error (MAE), Pearson correlation coefficient (PCC), and cosine similarity (CS). Our results demonstrate that scIALM accurately recovers the original data of the matrix with an error of 10e-4, and the mean value of the four metrics reaches 4.5072 (MSE), 0.765 (MAE), 0.8701 (PCC), 0.8896 (CS). In addition, at 10%-50% random masking noise, scIALM is the least sensitive to the masking ratio. For downstream analysis, this study uses adjusted rand index (ARI) and normalized mutual information (NMI) to evaluate the clustering effect, and the results are improved on three datasets containing real cluster labels.

NTCdb: Single-cell transcriptome database of human inflammatory-associated diseases

With both the advancement of technology and the decline in costs, single-cell transcriptomics sequencing has become widespread in the biomedical area in recent years. It can facilitate the pathogenic characteristics at the single-cell level, which will assist clinical researchers in exploring the mechanism of diseases. As a result, single-cell transcriptome data based on clinical samples grew exponentially. However, there is still a lack of a comprehensive database about immunocytes in inflammatory-associated diseases. To address this deficiency, we propose a human inflammatory-associated disease-based single-cell transcriptome database, NTCdb (www.ntcdb.org.cn). NTCdb integrates the open-source data of 1,023,166 cells derived from 11 tissues of 17 inflammatory-associated diseases in a uniform pipeline. It provides a set of analyzing results, including cell communication analysis, enrichment analysis, and Pseudo-Time analysis, to obtain various characteristics of immune cells in inflammatory-associated disease. Taking COVID-19 as a case study, NTCdb displays important information including potentially significant functions of certain cells, genes, and signaling pathways, as well as the commonalities of specific immunocytes between different inflammatory-associated disease.

Advancing skeletal health and disease research with single-cell RNA sequencing

Orthopedic conditions have emerged as global health concerns, impacting approximately 1.7 billion individuals worldwide. However, the limited understanding of the underlying pathological processes at the cellular and molecular level has hindered the development of comprehensive treatment options for these disorders. The advent of single-cell RNA sequencing (scRNA-seq) technology has revolutionized biomedical research by enabling detailed examination of cellular and molecular diversity. Nevertheless, investigating mechanisms at the single-cell level in highly mineralized skeletal tissue poses technical challenges. In this comprehensive review, we present a streamlined approach to obtaining high-quality single cells from skeletal tissue and provide an overview of existing scRNA-seq technologies employed in skeletal studies along with practical bioinformatic analysis pipelines. By utilizing these methodologies, crucial insights into the developmental dynamics, maintenance of homeostasis, and pathological processes involved in spine, joint, bone, muscle, and tendon disorders have been uncovered. Specifically focusing on the joint diseases of degenerative disc disease, osteoarthritis, and rheumatoid arthritis using scRNA-seq has provided novel insights and a more nuanced comprehension. These findings have paved the way for discovering novel therapeutic targets that offer potential benefits to patients suffering from diverse skeletal disorders.

Single-cell sequencing advances in research on mesenchymal stem/stromal cells

Mesenchymal stem/stromal cells (MSCs), originating from the mesoderm, represent a multifunctional stem cell population capable of differentiating into diverse cell types and exhibiting a wide range of biological functions. Despite more than half a century of research, MSCs continue to be among the most extensively studied cell types in clinical research projects globally. However, their significant heterogeneity and phenotypic instability have significantly hindered their exploration and application. Single-cell sequencing technology emerges as a powerful tool to address these challenges, offering precise dissection of complex cellular samples. It uncovers the genetic structure and gene expression status of individual contained cells on a massive scale and reveals the heterogeneity among these cells. It links the molecular characteristics of MSCs with their clinical applications, contributing to the advancement of regenerative medicine. With the development and cost reduction of single-cell analysis techniques, sequencing technology is now widely applied in fundamental research and clinical trials. This study aimed to review the application of single-cell sequencing in MSC research and assess its prospects.

Advancements in Single-Cell RNA Sequencing Research for Neurological Diseases

Neurological diseases are a major cause of the global burden of disease. Although the mechanisms of the occurrence and development of neurological diseases are not fully clear, most of them are associated with cells mediating neuroinflammation. Yet medications and other therapeutic options to improve treatment are still very limited. Single-cell RNA sequencing (scRNA-seq), as a delightfully potent breakthrough technology, not only identifies various cell types and response states but also uncovers cell-specific gene expression changes, gene regulatory networks, intercellular communication, and cellular movement trajectories, among others, in different cell types. In this review, we describe the technology of scRNA-seq in detail and discuss and summarize the application of scRNA-seq in exploring neurological diseases, elaborating the corresponding specific mechanisms of the diseases as well as providing a reliable basis for new therapeutic approaches. Finally, we affirm that scRNA-seq promotes the development of the neuroscience field and enables us to have a deeper cellular understanding of neurological diseases in the future, which provides strong support for the treatment of neurological diseases and the improvement of patients' prognosis.

Deciphering cutaneous melanoma prognosis through LDL metabolism: Single-cell transcriptomics analysis via 101 machine learning algorithms

Cutaneous melanoma poses a formidable challenge within the field of oncology, marked by its aggressive nature and capacity for metastasis. Despite extensive research uncovering numerous genetic and molecular contributors to cutaneous melanoma development, there remains a critical knowledge gap concerning the role of lipids, notably low-density lipoprotein (LDL), in this lethal skin cancer. This article endeavours to bridge this knowledge gap by delving into the intricate interplay between LDL metabolism and cutaneous melanoma, shedding light on how lipids influence tumour progression, immune responses and potential therapeutic avenues. Genes associated with LDL metabolism were extracted from the GSEA database. We acquired and analysed single-cell sequencing data (GSE215120) and bulk-RNA sequencing data, including the TCGA data set, GSE19234, GSE22153 and GSE65904. Our analysis unveiled the heterogeneity of LDL across various cell types at the single-cell sequencing level. Additionally, we constructed an LDL-related signature (LRS) using machine learning algorithms, incorporating differentially expressed genes and highly correlated genes. The LRS serves as a valuable tool for assessing the prognosis, immunity and mutation status of patients with cutaneous melanoma. Furthermore, we conducted experiments on A375 and WM-115 cells to validate the function of PPP2R1A, a pivotal gene within the LRS. Our comprehensive approach, combining advanced bioinformatics analyses with an extensive review of current literature, presents compelling evidence regarding the significance of LDL within the cutaneous melanoma microenvironment.

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.

Application of Single Cell Gene Expression Technologies to Neurotoxicology

Neurotoxicological research faces the challenge of linking biological changes resulting from exposures to neuronal function. An additional challenge is understanding cell-type specific differences and selective vulnerabilities of distinct neuronal populations to toxic insults. Single cell RNA-sequencing (scRNA-seq) allows for measurement of the transcriptome of individual cells. This makes it a valuable tool for validating and characterizing cell types present in multicell type samples in complex tissue or cell culture models, but also for understanding how different cell types respond to toxic insults. Pathway analysis of differentially expressed genes can provide in depth insights into underlying cell type-specific mechanisms of neurotoxicity. Toxicological data often has to be translated to outcomes for human health which requires an understanding of inter-species differences. Transcriptomic data aids in understanding these differences, including understanding developmental timelines of different species. We believe that scRNA-seq holds exciting promises for future neurotoxicological research.

Single-cell RNA sequencing to understand host-virus interactions

Single-cell RNA sequencing (scRNA-seq) has allowed for the profiling of host and virus transcripts and host-virus interactions at single-cell resolution. This review summarizes the existing scRNA-seq technologies together with their strengths and weaknesses. The applications of scRNA-seq in various virological studies are discussed in depth, which broaden the understanding of the immune atlas, host-virus interactions, and immune repertoire. scRNA-seq can be widely used for virology in the near future to better understand the pathogenic mechanisms and discover more effective therapeutic strategies.

Urinary SPP1 has potential as a non-invasive diagnostic marker for focal segmental glomerulosclerosis

Focal segmental glomerulosclerosis (FSGS) is a type of chronic glomerular nephropathy showing characteristic glomerular sclerosis, diagnosed by kidney biopsy. However, it is difficult and expensive to monitor disease progression with repeated renal biopsy in clinical practice, and thus here we explored the feasibility of urine biomarkers as non-invasive diagnostic tools. We downloaded scRNA-seq datasets of 20 urine cell samples and 3 kidney tissues and obtained two gene lists encoding extracellular proteins for bioinformatic analysis; in addition, we identified key EP-Genes by immunohistochemical staining and performed bulk RNA sequencing with 12 urine samples. We report that urine cells and kidney cells were correlated. A total of 64 EP-Genes were acquired by intersecting genes of distal tubular cluster with extracellular proteins. Function enrichment analysis showed that EP-Genes might be involved in the immune response and extracellular components. Six key EP-Genes were identified and correlated with renal function. IMC showed that key EP-Genes were located mainly in tubules. Cross verification and examination of a urine RNAseq dataset showed that SPP1 had diagnostic potential for FSGS. The presence of urine SPP1 was primarily associated with macrophage infiltration in kidney, and the pathogenesis of FSGS may be related to innate immunity. Urinary cells seemed to be strongly similar to kidney cells. In summary, SPP1 levels reflect renal function and may have potential as a biomarker for non-invasive diagnosis of FSGS.

Unified Mouse and Human Kidney Single-Cell Expression Atlas Reveal Commonalities and Differences in Disease States

SIGNIFICANCE STATEMENT

Mouse models have been widely used to understand kidney disease pathomechanisms and play an important role in drug discovery. However, these models have not been systematically analyzed and compared. The authors characterized 18 different mouse kidney disease models at both bulk and single-cell gene expression levels and compared single-cell gene expression data from diabetic kidney disease (DKD) mice and from patients with DKD. Although single cell-level gene expression changes were mostly model-specific, different disease models showed similar changes when compared at a pathway level. The authors also found that changes in fractions of cell types are major drivers of bulk gene expression differences. Although the authors found only a small overlap of single cell-level gene expression changes between the mouse DKD model and patients, they observed consistent pathway-level changes.

BACKGROUND

Mouse models have been widely used to understand kidney disease pathomechanisms and play an important role in drug discovery. However, these models have not been systematically analyzed and compared.

METHODS

We analyzed single-cell RNA sequencing data (36 samples) and bulk gene expression data (42 samples) from 18 commonly used mouse kidney disease models. We compared single-nucleus RNA sequencing data from a mouse diabetic kidney disease model with data from patients with diabetic kidney disease and healthy controls.

RESULTS

We generated a uniformly processed mouse single-cell atlas containing information for nearly 300,000 cells, identifying all major kidney cell types and states. Our analysis revealed that changes in fractions of cell types are major drivers of differences in bulk gene expression. Although gene expression changes at the single-cell level were mostly model-specific, different disease models showed similar changes when compared at a pathway level. Tensor decomposition analysis highlighted the important changes in proximal tubule cells in disease states. Specifically, we identified important alterations in expression of metabolic and inflammation-associated pathways. The mouse diabetic kidney disease model and patients with diabetic kidney disease shared only a small number of conserved cell type-specific differentially expressed genes, but we observed pathway-level activation patterns conserved between mouse and human diabetic kidney disease samples.

CONCLUSIONS

This study provides a comprehensive mouse kidney single-cell atlas and defines gene expression commonalities and differences in disease states in mice. The results highlight the key role of cell heterogeneity in driving changes in bulk gene expression and the limited overlap of single-cell gene expression changes between animal models and patients, but they also reveal consistent pathway-level changes.

Dual-GCN-based deep clustering with triplet contrast for ScRNA-seq data analysis

Single-cell RNA sequencing (ScRNA-seq) technology reveals gene expression information at the cellular level. The critical tasks in ScRNA-seq data analysis are clustering and dimensionality reduction. Recent deep clustering algorithms are used to optimize the two tasks jointly, and their variations, graph-based deep clustering algorithms, are used to capture and preserve topological information in the process. However, the existing graph-based deep clustering algorithms ignore the distribution information of nodes when constructing cell graphs which leads to incomplete information in the embedding representation; and graph convolutional networks (GCN), which are most commonly used, often suffer from over-smoothing that leads to high sample similarity in the embedding representation and then poor clustering performance. Here, the dual-GCN-based deep clustering with Triplet contrast (scDGDC) is proposed for dimensionality reduction and clustering of scRNA-seq data. Two critical components are dual-GCN-based encoder for capturing more comprehensive topological information and triplet contrast for reducing GCN over-smoothing. The two components improve the dimensionality reduction and clustering performance of scDGDC in terms of information acquisition and model optimization, respectively. The experiments on eight real ScRNA-seq datasets showed that scDGDC achieves excellent performance for both clustering and dimensionality reduction tasks and is high robustness to parameters.

Bovine milk-derived cells express transcriptome markers of pluripotency and secrete bioactive factors with regenerative and antimicrobial activity

The bovine mammary stem/progenitor cell secretome stimulates regeneration in vitro and contains proteins associated with antimicrobial defense. This has led to the exploration of the secretome as a biologic treatment for mastitis, a costly inflammation of the udder commonly caused by bacteria. This study reports on a population of bovine mammary stem/progenitor cells isolated non-invasively from milk (MiDCs). MiDCs were characterized by immunophenotyping, mammosphere formation assays, and single cell RNA sequencing. They displayed epithelial morphology, exhibited markers of mammary stem/progenitor cells, and formed mammospheres, like mammary gland tissue-isolated stem/progenitor cells. Single cell RNA sequencing revealed two sub-populations of MiDCs: epithelial cells and macrophages. Functionally, the MiDC secretome increased fibroblast migration, promoted angiogenesis of endothelial cells, and inhibited the growth of mastitis-associated bacteria, including antibiotic-resistant strains, in vitro. These qualities of MiDCs render them a source of stem cells and stem cell products that may be used to treat diseases affecting the dairy industry, including mastitis.

Integration single-cell and bulk RNA-sequencing data to reveal senescence gene expression profiles in heart failure

Background

Heart failure (HF) represents one of healthcare's biggest challenges. Although rarely noticed, aging is a crucial risk factor for cardiovascular disease. Our study aims to reveal aging's role in HF by integrating single-cell RNA-sequencing (scRNA-seq) and bulk RNA-sequencing databases.

Methods

We collected HF heart sample data from the Gene Expression Omnibus database and senescence gene data from CellAge. The FindCluster () package was used for cell cluster analysis. Differentially expressed genes (DEG) were identified operating the FindMarkers function. Cell activity score calculation was performed using the AUCell package. UpSetR plotted the intersection between DEGs of active cell types, bulk data DEGs, and genes associated with aging. Using the DGIdb database gene-drug interaction data, we search for potential targeted therapeutics based on common senescence genes.

Results

The scRNA-seq data revealed myocardial heterogeneity in HF tissues. A series of crucial common senescence genes were found. The senescence gene expression profile hints at an intriguing connection between monocytes and HF. After analyzing the DEGs in the bulk dataset, the DEGs in scRNA-seq, the DEGs in each active cell type, and senescence genes, we identified ten genes as common senescence genes present in HF. Correlation analysis of transcriptomics, proteomics, and ceRNA was performed to provide ideas for future studies individually. Moreover, we discovered that common senescence genes and potential therapeutic drugs interact among different cell types. Further research is needed on the expression pattern of senescence genes and molecular regulation in HF.

Conclusions

In summary, we identified the functional significance of the senescence gene in HF using integrated data. It is possible that this more profound understanding of how senescence contributes to the development of HF will aid in unraveling the mechanisms that promote the disease and provide hints for developing therapeutics.

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.

Aberrant NAD synthetic flux in podocytes under diabetic conditions and effects of indoleamine 2,3-dioxygenase on promoting de novo NAD synthesis

Nicotinamide adenine dinucleotide (NAD) is an essential coenzyme in the kidney. The first step in de novo NAD synthesis is regulated by indoleamine 2,3-dioxygenase (IDO), a tryptophan-catabolizing enzyme. Here, we investigated NAD synthetic flux and NAD levels in podocytes under diabetic conditions. We also studied the effects of IDO overexpression on NAD synthetic flux and high glucose (HG)-induced podocyte injury. NAD synthetases in the de novo, Preiss-Handler and salvage pathways were analyzed using in vivo single-nucleus RNA sequencing datasets (GSE131882) of control and diabetic kidney disease (DKD). The mRNA levels of these NAD synthetases were measured in vitro in HG-treated podocytes. The effects of IDO on NAD synthesis were examined by transducing cultured podocytes with an adenovirus encoding IDO, and apoptosis, podocyte markers and mobility were investigated. Cellular transcriptome analysis revealed that control podocytes had relatively low levels of NAD synthetases. In DKD podocytes, de novo NAD synthetase levels were further downregulated. IDO levels were virtually undetectable and did not increase in DKD. In vitro experiments confirmed aberrant de novo NAD synthetic flux and decreased IDO levels in HG-treated podocytes. Overexpression of IDO promoted NAD de novo synthesis, reduced NAD-bypass metabolic enzyme, increased NAD content and recovered podocyte injury markers under diabetic conditions. Taken together, our findings suggest that the de novo NAD synthetic flux is aberrant in DKD, and IDO promotes de novo NAD synthesis and NAD levels, as well as alleviates injury in HG-treated podocytes.

Single-cell analysis highlights differences in druggable pathways underlying adaptive or fibrotic kidney regeneration

The kidney has tremendous capacity to repair after acute injury, however, pathways guiding adaptive and fibrotic repair are poorly understood. We developed a model of adaptive and fibrotic kidney regeneration by titrating ischemic injury dose. We performed detailed biochemical and histological analysis and profiled transcriptomic changes at bulk and single-cell level (> 110,000 cells) over time. Our analysis highlights kidney proximal tubule cells as key susceptible cells to injury. Adaptive proximal tubule repair correlated with fatty acid oxidation and oxidative phosphorylation. We identify a specific maladaptive/profibrotic proximal tubule cluster after long ischemia, which expresses proinflammatory and profibrotic cytokines and myeloid cell chemotactic factors. Druggability analysis highlights pyroptosis/ferroptosis as vulnerable pathways in these profibrotic cells. Pharmacological targeting of pyroptosis/ferroptosis in vivo pushed cells towards adaptive repair and ameliorates fibrosis. In summary, our single-cell analysis defines key differences in adaptive and fibrotic repair and identifies druggable pathways for pharmacological intervention to prevent kidney fibrosis.

New insights from the single-cell level: Tumor associated macrophages heterogeneity and personalized therapy

Tumor-associated macrophages (TAMs) are important immune cells in the tumor microenvironment, and their invasion in tumors is closely related to poor prognosis. Although TAMs are recognized as therapeutic targets, their heterogeneity makes studying tumor mechanism and developing drugs targeting TAMs difficult. The study of TAMs heterogeneity can be used to analyze the mechanism of tumor progression and drug resistance, and may provide possible treatment strategies for cancer patients. Single-cell RNA sequencing (scRNA-seq) can reveal the RNA expression profile for each TAM to distinguish heterogeneity, thereby providing a more efficient detection method and more accurate information for TAM-related studies. In this review, by summarizing the research progress in macrophage heterogeneity and other aspects of scRNA-seq over the past five years, we introduced the development of scRNA-seq technology and its application status in solid tumors, analyzed the advantages and selections of scRNA-seq in TAMs, and summarized the detailed specific research fields. To explore the mechanism of tumor progression and drug intervention from single cell level will provide new perspective for personalized treatment strategies targeting macrophages.

Microfluidics Facilitates the Development of Single-Cell RNA Sequencing

Single-cell RNA sequencing (scRNA-seq) technology provides a powerful tool for understanding complex biosystems at the single-cell and single-molecule level. The past decade has been a golden period for the development of single-cell sequencing, with scRNA-seq undergoing a tremendous leap in sensitivity and throughput. The application of droplet- and microwell-based microfluidics in scRNA-seq has contributed greatly to improving sequencing throughput. This review introduces the history of development and important technical factors of scRNA-seq. We mainly focus on the role of microfluidics in facilitating the development of scRNA-seq technology. To end, we discuss the future directions for scRNA-seq.

Ceramides and Acute Kidney Injury

Altered lipid metabolism is a characteristic feature and potential driving factor of acute kidney injury (AKI). Of the lipids that accumulate in injured renal tissues, ceramides are potent regulators of metabolism and cell fate. Up-regulation of ceramide synthesis is a common feature shared across several AKI etiologies in vitro and in vivo. Furthermore, ceramide accumulation is an early event in the natural history of AKI that precedes cell death and organ dysfunction. Emerging evidence suggests that inhibition of ceramide accumulation may improve renal outcomes in several models of AKI. This review examines the landscape of ceramide metabolism and regulation in the healthy and injured kidney. Furthermore, we discuss the body of literature regarding ceramides as therapeutic targets for AKI and consider potential mechanisms by which ceramides drive kidney pathogenesis.

Single cell regulatory landscape of the mouse kidney highlights cellular differentiation programs and disease targets

Determining the epigenetic program that generates unique cell types in the kidney is critical for understanding cell-type heterogeneity during tissue homeostasis and injury response. Here, we profile open chromatin and gene expression in developing and adult mouse kidneys at single cell resolution. We show critical reliance of gene expression on distal regulatory elements (enhancers). We reveal key cell type-specific transcription factors and major gene-regulatory circuits for kidney cells. Dynamic chromatin and expression changes during nephron progenitor differentiation demonstrates that podocyte commitment occurs early and is associated with sustained Foxl1 expression. Renal tubule cells follow a more complex differentiation, where Hfn4a is associated with proximal and Tfap2b with distal fate. Mapping single nucleotide variants associated with human kidney disease implicates critical cell types, developmental stages, genes, and regulatory mechanisms. The single cell multi-omics atlas reveals key chromatin remodeling events and gene expression dynamics associated with kidney development.