maftools - An R package to summarize, analyze and visualize MAF files.
With advances in Cancer Genomics, Mutation Annotation Format (MAF) is being widley accepted and used to store variants detected. The Cancer Genome Atlas Project has seqenced over 30 different cancers with sample size of each cancer type being over 200. The resulting data consisting of genetic variants is stored in the form of Mutation Annotation Format. This package attempts to summarize, analyze, annotate and visualize MAF files in an efficient manner either from TCGA sources or any in-house studies as long as the data is in MAF format. Maftools can also handle ICGC Simple Somatic Mutation format.
maftools is on
Please cite the below if you find this tool useful for you.
Mayakonda, A. and H.P. Koeffler, Maftools: Efficient analysis, visualization and summarization of MAF files from large-scale cohort based cancer studies. bioRxiv, 2016. doi: http://dx.doi.org/10.1101/052662
MAF field requirements.
MAF files contains many fields ranging from chromosome names to cosmic annotations. However, most of the analysis in maftools uses following fields. Please stick to MAF specifications for better results.
- Mandatoty fields: Hugo_Symbol, Chromosome, Start_Position, End_position, Reference_Allele, Tumor_Seq_Allele2, Variant_Classification, Variant_Type and Tumor_Sample_Barcode.
- Recommended optional fields: non MAF specific fields containing vaf and amino acid change information. Complete specififcation of MAF files can be found on NCI TCGA page.
NOTE: If you have variants stored as VCFs or as an MAF like tab seperated format, convert them to MAF using vcf2maf/maf2maf. Merge MAFs from all samples into a single MAF before processing with maftools.
Vignette and a case study.
A complete documentation of maftools using TCGA LAML1 as a case study can be found here.
Simple usage: Just read maf file using
read.maf and pass the resulting maf object to any one of the function for plotting or analysis or set operations.
Stuffs maftools can do.
- Detect cancer driver genes based on positional clustering of variants2.
- Detect Mutually exclusive set of genes3.
- Compare two MAF files (cohorts) to detect differentially mutated genes.
- Add pfam domains and summarize.
- Extract mutational signatures and compare them to validated signatures4.
- APOBEC Enrichment score estimation5.
- Tumor heterogenity and MATH (Mutant-Allele Tumor Heterogeneity) score estimation6.
- Read and summarize GISTIC results.
- Pan-cancer analysis/comparisison
- Survival analysis
- Compare mutation load against all 33 TCGA cohorts
- Rich Visualizations
- Make oncoplots7.
- Make lollipop plots.
- Map variants on copy number (CBS) segments
- Forest plots
- Plot Transitions and Transversions.
- Plot maf summary.
- Rainfall plots and change point detection
- Annotate variants locally using Oncotator API.
- Convert Annovar annotations into MAF.
- Convert ICGC simple somatic mutation format into MAF.
Easy way: Install from Bioconductor.
if (!require("BiocManager")) install.packages("BiocManager") BiocManager::install("maftools")
Install from Github for updated features (some of functions from here may not be available on Bioconductor release branch).
#Install Bioconductor dependencies. if (!requireNamespace("BiocManager", quietly=TRUE)) install.packages("BiocManager") BiocManager::install("ComplexHeatmap") BiocManager::install("VariantAnnotation") BiocManager::install("Biostrings") #Install maftools from github repository. library("devtools") install_github(repo = "PoisonAlien/maftools")
For full documentation please refer to vignette.
- Cancer Genome Atlas Research, N., Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia. N Engl J Med, 2013. 368(22): p. 2059-74.
- Tamborero, D., A. Gonzalez-Perez, and N. Lopez-Bigas, OncodriveCLUST: exploiting the positional clustering of somatic mutations to identify cancer genes. Bioinformatics, 2013. 29(18): p. 2238-44.
- Leiserson, M.D., Wu, H.T., Vandin, F. & Raphael, B.J. CoMEt: a statistical approach to identify combinations of mutually exclusive alterations in cancer. Genome Biol 16, 160 (2015).
- Alexandrov, L.B., et al., Signatures of mutational processes in human cancer. Nature, 2013. 500(7463): p. 415-21.
- Roberts SA, Lawrence MS, Klimczak LJ, et al. An APOBEC Cytidine Deaminase Mutagenesis Pattern is Widespread in Human Cancers. Nature genetics. 2013;45(9):970-976. doi:10.1038/ng.2702.
- Mroz, E.A. & Rocco, J.W. MATH, a novel measure of intratumor genetic heterogeneity, is high in poor-outcome classes of head and neck squamous cell carcinoma. Oral Oncol 49, 211-5 (2013).
- Gu, Z., Eils, R. & Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics (2016).