Question: Weird Kmer at 5' of read 1 in paired-end RNA-seq
gravatar for izzy.yichao.cai
3.4 years ago by
izzy.yichao.cai170 wrote:

Hi all,

I've come across a weird issue with my RNA-seq data.

Some details about the data:

There 4 libraries of RNA, all total RNA-seq in wild-type and mutant cell with 2 replicate. Total 2 wild-type(Lane1, Lan2) + 2 mutant(Lane3, Lane4).

  1. Platform: Ilumina Truseq
  2. Prep kit: TruSeq LT Kits and TruSeq v1/v2 Kits
  3. Adaptor: Lane1--TruSeq Adapter, Index 4; Lane2--TruSeq Adapter, Index 6; Lane3--TruSeq Adapter, Index 12; Lane4--TruSeq Adapter, Index 12.

After I got the RNA-seq data, I ran FastQC.

In the result I found:

  1. Failed Per base sequence content --- might due to random hexamer priming bias (Commonly seen in RNA-seq data)
  2. Failed Kmer content --- Both in R1 and R2. The weird thing is that the over-represented Kmer are perfectly consistent at the 5' of R1. There is a perfect 12 bp sequence that at 5' of R1 : CAGACCGCGTTC. This only appears on R1 not in R2, and it is all across 4 libraries. Although R2 also has kmer at 5', their order is not consistent and the peak is somewhat random.  This how it looks like:

Lane1-R1 Lane1_R1

Lane1-R1 enter image description here

Lane2-R1 enter image description here

What I've done:

  1. Blast --- no result
  2. Check the adaptor sequence --- no match

Here is the all the image: (I cannot upload my image to free image hosting website. My images are saved via right click save as from the browser. This is weird.)

Similar issue(Despite that he managed to find the source of the kmer):

Anyone has any idea about what is going wrong here? Thank you so much for you input!

sequencing rna-seq next-gen • 1.4k views
ADD COMMENTlink modified 3.4 years ago • written 3.4 years ago by izzy.yichao.cai170

There is probably nothing wrong. This may represent signature pattern seen at the beginning of RNAseq reads. You should read this blog post by Dr. Simon Andrews, author of FastQC.

ADD REPLYlink modified 3.4 years ago • written 3.4 years ago by genomax89k

It's spelled 'weird' (sorry, but it drives me nuts...)

ADD REPLYlink written 3.4 years ago by cschu1812.5k

LOL. Noticed that already. Thanks for the correction!

ADD REPLYlink written 3.4 years ago by izzy.yichao.cai170

There's at least one more instance in your post ;P

ADD REPLYlink written 3.4 years ago by cschu1812.5k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 620 users visited in the last hour