5
$\begingroup$

I'm working with rna-seq samples. I see 5' bias and also 3' bias in the per-base sequence content plot. From this link I see that the bias at the start of the sequences appears to be the result of biased selection of fragments from the library. And this won't effect any downstream analysis. But In all of my samples I also see bias at the end of the sequences. What can be the problem? Samples are strand-specific sequenced with poly-A selection protocol.

In the Transcript coverage profile I see 5'-3' bias for all my samples is: 0.9

Showing the per-base sequence plot for one of the sample: For all the samples I see that 5' and 3' bias.

When I looked at Adapter content this is what I see:

No samples found with any adapter contamination > 0.1%

enter image description here

$\endgroup$
5
$\begingroup$

I expect there was a sequencing problem during the last base, where some of the reagents were running low on the sequencer. This won't pose any real problem, RNAseq aligners like STAR will just soft-clip the last base or two if they're mismatches.

It's common to see a bit of bias toward the 5' or 3' ends in RNAseq, mostly due to whether poly-A selection was done or there was a bit of degradation. The primary thing to ensure is that the bias is similar across samples/groups. If that's the case then you don't have to worry that your downstream analyses are going be affected by some sort of per-group bias.

$\endgroup$
  • $\begingroup$ I'm using hisat2 and yes I see the bias is similar across samples. Thanks for the answer. $\endgroup$ – user3351523 Feb 18 '18 at 14:39
  • 1
    $\begingroup$ Ah, for hisat2 you might get a bit better results if you trim off the last 1-2 bases. This is because hisat2 does global alignments (i.e., it doesn't do soft-clipping). Or at least that was the case last time I checked. $\endgroup$ – Devon Ryan Feb 18 '18 at 14:41
  • $\begingroup$ Cool. I will consider this and will trim few bases in the 3' end. Does it need any trimming in the 5' end also? $\endgroup$ – user3351523 Feb 18 '18 at 14:46
  • 1
    $\begingroup$ Nope, the 5’ end should be ok. $\endgroup$ – Devon Ryan Feb 18 '18 at 17:13
  • 1
    $\begingroup$ I thought tha hisat2 had a --local option, but when I run hisat2 -h | grep local, I only get --ma <int> match bonus (0 for --end-to-end, 2 for --local). No --local proper. The help does not look coherent. $\endgroup$ – bli Feb 19 '18 at 9:34

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.