Question: Reproducing this picture
gravatar for F
12 days ago by
F3.3k wrote:


I have 21 .vcf files from tumour Vs normal samples. I likely annotated them individually by Ensembl Variant Effect Predictor (VEP) and snpEff. Now, I have 21 .txt files for each of these tools as results. For example this is results of one .vcf in Vep

I don't know how people produce this picture from their results

enter image description here

I found this post

How to create a mutation landscape (waterfall) plot with GenVisR

But, his input for tutorial is not like my Vep or snpEff

Any suggestion please?

sequencing snp R genome • 236 views
ADD COMMENTlink modified 12 days ago by Chris Miller20k • written 12 days ago by F3.3k
gravatar for Chris Miller
12 days ago by
Chris Miller20k
Washington University in St. Louis, MO
Chris Miller20k wrote:

GenVisR is a good R package for this, as is ProteinPaint from St Jude:

ADD COMMENTlink written 12 days ago by Chris Miller20k
gravatar for Friederike
12 days ago by
United States
Friederike2.5k wrote:

Did you check the help for waterfall or lolliplot? I haven't used those packages, but it seems to me you will need the genomic coordinates of both the gene as well as the mutations. You will likely have to do some data wrangling to parse the location of the mutation of the file you linked into three distinct columns that R can work with (note how the file from the tutorial has separate columns for chromosome, start and stop) -- all the details should hopefully be explained in the help of both functions, and if not, please be more specific in terms of which details you don't understand.

ADD COMMENTlink written 12 days ago by Friederike2.5k
gravatar for RamRS
12 days ago by
Houston, TX
RamRS20k wrote:

Like Friederike says, these are called lollipop plots indeed. You can search online for various tools that plot them, such as pbnjay's mutsneedle, cBioPortal's MutationMapper, etc, but each has its own limitations.

A search for mutation plot on Google images picks up this biostars post: How To Create Mutation Diagram In R Or In Any Tools?

ADD COMMENTlink modified 12 days ago • written 12 days ago by RamRS20k

Sorry, both Vep and snpEff results don't have reference allele column

I am seeing in most of your kindly suggested tools for visualization we need these columns

Chromosome  Start_Position  End_Position    Reference_Allele    Variant_Allele

Even in Vep results start and end positions have been merged.

I am not sure how to deal with incompatibility in input files.

ADD REPLYlink modified 11 days ago • written 11 days ago by F3.3k

The VCF file that you used for VEP should have that information

ADD REPLYlink written 11 days ago by Friederike2.5k

Sorry for this silly question;

You please imagine I have called somatic indels from cancer Vs normal samples and I have such .vcf files, what is the next step? I googled a lot but I am getting more confused. If my question is the mutation underlying this type of cancer, what these vcf files would say? I saw people use MutsigCV but I don't know why. For example Vep and snpEff do good job in feature selection then why people use MutsigCV.

ADD REPLYlink modified 11 days ago • written 11 days ago by F3.3k

VEP and snpEff annotate variants. MutSigCV, AFAIK, picks significantly mutated genes, a totally different task. Also, if you're using GDC's MAF format, they refer to Alt alleles differently, Tumor_Seq_Allele1 and Tumor_Seq_Allele2. Ref allele is still called Reference_Allele though, so that should not be the problem.

VCF files WILL have the ref allele column, there can be no VCF file that does not have REF.

EDIT: This offshoot is not related to the original post (which deals with reproducing a lollipop plot), please search the forum for discussions related to VCF annotation/MutSigCV.

ADD REPLYlink modified 11 days ago • written 11 days ago by RamRS20k

Thank you, because I thought why I am trying produce a lollipop plot, eventually what is my goal that is why I got concerned

ADD REPLYlink written 11 days ago by F3.3k

Sorry, could I use dNdScv instead of MutSigCV for finding significant variants?

ADD REPLYlink written 9 days ago by F3.3k

I'm sorry, I don't know. MutSigCV is as far as my exposure goes.

ADD REPLYlink written 8 days ago by RamRS20k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1787 users visited in the last hour