Alternative DE methods on LUHMES CROP-seq Data

-- In-house scripts

This tutorial describes how to run alternative differential expression (DE) methods on LUHMES CROP-seq data. All file paths here are relative to this working directory: /project2/xinhe/yifan/Factor_analysis/LUHMES/.

0.1 DESeq2

0.1.1 R script

Location: run_deseq.R;

Command line usage example:

Rscript run_deseq.R --out_folder=processed_data/DESeq2/dev_top6k_negctrl/ --marker=ADNP --NTC_only=TRUE --permute=FALSE

marker: the gRNA perturbation condition to test the differential effects,
NTC_only: whether to compare with only nontargeting control cells or all cells in other conditions.

0.1.2 Run jobs at scale

Edit the Python script submit_DE_jobs_in_batches.py, which can submit sbatch jobs for all gRNA targets.

Make sure to modify the sbatch script defined in make_DEG_sbatch() and loop over all conditions.

To submit jobs, run the edited script: ./submit_DE_jobs_in_batches.py.

0.1.3 Output

Each gRNA condition vs negative control: processed_data/DESeq2/dev_top6k_negctrl/;

Each gRNA condition vs the rest of cells: processed_data/DESeq2/dev_top6k/;

0.2 MAST

Location: run_MASTcpmDetRate.R;

Command line usage example:

Rscript run_MASTcpmDetRate.R --out_folder=processed_data/MAST/dev_top6k_negctrl/ --marker=ADNP --NTC_only=TRUE --permute=FALSE

marker: the gRNA perturbation condition to test the differential effects,
NTC_only: whether to compare with only nontargeting control cells or all cells in other conditions.

0.2.1 Run jobs at scale

Edit the Python script submit_DE_jobs_in_batches.py, which can submit sbatch jobs for all gRNA targets.

Make sure to modify the sbatch script defined in make_DEG_sbatch() and loop over all conditions.

To submit jobs, run the edited script: ./submit_DE_jobs_in_batches.py.

0.2.2 Output

Each gRNA condition vs negative control: processed_data/MAST/dev_top6k_negctrl/;

Each gRNA condition vs the rest of cells: processed_data/MAST/dev_top6k/;

0.3 scMAGeCK

0.3.1 R script

Location: scmageck/run_scmageck.R;

Command line usage example (works on all gRNA conditions at once):

Rscript scmageck/run_scmageck.R --permute=FALSE

0.3.2 Output

Each gRNA condition vs negative control: scmageck/scmageck_lr.LUHMES.dev_res_top_6k.rds.