Repository logo
 

A step-by-step workflow for low-level analysis of single-cell RNA-seq data with Bioconductor.

Published version
Peer-reviewed

Type

Article

Change log

Authors

Lun, Aaron TL 
McCarthy, Davis J 
Marioni, John C 

Abstract

Single-cell RNA sequencing (scRNA-seq) is widely used to profile the transcriptome of individual cells. This provides biological resolution that cannot be matched by bulk RNA sequencing, at the cost of increased technical noise and data complexity. The differences between scRNA-seq and bulk RNA-seq data mean that the analysis of the former cannot be performed by recycling bioinformatics pipelines for the latter. Rather, dedicated single-cell methods are required at various steps to exploit the cellular resolution while accounting for technical noise. This article describes a computational workflow for low-level analyses of scRNA-seq data, based primarily on software packages from the open-source Bioconductor project. It covers basic steps including quality control, data exploration and normalization, as well as more complex procedures such as cell cycle phase assignment, identification of highly variable and correlated genes, clustering into subpopulations and marker gene detection. Analyses were demonstrated on gene-level count data from several publicly available datasets involving haematopoietic stem cells, brain-derived cells, T-helper cells and mouse embryonic stem cells. This will provide a range of usage scenarios from which readers can construct their own analysis pipelines.

Description

Keywords

Bioconductor, RNA-seq, Single cell, bioinformatics, workflow

Journal Title

F1000Res

Conference Name

Journal ISSN

2046-1402
1759-796X

Volume Title

5

Publisher

F1000 Research Ltd