SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Why is it necessary to visualize the RNA-seq data in a genome browser? pinki999 Bioinformatics 4 08-05-2013 08:13 AM
visualize gmap results Kennels Bioinformatics 0 02-14-2012 01:49 PM
Visualize pair-end RNA-Seq in UCSC browser oscarluoinau Bioinformatics 0 11-21-2011 01:24 AM
Best way to visualize raw reads, and RPKMs? Bradburne Bioinformatics 6 08-22-2010 11:27 PM
How to visualize bowtie output? caspase Illumina/Solexa 2 11-16-2009 10:15 PM

Reply
 
Thread Tools
Old 05-06-2010, 10:47 AM   #1
dmotti
Junior Member
 
Location: Miami

Join Date: Mar 2010
Posts: 2
Default How to visualize RPKM in a browser?

Hallo everyone,

I'm new in the forum. Recently I've been dealing with some NGS data coming from the Illumina platform. I used the classical Tophat-Cufflinks pipeline to analyze the data.
But I'm dealing with the issue of visualizing RPKM values in a genome browser.

UCSC does not support GTF files, nor GFF, I can upload the wig files, but that would help only with the raw reads.

In another thread on this forum someone suggested a modification of the GTF file that would fit in Gbrowse. I tried to do that and uploaded the result on Gbrowse in MGI. No output. the track is just not visualized (even if no error output is given).

In this very moment I'm somehow stuck there, can anyone help me with that?

D.
dmotti is offline   Reply With Quote
Old 05-06-2010, 11:30 AM   #2
Siva
Member
 
Location: Ames, IA

Join Date: Nov 2008
Posts: 57
Default

Quote:
Originally Posted by dmotti View Post
Hallo everyone,

I'm new in the forum. Recently I've been dealing with some NGS data coming from the Illumina platform. I used the classical Tophat-Cufflinks pipeline to analyze the data.
But I'm dealing with the issue of visualizing RPKM values in a genome browser.

UCSC does not support GTF files, nor GFF, I can upload the wig files, but that would help only with the raw reads.

In another thread on this forum someone suggested a modification of the GTF file that would fit in Gbrowse. I tried to do that and uploaded the result on Gbrowse in MGI. No output. the track is just not visualized (even if no error output is given).

In this very moment I'm somehow stuck there, can anyone help me with that?

D.
Hi
I may be wrong but I think since GFF or GTF files are tab delimited files containing details of assembled features (transcripts or whatever), you may not be able to visualize it in a browser at least not like an aligned file. You could feed the GFF/GTF coordinates to a statistical software and output a graph/chart.

Siva
Siva is offline   Reply With Quote
Old 05-06-2010, 02:35 PM   #3
dmotti
Junior Member
 
Location: Miami

Join Date: Mar 2010
Posts: 2
Default

Hi Siva,

Thanks for your answer. Indeed you're right, GTF and GFF files represent discrete elements, they would only show as a coverage on the browser.
however, there is the possibility to unable an abundance representation of the elements by setting a score for each of them that should be between 1:1000.
Each element at that point would be represented in a different color. I had to force my RPKM values into a scale of 1 to 1000, and considering I have a much wider range it does not represent exactly what I'm looking for.


In this thread:

http://seqanswers.com/forums/showthr...ght=.expr+file

vkrishna suggested a modified version of the GTF file that could be read by Gbrowse. I tried that and it did not work, but the kind of representation I'm looking for is exactly the one suggested by the user.



The ENCODE project presents several tracks derived from RNA Seq data that can be visualized in Genome Browser, so there must be a way to submit RPKMs to the browser.

Any suggestion?

D.
dmotti is offline   Reply With Quote
Old 09-30-2010, 05:56 PM   #4
jdanderson
Member
 
Location: Davis, CA

Join Date: Sep 2010
Posts: 45
Default

Hello dmotti,

Did you ever find a solution to the visualization of cufflinks output files? I am currently in the same boat. I was thinking about using DESeq, but it seems that it only wants the read counts, not any other type of value.

Any guidance would be most appreciated.

Regards,
Johnathon
jdanderson is offline   Reply With Quote
Old 08-28-2015, 10:07 AM   #5
dawn1313
Member
 
Location: Tucson

Join Date: Aug 2015
Posts: 21
Default

Hi,

It seems I haven't found the exact response to what kind of genome browser for cufflinks outputs. Basically I want to visualize the reads abundance (hopefully also as a log scale) across the genome from either the Tophat output 'accepted_hits.bam' or the Cufflinks output 'transcripts.gtf'.

Is there any suggestion? Thank you very much!
Dawn
dawn1313 is offline   Reply With Quote
Reply

Tags
cufflinks, gbrowse, illumina, rpkm, ucsc

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 09:52 PM.


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