noamteyssier / crispr_screen

Efficient Differential Expression and Gene Aggregation for CRISPR screens.
https://noamteyssier.github.io/crispr_screen/
MIT License
6 stars 0 forks source link

CRISPR Screen

MIT licensed actions status codecov DOI

This tool is recreation of the MAGeCK Differential Expression algorithm described in the original paper with options included. It also extends the algorithm and provides an option for the MAGeCK-INC pipeline which performs a Mann-Whitney U-Test on the Non-Targeting Controls for gene score aggregation as opposed to the traditional alpha-RRA algorithm. The goal of this project is to have an efficient and clean implementation of a CRISPR-screening platform that is as easy to use as possible without sacrificing efficiency and interpretability.

Installation

Installing Rust

If you do not have the rust toolchain installed you can do so with the following command

curl --proto '=https' --tlsv1.2 -sSf https://sh.rustup.rs | sh

Installing CRISPR Screen

cargo install crispr_screen

Usage

This tool is meant to run from the commandline and expects at minimum 3 arguments:

  1. an sgRNA-Gene count table
  2. the control label(s)
  3. the treatment label(s).

The following list of examples is not exhaustive and if you are confused by any parameters I recommend running the help menu:

crispr_screen --help

Basic Run with a Single Control and a Single Treatment

Here is the command layout for a basic run with only these arguments:

crispr_screen test -i <count_table> -c <control_label> -t <treatment_label>

Basic Run with a Single Control and Multiple Treatments

If you have multiple treatment and/or control labels you can provide multiple arguments to the same flag. In this case I am providing two treatment labels for the same -t flag.

crispr_screen test -i <count_table> -c <control_label> -t <treatment_label_a> <treatment_label_b>

Basic Run with a Specified Output Prefix

By default this will write the sgRNA-level and Gene-level results to results.sgrna_results.tab and results.gene_results.tab respectively, but you can specificy your own output prefix (i.e. replacing results) with the -o flag.

crispr_screen test -i <count_table> -c <control_label> -t <treatment_label> -o my_prefix

Alternative Normalization Strategies

Currently there are two normalization strategies provided by this tool - MedianRatio and Total. You can specify which one you'd like to use with the -n flag. The two accepted options are median-ratio and total.

crispr_screen test -i <count_table> -c <control_label> -t <treatment_label> -n total

Alternative Gene-Aggregation Strategies

Currently there are two gene aggregation strategies provided by this tool - AlphaRRA and INC. By default the gene aggregation strategy run is AlphaRRA. However, you can specify which one you'd like to use with the -g flag. The two accepted options are rra and inc.

crispr_screen test -i <count_table> -c <control_label> -t <treatment_label> -g inc

References