A scalable SCENIC workflow for single-cell gene regulatory network analysis
Top Cited Papers
- 19 June 2020
- journal article
- research article
- Published by Springer Science and Business Media LLC in Nature Protocols
- Vol. 15 (7), 2247-2276
- https://doi.org/10.1038/s41596-020-0336-2
Abstract
This protocol explains how to perform a fast SCENIC analysis alongside standard best practices steps on single-cell RNA-sequencing data using software containers and Nextflow pipelines. SCENIC reconstructs regulons (i.e., transcription factors and their target genes) assesses the activity of these discovered regulons in individual cells and uses these cellular activity patterns to find meaningful clusters of cells. Here we present an improved version of SCENIC with several advances. SCENIC has been refactored and reimplemented in Python (pySCENIC), resulting in a tenfold increase in speed, and has been packaged into containers for ease of use. It is now also possible to use epigenomic track databases, as well as motifs, to refine regulons. In this protocol, we explain the different steps of SCENIC: the workflow starts from the count matrix depicting the gene abundances for all cells and consists of three stages. First, coexpression modules are inferred using a regression per-target approach (GRNBoost2). Next, the indirect targets are pruned from these modules using cis-regulatory motif discovery (cisTarget). Lastly, the activity of these regulons is quantified via an enrichment score for the regulon’s target genes (AUCell). Nonlinear projection methods can be used to display visual groupings of cells based on the cellular activity patterns of these regulons. The results can be exported as a loom file and visualized in the SCope web application. This protocol is illustrated on two use cases: a peripheral blood mononuclear cell data set and a panel of single-cell RNA-sequencing cancer experiments. For a data set of 10,000 genes and 50,000 cells, the pipeline runs in <2 h.This publication has 52 references indexed in Scilit:
- STAR: ultrafast universal RNA-seq alignerBioinformatics, 2012
- An integrated encyclopedia of DNA elements in the human genomeNature, 2012
- i-cisTarget: an integrative genomics method for the prediction of regulatory features and cis-regulatory modulesNucleic Acids Research, 2012
- Inferring Regulatory Networks from Expression Data Using Tree-Based MethodsPLOS ONE, 2010
- A global network of transcription factors, involving E2A, EBF1 and Foxo1, that orchestrates B cell fateNature Immunology, 2010
- EnsemblCompara GeneTrees: Complete, duplication-aware phylogenetic trees in vertebratesGenome Research, 2008
- Quantifying similarity between motifsGenome Biology, 2007
- Gene prioritization through genomic data fusionNature Biotechnology, 2006
- Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profilesProceedings of the National Academy of Sciences of the United States of America, 2005
- Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction NetworksGenome Research, 2003