plotting gene structure in R
2
0
Entering edit mode
6 months ago
chictu ▴ 10

I have a list of amplified genes with genomic intervals and would like to visualise the specific part of gene that is affected. Is there a way to plot a gene and highlight the region of interest? I was thinking a plot with introns and exons of the genes and genomic interval at the bottom with a box highlighting the region of interest but I have no idea how to go about it. Any guidance would be greatly appreciated.

Here is an example of an interval which contains the first few exons of EGFR: chr7:55234427-53437292 (hg19)

visualisation gene-structure • 796 views
ADD COMMENT
2
Entering edit mode
6 months ago
biofalconch ★ 1.1k

Maybe give ggtranscript a go

ADD COMMENT
2
Entering edit mode
6 months ago

ggtranscript seems to be a quite neat new package. The alternative, that is around a bit longer, is Gviz, but you can't straightforwardly combine it with ggplot. To highlight the region of interest, you could either use the Annotation track or the Highlight track.

ADD COMMENT
0
Entering edit mode

Thank you for the suggestions! I will give these a try.

ADD REPLY

Login before adding your answer.

Traffic: 1850 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