Knowledge base

Sushi package draws gene structure diagram

Introduction

The plotGenes of the Sushi package can plot gene structures, such as the visualization of several transcriptome structures.

Code explanation

Input
  • bed format data frame
head(Sushi_genes.bed)

 The structure is as follows

chrom    start     stop  gene score strand1 
chr15 73017309 73017438 COX5A     .     - 
chr15 72999672 72999836 COX5A     .     - 
chr15 73003042 73003164 COX5A     .     - 
chr15 73006160 73006281 COX5A     .     - 
chr15 73008510 73008626 COX5A     .     -
Sushi_transcripts.bed[1:20,]

  The structure is as follows

  chrom start stop gene score strand type
 chr15 73062668 73062770 ENST00000362710 0.000000 -1 exon
 chr15 73097788 73097929 ENST00000361900 0.000000 1 exon
 chr15 73097264 73097365 ENST00000361900 0.000000 1 exon
 chr15 73095987 73096143 ENST00000361900 0.000000 1 exon
 chr15 73092071 73092199 ENST00000361900 0.000000 1 exon
 chr15 73091234 73091240 ENST00000361900 0.000000 1 exon
 chr15 73017309 73017408 ENST00000322347 31.488695 -1 exon
 chr15 73006160 73006281 ENST00000322347 31.488695 -1 exon 
Code example
library('Sushi')  # Load Sushi包

Sushi_data = data(package ='Sushi') # List Sushi own data
data(list = Sushi_data$results[,3]) # Load data

chrom = "chr15"
chromstart  = 72998000
chromend = 73020000

 pg = plotGenes(Sushi_genes.bed,chrom,chromstart,chromend,
   types=Sushi_genes.bed$type,maxrows=1,bheight=0.2,
   plotgenetype="arrow",bentline=FALSE,labeloffset=.4,
   fontsize=1.2,arrowlength = 0.025,labeltext=TRUE)
   # Because our typse input is Sushi_genes.bed$type, which is actually NULL, the default "exon" will be used, and exon and utr are optional

 labelgenome( chrom, chromstart,chromend,n=3,scale="Mb")  #Add genome coordinates

plotgenetype="box" can display the area whose type is exon and utr in the form of box

library('Sushi')  # Load Sushi bag

Sushi_data = data(package ='Sushi') # List Sushi's own data
data(list = Sushi_data$results[,3]) # Load data

chrom = "chr15"
chromstart = 72965000
chromend = 72990000

pg = plotGenes(Sushi_transcripts.bed,chrom,chromstart,chromend ,
types = Sushi_transcripts.bed$type,
colorby=log10(Sushi_transcripts.bed$score+0.001),
colorbycol= SushiColors(5),colorbyrange=c(0,1.0),
labeltext=TRUE,maxrows=50,height=0.4,plotgenetype="box")
# colorby is the column that specifies the color change, and colorbycol is the color palette that specifies the color change

labelgenome( chrom, chromstart,chromend,n=3,scale="Mb") #Add genome coordinates

addlegend(pg[[1]],palette=pg[[2]],title="log10(FPKM)",side="right",
bottominset=0.4,topinset=0,xoffset=-.035,labelside="left",
width=0.025,title.offset=0.055)  #Add legend

Reference

https://www.bioconductor.org/packages/release/bioc/vignettes/Sushi/inst/doc/Sushi.pdf

62 thoughts on “Sushi package draws gene structure diagram

  1. You have made some good points there. I checked on the internet
    for more info about the issue and found most individuals
    will go along with your views on this website.

  2. I’ve been exploring for a little bit for any high-quality articles or
    weblog posts on this kind of space . Exploring in Yahoo I finally stumbled upon this web site.
    Reading this info So i am happy to express that I’ve a very just right uncanny feeling I came upon just
    what I needed. I such a lot for sure will make sure to don?t put out of your mind this website and give it a
    look regularly.

  3. Pretty component of content. I simply stumbled upon your blog and in accession capital to assert that I
    get in fact enjoyed account your blog posts. Any way I’ll be subscribing in your
    feeds and even I fulfillment you get right of entry to consistently
    fast.

  4. May I just say what a comfort to find somebody who truly understands what they are talking about online.
    You certainly know how to bring a problem to light and make it important.
    A lot more people need to read this and understand this side of
    your story. I was surprised you are not more popular because you most certainly have the gift.

  5. Excellent post. I was checking continuously this blog and I am impressed!
    Very useful info specifically the last part 🙂 I care for such information much.
    I was looking for this certain information for a very
    long time. Thank you and best of luck.

  6. Howdy! Do you know if they make any plugins to help with Search Engine
    Optimization? I’m trying to get my blog to rank for some targeted keywords but I’m not seeing
    very good success. If you know of any please share. Many thanks!

  7. Thank you, I’ve recently been looking for information approximately this topic for a long time
    and yours is the greatest I’ve found out so far. However, what
    about the conclusion? Are you sure concerning the supply?

  8. My programmer is trying to persuade me to move to .net from PHP.
    I have always disliked the idea because of the costs.
    But he’s tryiong none the less. I’ve been using Movable-type on a variety of websites for about a year
    and am worried about switching to another platform.
    I have heard great things about blogengine.net.
    Is there a way I can transfer all my wordpress content into it?
    Any help would be really appreciated!

  9. Hey there I am so delighted I found your weblog, I really found you
    by accident, while I was researching on Yahoo for something else,
    Regardless I am here now and would just like to say thanks a lot for a remarkable post and a all round exciting blog (I also love the theme/design), I don’t have time to read through it all at the moment but I
    have saved it and also included your RSS feeds, so
    when I have time I will be back to read much more, Please do keep up the awesome jo.

  10. Thank you for any other informative website. Where else may just I am getting that kind of info written in such a perfect way?
    I have a mission that I’m simply now working on, and I’ve been on the look out for
    such information.

  11. Heya i am for the primary time here. I found this board and I to find It truly helpful & it helped
    me out a lot. I hope to provide one thing again and help others such as you
    helped me.

  12. It is appropriate time to make some plans for the future and it’s time to be happy.
    I have read this publish and if I may just I desire to recommend you
    few attention-grabbing things or advice. Perhaps you can write next articles relating to this article.

    I want to learn even more issues about it!

  13. Pretty section of content. I just stumbled upon your web site and in accession capital to assert that I acquire in fact enjoyed account your blog
    posts. Any way I’ll be subscribing to your feeds and even I achievement you access consistently fast.

  14. It is the best time to make a few plans for
    the future and it’s time to be happy. I’ve read this put up and if I may I wish to
    suggest you some attention-grabbing things or advice.
    Maybe you could write subsequent articles relating to this article.

    I desire to learn even more things about it!

  15. Hi! This is kind of off topic but I need some help from an established blog.

    Is it hard to set up your own blog? I’m not very techincal but I can figure things out pretty quick.
    I’m thinking about creating my own but I’m not sure where to begin. Do you have any ideas or suggestions?
    Thank you

  16. Attractive part of content. I just stumbled upon your website and in accession capital to assert
    that I acquire in fact loved account your blog posts.
    Any way I’ll be subscribing to your augment or even I fulfillment you get admission to persistently fast.

  17. Hey there, I think your site might be having browser compatibility issues.
    When I look at your blog site in Opera, it looks fine but when opening
    in Internet Explorer, it has some overlapping. I just wanted to give you a
    quick heads up! Other then that, amazing blog!

  18. My partner and I absolutely love your blog and find the majority
    of your post’s to be just what I’m looking for.
    can you offer guest writers to write content for yourself?
    I wouldn’t mind creating a post or elaborating on most of the subjects you write regarding here.
    Again, awesome web site!

  19. Hey there I am so glad I found your blog page, I really found you by
    error, while I was browsing on Google for something else,
    Nonetheless I am here now and would just like to say cheers for a fantastic post and a all round entertaining blog
    (I also love the theme/design), I don’t have time
    to read through it all at the moment but I have bookmarked it and also included your RSS feeds, so when I have time I will
    be back to read much more, Please do keep up the great job.

  20. I blog quite often and I genuinely appreciate your information. This article has truly peaked my interest.

    I am going to bookmark your site and keep checking for new information about once per week.
    I opted in for your Feed too.

  21. hello!,I really like your writing very a lot! share we communicate extra about your post
    on AOL? I require an expert in this house to solve my problem.
    Maybe that’s you! Taking a look ahead to see you.

  22. Hi there, You have done an ihcrеdible job.
    I’ll certainly digg it and personally suggest tο mmy friends.

    I’m confiden thy will be benefited from this sitе.

  23. Hey I know this is off topic but I was wondering if you knew of any widgets I could add to my blog that automatically
    tweet my newest twitter updates. I’ve been looking for a plug-in like this for
    quite some time and was hoping maybe you would have some experience with something like this.
    Please let me know if you run into anything.
    I truly enjoy reading your blog and I look forward to your new updates.

  24. Please let me know if you’re looking for a writer for your site.
    You have some really good posts and I think I would be
    a good asset. If you ever want to take some of the load off, I’d love to write some
    articles for your blog in exchange for a link back to
    mine. Please send me an email if interested. Regards!

  25. I got this website from my buddy who shared with me concerning this
    web site and now this time I am visiting this site and reading very
    informative articles or reviews at this place.

  26. Normally I do not learn article on blogs, however I would like
    to say that this write-up very pressured me to try and do so!

    Your writing taste has been surprised me. Thanks, quite nice article.

Leave a Reply

Your email address will not be published. Required fields are marked *