search

aertslab / SCENIC

SCENIC is an R package to infer Gene Regulatory Networks and cell types from single-cell RNA-seq data.
http://scenic.aertslab.org
GNU General Public License v3.0
420 stars 95 forks source link

calculate TF using the differential expression genes #162

Closed wtj505531480 closed 4 years ago

wtj505531480 commented 4 years ago

Thanks for the SCENIC, could I use the differential expression genes as the input of the SCENIC?thank you

s-aibar commented 4 years ago

Yes, there are two main options that I can think of:

A) You could filter the input matrix to include only the differentially expressed genes (but beware that that will likely reduce a lot the sensitivity of the networks, see the Nature Methods 2017 paper).

B) If you are only intererested on potential regulators, you can just use RcisTarget on the differentially expressed genes. That will provide potential TFs/motifs as regulators for those genes (and therefore, the cell type or state). That analysis might already be enough, but if you want to get a more complete "SCENIC-like" analysis, you also check which of those genes are regulated by each TF (you can see the code in the detailed vignette for Step 2). That will be kind of a manual regulon/network, and of course, you can then score it on your cells with AUCell (detailed vignette for step 3).

wtj505531480 commented 4 years ago

Thanks for your reply !!