30 releases
0.3.19 | Oct 30, 2024 |
---|---|
0.3.13 | Sep 10, 2024 |
0.2.8 | Jul 18, 2024 |
0.2.6 | Mar 25, 2024 |
0.1.14 | Feb 15, 2023 |
#102 in Biology
39 downloads per month
325KB
4.5K
SLoC
CRISPR Screen
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:
- an sgRNA-Gene count table
- the control label(s)
- 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
Dependencies
~37–70MB
~1M SLoC