We are designing a CRISPR/Cas9 experiment and thinking of the down-stream data analyses.

Are there any R packages for analysing raw NGS read count data from pooled genetic screens using CRIPSR/Cas9 to disrupt gene expression in a population of cells? I guess we will need to start with basics, i.e. sequences processing, data exploration, visualisation etc.

  • 3
    $\begingroup$ How is this different from other expression analysis? I.e. why can’t you use one of the standard workflows for differential expression analysis? $\endgroup$ – Konrad Rudolph Jun 5 '17 at 17:26
  • $\begingroup$ I would just like to highlight that using CRISPR/Cas9 for editing and just estimating the expression might not be a great idea per se. If you have really a lot number of samples you might be interested first to check the mutational burden introduced by this editing. Then preferably go with expression analysis. I would go with the below tools as mentioned as well but will not do anything without even considering if my editing is throwing a lot of SNVs or not. If its not throwing then you should be fine doing the downstream analysis with linear mixed effect. $\endgroup$ – ivivek_ngs Jun 6 '17 at 14:19

Assuming you mean CRISPR screens targeting many loci (for example using the GeCKO library) there is an R package here.

It uses a linear mixed effect model to compare guide counts before and after a selection step, allowing for multiple guides/gene and multiple replicates. It can also do the initial read alignment AFAIK and the author should be quite responsive in case you have questions.

Another option could be the python program MAGeCK


From my experience, there is few about CRISPR pooled screenings in R. One of them is ScreenBEAM. It uses a linear mixture model to determine significance at gene level by considering all the guides together at the same step in the analysis. However, it does not work very well for me.

There are other ways of analyze these screenings outside of R that are easy to use such the suggested MAGeCK or others such HiTSelect, BAGEL or CRISPR-Analyzer.

CRISPR-Analyzer is an online tool that could be useful and that recapitulates some of the mentioned methods in a single analysis. However, if you want to work from command line I would recommend you BAGEL for gene essentiality experiments and MAGeCK for other experimental setup (i.e, control vs treatment)

Sorry for not putting all the links but I can only put two at the moment!


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.