SEQanswers

Go Back   SEQanswers > Applications Forums > De novo discovery



Similar Threads
Thread Thread Starter Forum Replies Last Post
Trinity Assembly, Isoforms vs. Paralogs? BGould Bioinformatics 1 01-15-2014 12:51 AM
de novo assembly using Trinity versus Velvet-Oases Nol De novo discovery 8 10-26-2013 11:56 AM
get wrong when I do de novo assembly by trinity for two 22.6G reads file taoxiang180 Bioinformatics 4 08-20-2012 09:52 PM
Transcriptome assemblers - Oases, Trinity samanta General 0 10-04-2011 10:48 AM
Trinity - de novo transcriptome assembler samanta De novo discovery 2 07-29-2011 12:19 PM

Reply
 
Thread Tools
Old 10-08-2012, 10:24 PM   #1
upendra_35
Senior Member
 
Location: USA

Join Date: Apr 2010
Posts: 102
Default Trinity transcriptome assembly

After Trinity assembler finished its assembly i managed to calculate the basic statistics of the assembly which are as below

File Number Total Size Min Size Max Size Average Size Median Size N50 Trinity.fasta 158863 176660784 201 22887 1112.03 665 1863

Size@1Mbp Number@1Mbp Size@2Mbp Number@2Mbp Size@4Mbp Number@4Mbp Size@10Mbp Number@10Mbp 11440 65 8461 170 7088 430 5424 1417

Now my question is does these values look reasonable? Though N50 looks good i am worried about the number of transcripts that are less than 1kb (~ 60%) of the overall transcripts. Is this normal in Trinity?

Also how do people normally do downstream analysis after getting the assembly to select the best transcritps. I ask this because the number of Transcripts is way higher than expect number of genes in related species.

Thanks........
upendra_35 is offline   Reply With Quote
Old 10-15-2012, 02:02 AM   #2
mbayer
Member
 
Location: Dundee, Scotland

Join Date: Mar 2009
Posts: 29
Default

Hi,

personally I think that looks reasonable, assuming you have a eukaryotic organism -- the average gene length in eukaryotes is supposed to be in the 1,500 bp region. What organism is this, and do you know how common alternative splicing is in that species? That would obviously affect your number of transcripts relative to the number of genes. Also, given that most alternative splicing produces transcripts that are shorter than the full length mRNA, an average transcript length of 1,112 seems reasonable.

To evaluate your assembly I would run something on the transcripts that predicts proteins, like getorf from the EMBOSS suite of tools, then select the longest predicted protein and BLAST this against related species. This will give you an idea of how good your assembly is.

cheers

Micha
mbayer is offline   Reply With Quote
Old 10-15-2012, 12:00 PM   #3
upendra_35
Senior Member
 
Location: USA

Join Date: Apr 2010
Posts: 102
Default

Quote:
Originally Posted by mbayer View Post
Hi,

personally I think that looks reasonable, assuming you have a eukaryotic organism -- the average gene length in eukaryotes is supposed to be in the 1,500 bp region. What organism is this, and do you know how common alternative splicing is in that species? That would obviously affect your number of transcripts relative to the number of genes. Also, given that most alternative splicing produces transcripts that are shorter than the full length mRNA, an average transcript length of 1,112 seems reasonable.

To evaluate your assembly I would run something on the transcripts that predicts proteins, like getorf from the EMBOSS suite of tools, then select the longest predicted protein and BLAST this against related species. This will give you an idea of how good your assembly is.

cheers

Micha
Hi mbayer

Thanks a lot for your response.

The organism that i am working on is Brassica rapa plant which is very close to model plant Arabidopsis. When i got my final trinity output i didn't realize that they include alternative splicing transcripts. Anyway during last few days i learnt a lot of how to do downstream analysis and this is something i planned to do to get the final transcripts (i mean best transcripts)

1. expression based: after running the abundance estimation (bowtie-express), retain those that have some minimum FPKM value (such as 1).

2. run the ORF extraction pipeline included in Trinity (don't restrict it to complete ORFs, get both complete and partials) - retain those that encode long ORFs (eg. 200 aa)

3. blastx the trinity transcripts against RefSeq, retain those that have homology to known proteins (E<=1e-10)

Take the union of {1,2,3} above and call it 'best'.

Do you further comment on this?
upendra_35 is offline   Reply With Quote
Old 10-22-2012, 01:16 AM   #4
mbayer
Member
 
Location: Dundee, Scotland

Join Date: Mar 2009
Posts: 29
Default

Hi,

2 and 3 sound reasonable. As to point 1), I wouldn't exclude transcripts on the basis of being lowly expressed -- you may end up removing genuine transcripts from your final set. Remember that some transcripts really are expressed at very low levels, and also that Illumina sequencing contains an element of randomness which means that at the lower end of the expression range there may be transcripts that actually were present in the sample at very low levels but have not been caught by the sequencing/and or data analysis.

cheers

Micha
mbayer is offline   Reply With Quote
Old 10-23-2012, 07:51 AM   #5
Jan_R
Junior Member
 
Location: Seattle

Join Date: Jul 2011
Posts: 8
Default

Quote:
Originally Posted by upendra_35 View Post
The organism that i am working on is Brassica rapa plant which is very close to model plant Arabidopsis. When i got my final trinity output i didn't realize that they include alternative splicing transcripts.

Might be a bit late but anyway:
http://seqanswers.com/forums/showthread.php?t=19129

Most times the sequences in one comp(onent) do resemble each other. It is likley they are isoforms of the same gene.
At least that's my impression after sequence comparison by blastn and blastx. You can also put them in clustalw.
My experience is that we could easily amplify the few sequences we were interested in out of cDNA. So at least for these few examples trinity worked pretty nice for us.

A colleague of mine used different smaller sequences from one component to assemble larger fragments. The fragments he obtained from trinity all had similarity to a known gene but were much smaller than expected. He had to try two or three combinations of alignments before he got the expected fragment in his PCR.

So the number of components could give you rough impression about the total number of genes represented in your assembly.

Good luck with your assembly!
Jan_R is offline   Reply With Quote
Old 10-23-2012, 05:15 PM   #6
upendra_35
Senior Member
 
Location: USA

Join Date: Apr 2010
Posts: 102
Default

Quote:
Originally Posted by mbayer View Post
Hi,

2 and 3 sound reasonable. As to point 1), I wouldn't exclude transcripts on the basis of being lowly expressed -- you may end up removing genuine transcripts from your final set. Remember that some transcripts really are expressed at very low levels, and also that Illumina sequencing contains an element of randomness which means that at the lower end of the expression range there may be transcripts that actually were present in the sample at very low levels but have not been caught by the sequencing/and or data analysis.

cheers

Micha
Does that mean i should have more relaxed FPKM cut-off? What do you think would be ideal by the way?
upendra_35 is offline   Reply With Quote
Old 10-23-2012, 05:19 PM   #7
upendra_35
Senior Member
 
Location: USA

Join Date: Apr 2010
Posts: 102
Default

Quote:
Originally Posted by Jan_R View Post
Might be a bit late but anyway:
http://seqanswers.com/forums/showthread.php?t=19129

Most times the sequences in one comp(onent) do resemble each other. It is likley they are isoforms of the same gene.
At least that's my impression after sequence comparison by blastn and blastx. You can also put them in clustalw.
My experience is that we could easily amplify the few sequences we were interested in out of cDNA. So at least for these few examples trinity worked pretty nice for us.

A colleague of mine used different smaller sequences from one component to assemble larger fragments. The fragments he obtained from trinity all had similarity to a known gene but were much smaller than expected. He had to try two or three combinations of alignments before he got the expected fragment in his PCR.

So the number of components could give you rough impression about the total number of genes represented in your assembly.

Good luck with your assembly!
Thanks for the info. Yes i did worked out on this a bit but in the manner that your colleague work. I have a feeling by myself that the Trinity transcripts are smaller when compared to known genes and probably i might have do the same as your colleague. Do you know if he had written any script to do this? Thanks in advance...
upendra_35 is offline   Reply With Quote
Old 10-24-2012, 02:29 AM   #8
Jan_R
Junior Member
 
Location: Seattle

Join Date: Jul 2011
Posts: 8
Default

Quote:
Originally Posted by upendra_35 View Post
Thanks for the info. Yes i did worked out on this a bit but in the manner that your colleague work. I have a feeling by myself that the Trinity transcripts are smaller when compared to known genes and probably i might have do the same as your colleague. Do you know if he had written any script to do this? Thanks in advance...
Seems to be Murphy's law that people are on vacation when one comes up with questions

He was specificaly interested in one specific gene so he didn't write a script for that. He just put them in a classic sequence alignment program (clustalw) together with a sequence from a closely related species.
After that he did an ordinary blast with the sequence of the related species and the reads. Best matching reads he also aligned onto the sequence to close gaps. So he more or less refined the work already done by Trinity.
I lack the experience to judge wether this was a wild guess but this approach worked out for him and he could amplify his gene at the second or third attempt.

But that is the part that bothers you when you already have a specific gene you're interested in.

I got small fragments assembled when the transcript just was not that abundant which makes perfectly sense.

What do you have in mind with your 'best of all times'-collection?
Jan_R is offline   Reply With Quote
Old 10-29-2012, 01:55 AM   #9
mbayer
Member
 
Location: Dundee, Scotland

Join Date: Mar 2009
Posts: 29
Default

Quote:
Originally Posted by upendra_35 View Post
Does that mean i should have more relaxed FPKM cut-off? What do you think would be ideal by the way?
I would not set a cut-off at all. Like I said, I don't think you can assume that a low expression level means there is an artefact -- it could be real.

Micha
mbayer is offline   Reply With Quote
Old 06-03-2015, 08:22 PM   #10
kurban910
Member
 
Location: urumqi

Join Date: Jul 2014
Posts: 58
Default

hi guys,
i assembled the raw reads and got trinity.fasta file, then i got the basic statistics by using TrinityStats.pl , but i also wanna get little bit more detialed statistical results like length distribution of transcripts in the fasta file with an image if its possible.
any suggesion for that?
thanks.
kurban910 is offline   Reply With Quote
Old 06-04-2015, 02:17 AM   #11
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,060
Default

Quote:
Originally Posted by kurban910 View Post
hi guys,
i assembled the raw reads and got trinity.fasta file, then i got the basic statistics by using TrinityStats.pl , but i also wanna get little bit more detialed statistical results like length distribution of transcripts in the fasta file with an image if its possible.
any suggesion for that?
thanks.
Try stats.sh from BBMap: http://seqanswers.com/forums/showthread.php?t=43529
GenoMax is offline   Reply With Quote
Old 06-04-2015, 09:34 AM   #12
kurban910
Member
 
Location: urumqi

Join Date: Jul 2014
Posts: 58
Default

thanks for the link @GenoMax
kurban910 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 01:43 AM.


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