Question: Packages in R for analyze premature stop codon and GC contents
gravatar for Kai_Qi
14 months ago by
Chicago, IL
Kai_Qi100 wrote:

Hi Everyone:

I am trying to go through analysis of alternative splicing from RNA-seq. So far I have got the coordinates of some intron, and I want to analyzed the feature of these intron: GC contents, whether and where it will generate premature stop codon, strength of 5' or 3' splice site.

Is there good package to use for these analysis?



rna-seq R • 306 views
ADD COMMENTlink modified 14 months ago by kristoffer.vittingseerup3.5k • written 14 months ago by Kai_Qi100
gravatar for kristoffer.vittingseerup
14 months ago by
European Union
kristoffer.vittingseerup3.5k wrote:

If you want to analyse the consequence of alternative splicing you need transcript level quantification of events (since the consequence of an intron retention depends on all of the connected up and down-stream exons). You can read more about different approaches to analysis of splicing here. Also if you want transcript level quantifications you can read more about such considerations here. My R package IsoformSwitchAnalyzeR will enable you to find such cases quite easily - both on a genome wide level and of individual genes - you can find examples of the analysis produced by the package here.

In terms of analysing biological sequences the Bioconductor package Biostrings is a very nice set of tools for doing so. You can find a nice example of how people use it to calculate GC content of CpG islands here.

ADD COMMENTlink written 14 months ago by kristoffer.vittingseerup3.5k

Thank you very much for your advice. I am going to give it a try

ADD REPLYlink written 14 months ago by Kai_Qi100


I read about the manual, it is comprehensive but intimidating for me. I want to give it a try on the data I have. So far I have run RSEM to do the quantification of the isoform. I have a file called "rsem.isoforms.results", but I have no idea of how to use it. Do you have a step by step manual. The manual from the link import data from the package which I don't know how to use in my case. Thanks

ADD REPLYlink written 5 months ago by Kai_Qi100

IsoformSwitchAnalyzeR can work directly from the "rsem.isoforms.results". I suggest you take a look at the quick start part of the vignette which should describe exactly what you need including code you can use.

ADD REPLYlink written 5 months ago by kristoffer.vittingseerup3.5k

OK. I will read it again to make it through.

ADD REPLYlink written 5 months ago by Kai_Qi100
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: 2543 users visited in the last hour