Small tutorial of WITER for test and annotate cancer-driver genes

Miaoxin Li (limx54@163.com)

 

Reference: https://pmglab.top/kggseq/doc10/UserManual.html

Input data:

1.      A somatic variant summary file of breast cancer  [compiled from tumor portal http://www.tumorportal.org/tumor_types?ttype=BRCA]

examples/hg19_BreastAdenocarcinoma.maf

Purpose: Identify cancer-driver genes of breast cancer 


Run the commands step by step to see what will happen

1.      Annotate sequence variants by gene models from RefGenes and GENCODE:

java -Xmx6g -jar kggseq.jar --maf-file examples/hg19_BreastAdenocarcinoma.maf --out test1 --excel --db-gene refgene,gencode --gene-feature-in 0,1,2,3,4,5,6,7

2.      Detect cancer-driver genes based on somatic mutation rates at non-synonymous variants by ITER

java -Xmx6g -jar kggseq.jar --maf-file examples/hg19_BreastAdenocarcinoma.maf --out test1 --excel --db-gene refgene,gencode --gene-feature-in 0,1,2,3,4,5,6,7 --qqplot --iter-gene-coding

 

3.      Detect cancer-driver genes based on somatic mutation rates at non-synonymous variants by WITER

java -Xmx6g -jar kggseq.jar --maf-file examples/hg19_BreastAdenocarcinoma.maf --out test1 --excel --db-gene refgene,gencode --gene-feature-in 0,1,2,3,4,5,6,7 --qqplot --witer-gene-coding --db-score dbnsfp --cancer-mut-predict --witer-gene-coding

 

4.      Detect cancer-driver genes based on somatic mutation rates at non-synonymous variants by WITER with a reference for background gene mutation

java -Xmx6g -jar kggseq.jar --maf-file examples/hg19_BreastAdenocarcinoma.maf --out test1 --excel --db-gene refgene,gencode --gene-feature-in 0,1,2,3,4,5,6,7 --qqplot --db-score dbnsfp --cancer-mut-predict --witer-gene-coding --ref-mut-file examples/Ovarian.gene.mutationburden.txt

5.      Detect cancer-driver genes based on somatic mutation rates at non-synonymous variants by WITER and annotate variants COSMIC

java -Xmx6g -jar kggseq.jar --maf-file examples/hg19_BreastAdenocarcinoma.maf --out test1 --excel --db-gene refgene,gencode --gene-feature-in 0,1,2,3,4,5,6,7 --qqplot --db-score dbnsfp --cancer-mut-predict --witer-gene-coding --cosmic-annot

 

6.      Detect cancer-driver genes based on somatic mutation rates at synonymous variants by ITER

java -Xmx6g -jar kggseq.jar --maf-file examples/hg19_BreastAdenocarcinoma.maf --out test1 --excel --db-gene refgene,gencode --gene-feature-in 0,1,2,3,4,5,6,7 --qqplot --response-gene-feature 7 --explanatory-gene-feature - --iter-gene-coding