scDown: A Pipeline for Single-Cell RNA-Seq Downstream Analysis

Abstract

Single-cell transcriptomics data are analyzed using two popular tools, Seurat and Scanpy. Multiple separate tools are used downstream of Seurat and Scanpy cell annotation to study cell differentiation and communication, including cell proportion difference analysis between conditions, pseudotime and trajectory analyses to study cell transition, and cell–cell communication analysis. To automate the integrative cell differentiation and communication analyses of single-cell RNA-seq data, we developed a single-cell RNA-seq downstream analysis pipeline called “scDown”. This R package includes cell proportion difference analysis, cell–cell communication analysis, pseudotime analysis, and RNA velocity analysis. Both Seurat and Scanpy annotated single-cell RNA-seq data are accepted in this pipeline. We applied scDown to a published dataset and identified a unique, previously undiscovered signature of neuronal inflammatory signaling associated with a rare genetic neurodevelopmental disorder. These findings were not identified with a simple implementation of Seurat differential gene expression analysis, illustrating the value of our pipeline in biological discovery. scDown can be broadly utilized in downstream analyses of scRNA-seq data, particularly in rare diseases.