SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
searching tool for ChIP-seq analysis of histone modification analysis tujchl Epigenetics 5 03-13-2013 09:00 AM
illumina smallRNA adapter sequence for downstram analysis + miRNA analysis steps ndeshpan Bioinformatics 2 06-14-2011 09:44 PM
how to mapp Next generation sequencing (NGS)data on Refrence genome ???? saurabh_igib Illumina/Solexa 3 04-16-2010 05:38 AM
ABRF2010-Next Gen Sequencing Analysis:Platform Independent Analysis You Can Use MichelleMD Events / Conferences 0 12-04-2009 06:35 PM

Reply
 
Thread Tools
Old 06-09-2011, 12:38 AM   #1
abpatel2288
Junior Member
 
Location: Anand, INDIA

Join Date: Jun 2011
Posts: 6
Exclamation Refrence for transcriptiome analysis

Should i use DNA ref or RNA ref. for mapping using gs reference mapper ??

What data base set should be used for finding expression RPKM, from transcriptome reads of cow?

Sequencing was carried by 454 gs-flx.
abpatel2288 is offline   Reply With Quote
Old 06-09-2011, 09:40 PM   #2
jmartin127
Member
 
Location: San Francisco

Join Date: Mar 2010
Posts: 15
Default

Since the reads you get from 454 will be DNA (the cDNA is sequenced), then you should align to a DNA reference. With 454 you may not have enough sequencing depth to get good RPKM measurements from a cow. How many reads do you have?

Quote:
Originally Posted by abpatel2288 View Post
Should i use DNA ref or RNA ref. for mapping using gs reference mapper ??

What data base set should be used for finding expression RPKM, from transcriptome reads of cow?

Sequencing was carried by 454 gs-flx.
jmartin127 is offline   Reply With Quote
Old 06-09-2011, 11:29 PM   #3
abpatel2288
Junior Member
 
Location: Anand, INDIA

Join Date: Jun 2011
Posts: 6
Default

Thanks dear
jmartin for your reply, i having around 529087 reads with minimum length of 40 after removing MID and other adapters..

as you told RPKM values are average. i have used CLC gemone work bench and NextGene software's for RPKM.

Results of both softwares are having much difference in RPKM values of same reads, i am totally puzzled why i am getting so much difference in value, and on which i should rely for further analysis .
abpatel2288 is offline   Reply With Quote
Old 06-10-2011, 06:38 AM   #4
jmartin127
Member
 
Location: San Francisco

Join Date: Mar 2010
Posts: 15
Default

I see what you mean, I cannot say for sure why the values are so different from the two software programs. Are you aligning to the genome or to the transcripts themselves to calculate RPKM? Usually I work with Illumina reads, so the tools I typically use wouldn't be applicable to you unfortunately.

Quote:
Originally Posted by abpatel2288 View Post
i have used CLC gemone work bench and NextGene software's for RPKM. Results of both softwares are having much difference in RPKM values of same reads, i am totally puzzled why i am getting so much difference in value, and on which i should rely for further analysis .
jmartin127 is offline   Reply With Quote
Old 06-10-2011, 07:43 PM   #5
abpatel2288
Junior Member
 
Location: Anand, INDIA

Join Date: Jun 2011
Posts: 6
Default

I am aligning reads to genome... Is there any other option which i would check ?
abpatel2288 is offline   Reply With Quote
Old 06-13-2011, 06:36 AM   #6
jmartin127
Member
 
Location: San Francisco

Join Date: Mar 2010
Posts: 15
Default

Often when I've done this type of analysis, I've aligned my reads to the spliced transcript sequences. For most organisms this will be listed as the "XXX_cDNA.fasta", for example. To get the expression values you can basically count the number of reads aligning to each transcript, and then divide by the gene length and the number of reads in your sample to get RPKM (you would also have to multiply by 1000*1000000, since it is Reads Per Kilobase of exon model per Million mapped reads).

Quote:
Originally Posted by abpatel2288 View Post
I am aligning reads to genome... Is there any other option which i would check ?
jmartin127 is offline   Reply With Quote
Old 06-13-2011, 06:41 AM   #7
vineeth_s
Junior Member
 
Location: Germany

Join Date: Jun 2011
Posts: 9
Default Another way

If you have done RNA-Seq (which is what you have I would guess, considering you say transcriptome analysis), perhaps you could give Trinity, from the Broad Institute, a run

Trinity is a de novo transcriptome assembler, voiding the need for you to have a reference

And here's the paper

If you do use it, it would be great to hear what your experiences with it were
vineeth_s is offline   Reply With Quote
Old 06-13-2011, 07:44 PM   #8
abpatel2288
Junior Member
 
Location: Anand, INDIA

Join Date: Jun 2011
Posts: 6
Default

will it work with 454, GS FLX
abpatel2288 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 02:25 AM.


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