Command line Rscript to perform differential expression analysis using DESeq2, edgeR and limma-voom quick and reproducible. Main use is for human samples, but this could easily get adapted if desired.
- Performs counting using featureCounts (bam input)
- Can handle salmon quantification
- Allows specification of covariates
- Rigorous checking of input data
- Creates various plots
- Creates detailed tables and lists of differentially expressed genes
DEA.R [sampleInfoFile] [annotation.gtf]
with:
- sampleInfoFile: file specifying samples and covariates (see further)
- annotation.gtf: An annotation file in gtf format matching the reference genome used for alignment
The sample info file contains all information required for the script.
These fields are mandatory present in the file, header is case sensitive
- 'file': exact path to bamfile or quant.sf to be used for counting
- 'sample': identifier of the sample to be used
- 'condition': main factor on which the differential expression should be performed.
The reference level for the field condition should be "CON" (case sensitive).
If more than 2 groups are present in this field all contrasts with the "CON" group will be tested. - 'sequencing': containing mode of sequencing: either PE or SE
- 'strandedness' type of strandedness of data: either unstranded, stranded or reverse
- 'gender': with possible values "m", "f" and "u"
can be specified in the sample info file, e.g.:
- 'library_prep_by': Bob, Alice
This tool comes with a bash script "MakeTest.sh" which will download data from SRA, perform alignment and create a sample info file which you can use directly for testing the DEA.R script. Only one million reads per sample are downloaded to reduce time and storage space used.
The data used is from RNA-seq on airway smooth muscle cells treated with dexamethasone treated with dexamethasone,
- Citation: Himes BE, Jiang X, Wagner P, Hu R, Wang Q, Klanderman B, Whitaker RM, Duan Q, Lasky-Su J, Nikolos C, Jester W, Johnson M, Panettieri R Jr, Tantisira KG, Weiss ST, Lu Q. 'RNA-Seq Transcriptome Profiling Identifies CRISPLD2 as a Glucocorticoid Responsive Gene that Modulates Cytokine Function in Airway Smooth Muscle Cells.' PLoS One. 2014 Jun 13;9(6):e99625.
- PMID: 24926665.
- GEO: GSE52778.
Dependencies for this bash script are fastq-dump, STAR, wget and gunzip
It will use about 33 Gb of space on your hard disk, so make sure sufficient space is left.
The only argument this script takes is the number of processes which it can use,
e.g. for running twelve processes you use bash MakeTest.sh 12
Make sure to adjust this accordingly to your available computer architecture
Command line interface to Enrichr. I am not affiliated to the authors of Enrichr but just use their API. Input is a list of gene symbols, either using a file or on stdin.
enrichr-cli.py [-h] [-g GENES] [-d [DATABASES [DATABASES ...]]] [-w] [-p PREFIX] [-o OUTDIR]
mandatory argument:
-g GENES, --genes GENES
A genelist to be queried, either a list in a file or '-' for a list of genes on stdin
optional arguments:
-h, --help show this help message and exit
-d [DATABASES [DATABASES ...]], --databases [DATABASES [DATABASES ...]]
Databases to query, omit to use a default set.
-w, --which List databases which can be queried and quit.
-p PREFIX, --prefix PREFIX
Fixed prefix to name the output files with
-o OUTDIR, --outdir OUTDIR
Output directory to store files in. Will be created if it doesn't exist.
Script to create violin plots of expression data
SimpleViolin.R countsfile.txt sample_info_file.txt targets.
Targets are either
- A comma separated list of gene symbols
- A differential expression result, of which the top 10 will be taken