This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Packages in R for analyze premature stop codon and GC contents

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?

Thanks

Cai

r rna-seq
• 1,875 views • link
updated by Kristoffer Vitting-Seerup Kristoffer Vitting-Seerup • written by Kai_Qi Kai_Qi

1 answer

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.

• 0 views • link
written by Kristoffer Vitting-Seerup Kristoffer Vitting-Seerup

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

• 0 views • link
written by Kai_Qi Kai_Qi

Hello:

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

• 0 views • link
written by Kai_Qi Kai_Qi

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.

• 0 views • link
written by Kristoffer Vitting-Seerup Kristoffer Vitting-Seerup

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

• 0 views • link
written by Kai_Qi Kai_Qi

Log in to answer this question.

More posts like this