I Need To Plot 73 Mutations On A Large Gene With >60 Exons
5
5
Entering edit mode
11.9 years ago
Biomed 4.9k

I need to plot 73 mutations on a large gene for a publication. Th mutations of course must be correctly placed on the gene according to their genomic coordinates. Is there a tool, library , R module etc that I can use for that or do I need to draw it by hand in a drawing program. Some examples are http://www.nature.com/ng/journal/v43/n5/fig_tab/ng.810_F3.html and http://www.nature.com/ng/journal/v42/n11/fig_tab/ng.677_F1.html

Thank you

mutation gene plot • 6.0k views
ADD COMMENT
4
Entering edit mode
11.9 years ago

You can also try Fancy Gene.

ADD COMMENT
3
Entering edit mode
11.9 years ago
Bill Pearson ★ 1.0k

I strongly suggest you learn a little bit of SVG to do this. SVG images are easily edited with Adobe Illustrator, and can be visualized in many web browsers.

Basically, you need to set up a coordinate system for your genes/exons and mutations. Learn how to draw some boxes and lines with SVG and you will quickly see how easy it can be.

ADD COMMENT
2
Entering edit mode
11.9 years ago

if your protein is available as a genbank file, save it as XML and process it with XSLT and the following stylesheet: https://github.com/lindenb/xslt-sandbox/blob/master/stylesheets/bio/ncbi/gb2svg.xsl

xsltproc gb2svg.xsl protein.xml > protein.svg

it should prepare a SVG file that you can later modify with inkscape

ADD COMMENT
0
Entering edit mode

I couldn't get it to work. Nothing happens when I run the above command with the genbank xml of the human protein.

ADD REPLY
1
Entering edit mode
11.9 years ago

You can see an example of how this is done in Alexa-seq using R code. The figure shows LAMA3 with 77 exons. There are two representations of the gene, one that shows the overall structure of exons and their proportional distance and the other has a sort of zoomed in view with exons depicted on log2 scale and introns on log10 scale. In this case some alternative splicing information is marked up. But instead you could indicate mutations. Unfortunately this is not part of an existing R package and is just done with basic R plotting functions in an R script of about 300 lines. But, maybe it will give you some ideas.

Alexa-seq example

ADD COMMENT
0
Entering edit mode
11.9 years ago
Bjoernsen ▴ 40

If it is an eukaryotic protein and if there is a genome available (at diArk), you can try WebScipio.

ADD COMMENT

Login before adding your answer.

Traffic: 2369 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6