This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Alignng RNAseq data for gene prediction

Hi,

I have rna-seq data for five tissue samples and I am planning to use these for gene prediction. Should I align them all together in one step to the genome or each tissue separately and later merge the bams??

Thanks, Karthic

gene prediction rna-seq

Alternative plan:

  • make a genome-guided transcriptome assembly (e.g. with Trinity)
  • use the generated transcripts as evidence for gene prediction

I have always known Trinity as de-novo transcriptome assembler.

Thanks. I will go through that.

Should I align them all together in one step

If you have a reference available then you are likely not doing gene prediction. And if you don't then you should be assemblig the data as suggested by @Michael.

We have assembled the genome and there is no other annotation available for this species. We have rnaseq and isoseq for some tissues. Currenlty figuring out the way I should prepare the files.

Was the genome assembly done independent of the RNAseq data? What do you mean by "prepare the files"?

Yes, independent of the RNAseq data. I mean gathering the evidence for the gene prediction by utilizing the RNAseq and isoseq data.

The RNAseq data should be assembled with trinity and transcripts to be given as input to tools like augustus, genemark etc or they should be mapped to genome with tools like hisat2/tophat and generate models with stringtie/cufflinks and later given as input to augustus, genemark etc.

1 answer

I think you can map everything at once. If you are worried that you might loose isoform information - you can also do it separately for each tissue and merge only at the end. I had pretty good experience with GAWN for that type of analysis - check it out.

Thank you. I will check it out.

Log in to answer this question.