This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to visualize gene structure of a (prokaryotic) locus?

Hello, guys!

I was wondering, how an image like this can be automatically (programmatically) generated from the initial data in a text file?

enter image description here

The figure is taken from https://www.ncbi.nlm.nih.gov/pubmed/29168502

Could anyone recommend a good R package for this?

Thanks in advance, Ivan

r genome gene visualisation python
• 4,434 views • link
updated by h.mon h.mon • written by ivan.antonov ivan.antonov

You can try EasyFig. It does pretty nice figures.

That looks similar to the output of MultiGeneBlast, but I’d wager it was probably made by hand.

• 1 views • link
written by Joe Joe

Thanks for the response! It may be difficult to accurately visualize the gene lengths and intergenenic distances by hand, though..

• 1 views • link
written by ivan.antonov ivan.antonov

what would be your input ? a GFF ? a genbank file ?

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Yeah, something like GFF format. I can prepare a list of genomic coordinates for the genes of interest in any format...

• 1 views • link
written by ivan.antonov ivan.antonov

3 answers

I just found the gggenes R package, it seems it was designed with your question in mind. See a small example in my answer to Genome Diagram in R.

• 1 views • link
written by h.mon h.mon

This is just perfect -- thanks a lot!

• 1 views • link
written by ivan.antonov ivan.antonov

You might also be interested in the forum I threw together a while back.

Many of the tools in that thread will do what you need.

Interesting/useful/beautiful bio sequence viz tools

• 1 views • link
written by Joe Joe

Thanks -- DNAFeaturesViewer and GenomeDiagrams python libs look like what I need:)

• 1 views • link
written by ivan.antonov ivan.antonov

Use the XML genbank output of NCBI Eutilities for your genome and transform it using a XSLT stylesheet to SVG. Here is an example using two accessions in the range 1000-50000.

you can then edit it with inkscape.

• 1 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

Log in to answer this question.

More posts like this