SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Viral genome assembly from RNA-seq data mcastro RNA Sequencing 2 11-15-2015 05:21 AM
extracting seq from viral genome honey Bioinformatics 3 12-20-2012 03:44 AM
viral genome alignment mathew Bioinformatics 8 12-03-2012 08:33 AM
annotation for viral genome honey Bioinformatics 3 10-07-2012 09:18 PM
how to optimize viral genome assembly in WGS 7.0? rexxi Bioinformatics 1 08-01-2012 11:01 PM

Reply
 
Thread Tools
Old 06-11-2013, 01:17 AM   #1
vijesh
Junior Member
 
Location: kerala

Join Date: Jun 2013
Posts: 7
Question Viral genome vs transcripts

i have a fasta file containing all the virus genome which will affect cassava. i have transcript file also in fasta. in transcript file each transcript have 100 nucleotide length. I want to know whether these transcripts are inside the virus genom. any one help me pls?
vijesh is offline   Reply With Quote
Old 06-11-2013, 03:14 AM   #2
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

One relatively straight-forward way to attack this problem would be to make a combined fasta file for the viral genomes and cassava and then map the transcripts/transcript fragments against it. You can then get an idea how likely each fragment is to come from the host or virus. I assume that tophat or something like that would be an appropriate aligner, since presumably you have a mixture of host (i.e. spliced) and viral (I assume transcripts would be single exon, but I've only made viruses, not studied them) reads and not doing so might bias things.

I should note that I've never done what you're trying to do, but if no one replies with a better idea then this is enough to get you started.
dpryan is offline   Reply With Quote
Old 06-11-2013, 08:51 PM   #3
vijesh
Junior Member
 
Location: kerala

Join Date: Jun 2013
Posts: 7
Default doubt regarding tophat

Quote:
Originally Posted by dpryan View Post
One relatively straight-forward way to attack this problem would be to make a combined fasta file for the viral genomes and cassava and then map the transcripts/transcript fragments against it. You can then get an idea how likely each fragment is to come from the host or virus. I assume that tophat or something like that would be an appropriate aligner, since presumably you have a mixture of host (i.e. spliced) and viral (I assume transcripts would be single exon, but I've only made viruses, not studied them) reads and not doing so might bias things.

I should note that I've never done what you're trying to do, but if no one replies with a better idea then this is enough to get you started.
can we put the two files as input in tophat.means that genome file in fasta and also the transcript file?
vijesh is offline   Reply With Quote
Old 06-12-2013, 12:26 AM   #4
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

Quote:
Originally Posted by vijesh View Post
can we put the two files as input in tophat.means that genome file in fasta and also the transcript file?
Sort of, I can read a few different meanings into what you wrote so I'll write some pseudo-commands just to ensure that my meaning is clear. Let us assume that your Cassava genome file is called "Cassava.fa" with an annotation file "Cassava.gtf", your virus genome file is called "Viruses.fa", and the transcripts whose origin you're interested in determining are in a file called "reads.fa". The work flow would then be something like the following:

Code:
cat Cassava.fa Viruses.fa > CombinedGenome.fa
bowtie2-build CombinedGenome.fa Combined
tophat -G Cassava.gtf Combined reads.fa
You can then see which reads originated from which chromosome/contig or whichever virus by using samtools view on the resulting accepted_hits.bam file (you'll probably want to filter out reads that map equally to both Cassava and viral genomes, for which you'll probably have to write a program).

It would be good to compare the results with and without using the GTF annotation, off-hand I'm not entirely sure how or if that might bias things.
dpryan is offline   Reply With Quote
Old 06-14-2013, 12:11 AM   #5
vijesh
Junior Member
 
Location: kerala

Join Date: Jun 2013
Posts: 7
Red face viral genom vs transcript

i dont have the annotation file(gtf) .i have only the virus genome file (fasta) and transcript file(fasta).....
vijesh is offline   Reply With Quote
Old 06-14-2013, 12:20 AM   #6
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

You don't actually have to use the annotation file, I expect it might have just made things a bit more reliable (or maybe not, one would have to check).
dpryan 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:03 AM.


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