Question: Hisat2 build index for genome + transcriptome files together
0
gravatar for bioinfo17
3.2 years ago by
bioinfo1730
bioinfo1730 wrote:

Hi,

Could anyone please explain how to use Hisat2 when both the genome and transcriptome is available for the same reference strain. The genome is annoated using maker and hence have the gene models in gff format. The transcriptome is also available but is present only as an assembled fasta format. Is it OK to merge both the genome + transcriptome files together and build an index?

Thanks for your time.

rna-seq • 1.8k views
ADD COMMENTlink modified 3.2 years ago by Devon Ryan96k • written 3.2 years ago by bioinfo1730
1
gravatar for Devon Ryan
3.2 years ago by
Devon Ryan96k
Freiburg, Germany
Devon Ryan96k wrote:

Merging the genome and transcriptome would result in largely useless output (you'd have mostly multimappers). Better would be to either align against the transcriptome or against the genome (having created the index with the GFF file).

ADD COMMENTlink written 3.2 years ago by Devon Ryan96k

Thanks Devon for your quick reply. The % alignment rate was very low (~2-20%) when the sample reads were mapped against the individual genome fasta file, however, the alignment % significantly increased (~60-80%) when the sample reads were mapped against the transcriptome alone. I am unsure which one to use for downstream differential expression analysis using Stringtie etc. ps: I did not use the gff for building index, only used either the genome or transcriptome fasta file for creating index.

Any suggestions please? Many thanks

ADD REPLYlink written 3.2 years ago by bioinfo1730

It sounds like the genome assembly isn't very good. Align against the transcriptome and don't use stringTie, but rather salmon.

ADD REPLYlink written 3.2 years ago by Devon Ryan96k

Will give Salmon a try, thank you!

ADD REPLYlink written 3.2 years ago by bioinfo1730
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: 1566 users visited in the last hour