Question: differential expression analysis without genome
0
gravatar for mxlsherry1992
8 months ago by
mxlsherry199230 wrote:

Dear all,

When I did the differential expression analysis for the RNA seq data (no reference genome), after Trinity assembly to got Trinity.fasta, the next step should be align_and_estimate_abundance.pl to map the reads for each samples to the Trinity.fasta file. But the question is that is we need to do cd-hit to remove the redundancy for the Trinity.fasta first, and then do align_and_estimate_abundance.pl ?

Thank you!

rna-seq • 181 views
ADD COMMENTlink modified 8 months ago by h.mon29k • written 8 months ago by mxlsherry199230
0
gravatar for h.mon
8 months ago by
h.mon29k
Brazil
h.mon29k wrote:

But the question is that is we need to do cd-hit to remove the redundancy for the Trinity.fasta first

One could, but there are probably better options to reduce redundancy from the assembly - have a look at the Trinity FAQ. I would consider first:

Using the supertranscripts method, which will produce a genome-like gene representation of the transcriptome assembly, you can then follow up with Differential Transcript Usage via SuperTranscripts.

You can also use the Trinity.fasta.gene_trans_map generated by Trinity to get "gene" counts in addition to the transcript counts.

After quantifying transcript abundance, filtering transcripts with low counts (can be applied in conjunction with the two methods above).

ADD COMMENTlink written 8 months ago by h.mon29k

Thanks for the helpful reply!!

ADD REPLYlink written 8 months ago by mxlsherry199230
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1812 users visited in the last hour