SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
can small rnaseq data be analyzed like rnaseq data? PFS Bioinformatics 5 05-02-2017 09:16 AM
Basic steps in prokaryotic RNAseq data analysis Muhidini Bioinformatics 0 05-09-2012 07:27 AM
[NGS - analysis of gene expression data] Machine Learning + RNAseq data Chuckytah Bioinformatics 7 03-05-2012 04:16 AM
Companies offering NGS data analysis including RNAseq agseq Bioinformatics 0 05-05-2011 07:24 AM
ERANGE 3.1 RNAseq analysis Xi Wang Bioinformatics 2 02-03-2010 07:03 PM

Reply
 
Thread Tools
Old 05-24-2012, 08:45 AM   #1
chknbio
Member
 
Location: Baltimore

Join Date: May 2012
Posts: 14
Default RNAseq data analysis

I recently started some RNAsequencing of tumor samples. My main question is regarding the differential gene expression among samples and controls. I have successful utilized tophat, cufflinks, and cuffdiff. I am new to bioinformatics and RNAseq analysis. My experience utilizing Linux terminals for analysis is only recent. I have seen a few tutorials online for RNAseq, but some of them only take you through the cufflinks data retrieval. I am looking for what next....the steps and analysis after cuffdiff. Are there programs that exist or tutorials that exist for running a script in python or R for analysis of the data? How do others go about analyzing RNAseq data for differential expression among samples?

Any advice on the downstream analysis of RNAseq data after the tuxedo suite would be greatly appreciated.

Thanks!
chknbio is offline   Reply With Quote
Old 05-24-2012, 09:29 AM   #2
ETHANol
Senior Member
 
Location: Western Australia

Join Date: Feb 2010
Posts: 308
Default

cumbRbund is the next step in the tuxedo suite. After that it all depends on what you want to do.
__________________
--------------
Ethan
ETHANol is offline   Reply With Quote
Old 05-24-2012, 09:30 AM   #3
usad
Member
 
Location: aachen

Join Date: Sep 2009
Posts: 53
Default

Hi

whilst you get much further with the command line, you can always try out our tool RobiNA http://mapman.gabipd.org/web/guest/robin which will give you a graphical user interface and run all necessary steps in the background to call differential expression. It also comes with a long manual.

Apart from that you might check out the tutorials in the seqanswers for running your data through.

Cheers,
björn
usad is offline   Reply With Quote
Old 05-24-2012, 09:37 AM   #4
chknbio
Member
 
Location: Baltimore

Join Date: May 2012
Posts: 14
Default

Quote:
Originally Posted by ETHANol View Post
cumbRbund is the next step in the tuxedo suite. After that it all depends on what you want to do.
yes, I have used cummeRbund to plot some of the graphs. Which look nice, but I am more curious in pulling out differential expressed genes. I would like to look at the data globally vs. picking and choosing genes. From my use of cummeRbund, I was able to obtain plots which show differences and similarities between samples however, I wasn't able to see where those differences were located (i.e., specific genes or novel transcripts) - just dots on a graph. Maybe I am missing the feature that allows one to export the transcripts or genes that are upregulated or downregulated x-fold from the control or other samples.
chknbio is offline   Reply With Quote
Old 05-24-2012, 09:50 AM   #5
chknbio
Member
 
Location: Baltimore

Join Date: May 2012
Posts: 14
Default

Quote:
Originally Posted by usad View Post
Hi

whilst you get much further with the command line, you can always try out our tool RobiNA http://mapman.gabipd.org/web/guest/robin which will give you a graphical user interface and run all necessary steps in the background to call differential expression. It also comes with a long manual.

Apart from that you might check out the tutorials in the seqanswers for running your data through.

Cheers,
björn
Hi Bjorn,

Thank you for the reply. And for the link on this tool. I already have the cuffdiff output and I would use command line but I am wondering where to begin. How do I retrieve all of the differential regulated genes or transcripts?

-Tom
chknbio is offline   Reply With Quote
Old 05-24-2012, 10:56 AM   #6
chknbio
Member
 
Location: Baltimore

Join Date: May 2012
Posts: 14
Default

This is a very good tutorial here, but it still leaves some questions unanswered.
http://vallandingham.me/RNA_seq_diff...xpression.html

How do you determine what cuffdiff lists as noteworthy genes for upregulation or downregulated? I am only seeing the FPKM for each gene and the lo and high confidence. How can I see the genes that are predicted to be noteworthy as mentioned in this tutorial?

This author took the log2 of all FPKM values? Isn't the FPKM a value of within samples expression compared to overall transcript numbers?

I would also like to know what scripts or programs did he run in R?

Any guidance or other tutorials that will lead me forward after obtaining cuffdiff file would be helpful. I have four samples and run cuffdiff with tophat bam files for each and the chicken annotation gtf.

Thanks!
chknbio is offline   Reply With Quote
Old 05-24-2012, 11:40 AM   #7
ETHANol
Senior Member
 
Location: Western Australia

Join Date: Feb 2010
Posts: 308
Default

There should be a file gene_exp.diff or something like that which has the results of the statistical test. Sort by qvalue and chose a cutoff you like. Is that what you are asking?

You can then take your significantly expressed gene list and feed that into goseq if you want to do gene ontology analysis.
__________________
--------------
Ethan
ETHANol is offline   Reply With Quote
Old 05-24-2012, 12:20 PM   #8
chknbio
Member
 
Location: Baltimore

Join Date: May 2012
Posts: 14
Default

Quote:
Originally Posted by ETHANol View Post
There should be a file gene_exp.diff or something like that which has the results of the statistical test. Sort by qvalue and chose a cutoff you like. Is that what you are asking?

You can then take your significantly expressed gene list and feed that into goseq if you want to do gene ontology analysis.
This is helpful. Thank you!
chknbio is offline   Reply With Quote
Old 05-25-2012, 12:44 AM   #9
EGrassi
Member
 
Location: Turin, Italy

Join Date: Oct 2010
Posts: 66
Default

The tuxedo procotol article has a section on what you want: http://www.nature.com/nprot/journal/...tml#/procedure check steps 16-18, there is some simple R code that you could follow and expand.
EGrassi is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:41 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO