oncodrive: Detect cancer driver genes based on positional clustering of...

Description Usage Arguments Details Value References See Also Examples

View source: R/oncodrive.R

Description

Clusters variants based on their position to detect disease causing genes.

Usage

1
2
3
4
5
6
7
8
9
oncodrive(
  maf,
  AACol = NULL,
  minMut = 5,
  pvalMethod = "zscore",
  nBgGenes = 100,
  bgEstimate = TRUE,
  ignoreGenes = NULL
)

Arguments

maf

an MAF object generated by read.maf

AACol

manually specify column name for amino acid changes. Default looks for field 'AAChange'

minMut

minimum number of mutations required for a gene to be included in analysis. Default 5.

pvalMethod

either zscore (default method for oncodriveCLUST), poisson or combined (uses lowest of the two pvalues).

nBgGenes

minimum number of genes required to estimate background score. Default 100. Do not change this unless its necessary.

bgEstimate

If FALSE skips background estimation from synonymous variants and uses predifined values estimated from COSMIC synonymous variants.

ignoreGenes

Ignore these genes from analysis. Default NULL. Helpful in case data contains large number of variants belonging to polymorphic genes such as mucins and TTN.

Details

This is the re-implimentation of algorithm defined in OncodriveCLUST article. Concept is based on the fact that most of the variants in cancer causing genes are enriched at few specific loci (aka hotspots). This method takes advantage of such positions to identify cancer genes. Cluster score of 1 means, a single hotspot hosts all observed variants. If you use this function, please cite OncodriveCLUST article.

Value

data table of genes ordered according to p-values.

References

Tamborero D, Gonzalez-Perez A and Lopez-Bigas N. OncodriveCLUST: exploiting the positional clustering of somatic mutations to identify cancer genes. Bioinformatics. 2013; doi: 10.1093/bioinformatics/btt395s

See Also

plotOncodrive

Examples

1
2
3
laml.maf <- system.file("extdata", "tcga_laml.maf.gz", package = "maftools")
laml <- read.maf(maf = laml.maf)
laml.sig <- oncodrive(maf = laml, AACol = 'Protein_Change', minMut = 5)

maftools documentation built on Feb. 6, 2021, 2 a.m.