Summarizing transcript level annotation from Trinity de novo assembly to gene level annotation
Entering edit mode
5 months ago
MB ▴ 40

Dear community,

I have a general question about differential expression analysis when one is working with non-model organisms. When genomic resources are absent, one of the first RNAseq analysis steps is the generation of the de novo transcriptome, which is used afterwards as reference for mapping/abundance estimation of the assembled transcripts. When these abundances are fed into differential expression analysis tools (e.g. DESEq2), it is recommended to use gene level estimates instead of using transcript/isoform abundances. When I am working with RNAseq data, I normally use: Trinity -> RSEM -> tximport/DESeq2, whereby I use the Trinity assumption about gene-isoform relation for summarizing abundances of gene level with tximport. This is of course not perfect, but it is a starting point. The problem comes, when you want to do functional enrichment: the transcriptome will be annotated on transcript level, but the differential gene expression is not on transcript level anymore, but for functional enrichment (e.g. GO term enrichment), you need annotation data and information about what is differentially expressed. For isoforms that are annotated to the same gene product (as in most cases), this is not problematic. But how to deal with isoforms from the same (Trinity) 'gene' which show different annotations (and would therefore get different GO terms)?

  • Going back to transcript level expression analysis (to avoid the annotation problem)?
  • Doing the annotation not on transcript level but select one representative isoform? If yes, how? Clustering, select the longest or most expressed isoform or other kinds of representatives?
  • Summarizing all isoform annotations (and therefore all GO terms) for one gene?

Are there any experiences/thoughts/recommendations about that? A link to overseen posts/literature would also be highly appreciated! Many thanks in advance, looking forward to your opinions/thoughts about that.

functional expression differential RNAseq transcriptome enrichment Trinity annotation genes • 355 views
Entering edit mode

Hello, I am wondering the exact same thing. I have been mulling over very similar options and sort of taking every approach, but I would love to hear from someone experienced with this situation.

Entering edit mode

MB did you get any resolution elsewhere in how to approach?


Login before adding your answer.

Traffic: 1080 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6