Question: Strand-option in hisat2
1
gravatar for vm.higareda
21 months ago by
vm.higareda20
vm.higareda20 wrote:

Hello

I am mapping my rna-seq reads using hisat2. My library is single-end strand-specific (dUTP) illumina-TrueSeq-Stranded.

I am confused about what option is correct. Hisat2 manual say:

"For single-end reads, use F or R. 'F' means a read corresponds to a transcript. 'R' means a read corresponds to the reverse complemented counterpart of a transcript"

is correct if I choose "F" option??. if I visualize my alignments in IGV for example, all my reads must be in "red color" ?

what means if i have blue and red?

hisat2 strand rna-seq • 1.9k views
ADD COMMENTlink modified 21 months ago • written 21 months ago by vm.higareda20

Yes, reads are arrow-shaped. and I have in both directions. if my library was dUTP, the arrows will point at the left? or right?enter image description here

I hope you can watch the image

ADD REPLYlink written 21 months ago by vm.higareda20

Please use ADD COMMENT or ADD REPLY to answer to previous reactions, as such this thread remains logically structured and easy to follow. I have now moved your post but as you can see it's not optimal. Adding an answer should only be used for providing a solution to the question asked.

ADD REPLYlink written 21 months ago by WouterDeCoster37k

Well, you selected a quite bad example to look at because this looks like overlapping genes on the (very gene-dense) chrMT. Select a "normal" locus and look again, for example, the GAPDH gene.

if my library was dUTP, the arrows will point at the left? or right?

That depends on the direction of your gene. If your gene's start is left and it's end is right, your reads will go from right to left. Reverse stranded.

ADD REPLYlink written 21 months ago by WouterDeCoster37k

What do you think about that

https://ibb.co/iFfbJQ

Thank you

ADD REPLYlink written 21 months ago by vm.higareda20

Oh, I wasn't aware that you are not working on human samples, so Gapdh1 turns out to be not such a good example as I thought. Next time you ask a question, don't forget to mention the organism you are working on.

Nevertheless, look at the left-most reads in the screenshot.

  • The arrows on the gene go from left to right
  • The arrows from the reads go from right to left

Same for Gapdh1:

  • The gene goes from right to left
  • The reads (red) go from left to right.

So, therefore, it's indeed reverse stranded.

ADD REPLYlink written 21 months ago by WouterDeCoster37k

Yes, I forgot say that is D. melanogaster.

I agree with you, when the gene goes from right to left, reads go from left to right, not only for Gapdh1. However as in the post, seem to be that option F or R do not have effect on the result.

I mapped my reads using F or R and the IGV visualizations are the same

What did you do at the end?.

I am still confused about what option I have to use in hisat2 to my rna-seq data, F or R?

On the other hand, if I count the reads with the XS attribute tag in the sam outputs using F or R options, the values are different, example:

hisat2 -q -p 30 --rna-strandness F -x ../../indice/drosophila -U 40trimo.fastq -S 40_F_strand.sam

grep -o XS:A:+ 40_F_strand.sam| wc -l

25656117

hisat2 -q -p 30 --rna-strandness R -x ../../indice/drosophila -U 40trimo.fastq -S 40_R_strand.sam

grep -o XS:A:+ 40_R_strand.sam| wc -l

41833200

I am worried because the next step is run cuffdiff with the hisat2 outputs, and their results completely different if I use "F" in the --rna-strandness parameter or if I use the default option ( not strandness). I have not try with the outputs of --rna-strandness F.

ADD REPLYlink modified 21 months ago • written 21 months ago by vm.higareda20

So apparently it does make a difference.
Since you reads are reversely stranded, you should use --rna-strandness R

ADD REPLYlink written 21 months ago by WouterDeCoster37k

Thank you for your comments I am very grateful.

ADD REPLYlink written 21 months ago by vm.higareda20

https://ibb.co/mJdDyQ

ADD REPLYlink modified 21 months ago • written 21 months ago by vm.higareda20

By the way, see also my reply here: A: Hisat R strandness aligning same strand as F strandness, how to do it correctly?

Perhaps things have changed, but this parameter used to do nothing in my tests.

ADD REPLYlink written 21 months ago by WouterDeCoster37k
0
gravatar for WouterDeCoster
21 months ago by
Belgium
WouterDeCoster37k wrote:

I'm not sure if there is a relation between colours and direction, and as far as I know, you can choose this colours how you want them. But you can see in which direction your reads go, with respect to genes. Reads are arrow-shaped, right? Don't know how your IGV is configured, obviously.

Anyway, reads generated using dUTP method are reverse stranded.

ADD COMMENTlink written 21 months ago by WouterDeCoster37k
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: 1762 users visited in the last hour