Question: Overrepresented sequence in RNA-seq only in R2
1
gravatar for benjamin.saintpierre
3.6 years ago by
benjamin.saintpierre20 wrote:

Hi everyone,

After running some RNA-seq experiment, we have observed an over-represented sequence which do not correspond to any adapter used. It clearly seems to be a contaminant.

But my main question is why do I get this sequence only in one strand?

As I get the sequence in 2 400 000 reads, I don't get its reverse-complemented sequence. I don't know why, because if I get a contaminant sequence, I should get half of the reads with the sequence, and half with its complement, shouldn't I?

rna-seq alignment sequence • 1.7k views
ADD COMMENTlink modified 2.5 years ago by Biostar ♦♦ 20 • written 3.6 years ago by benjamin.saintpierre20

Have you tried blasting it?

ADD REPLYlink written 3.6 years ago by WouterDeCoster43k

Yes, it's 60S Rrna human, so a contamination is obvious. But what i can't explain is why is there only this: R2adapter - Sequence of 300 bases - R1adapter and none of this: R1adapter - Sequence of 300 bases - R2adapter The adapters must be set in a random way, no?

ADD REPLYlink written 3.6 years ago by benjamin.saintpierre20

May be first you map your reads against rRNA and tRNA sequences (preferably via bowtie), and discard the mapable reads. If remaining reads still show over-represented sequences then there is something weird, it could also be the case of chimeric transcripts.

ADD REPLYlink written 3.6 years ago by Manvendra Singh2.1k

Actually, i've already done this mapping, and i know it's form this sequence: Chain 5, Structure Of The H. Sapiens 60s Rrna, 5070 bp. The problem is not what the sequence is, but why is it represented only in the 5'-3' strand, and not in the other way. I mean, if adapters are setting randomly around the rna fragment, why do i only observe R2adapter-Sequence of 300 bases-R1adapter, and not R1adapter-Sequence of 300 bases-R2adapter?

ADD REPLYlink written 3.6 years ago by benjamin.saintpierre20

Which library prep method was used for the RNA-seq?

ADD REPLYlink written 3.6 years ago by WouterDeCoster43k

Script Seq, Illumina

Thanks for your answers/time!

ADD REPLYlink written 3.6 years ago by benjamin.saintpierre20

This kit generates highly directional RNA-Seq libraries from rRNA-depleted or poly(A)+ RNA from any animal, plant, or bacterial species. Ribo-Zero ribosomal RNA (rRNA) removal kits are recommended for use with ScriptSeq v2 to achieve whole-transcriptome, diverse, stranded libraries.

So it's strand specific and rRNAs are supposed to be removed.

ADD REPLYlink written 3.6 years ago by WouterDeCoster43k

So if it is strand specific, it means that I will observe always get the R2adapter in 5' of my sequence and never in 3' ? So this: R1adapter-Sequence of 300 bases-R2adapter can't exist ?

ADD REPLYlink written 3.6 years ago by benjamin.saintpierre20

Yes, from the FAQ of Script-seq:

Does the ScriptSeq process generate directional libraries? Yes. ScriptSeq kits make directional libraries such that the sequence generated using the Illumina Read 1 sequencing primer is that of the sense strand of the original RNA. Using the Illumina Read 1 sequencing primer, the first nucleotide read always corresponds to the 5´ nucleotide of the original fragmented RNA molecule.

How does the ScriptSeq method generate directional libraries? After rRNA depletion (e.g., using Ribo-Zero Kits), the RNA is reverse-transcribed using random hexamers that include a 5´ tail containing a unique tagging sequence. After cDNA synthesis, a terminal tagging oligonucleotide (TTO) that has a blocked 3´ end is annealed to the 3´ end of the cDNA. Following extension with a DNA polymerase, a second, unique tagging sequence (complementary to the TTO) is added to the 3´ end of the cDNA. The resulting di-tagged cDNA is amplified by PCR using primers that anneal to the two different tagging sequences on the cDNA.

ADD REPLYlink modified 3.6 years ago • written 3.6 years ago by WouterDeCoster43k

You should not have RXadapter - Sequence of 300 bases - RYadapter in first place, as the first adapter will be sequenced by the primer anyway. Have you checked with the technician who did the experiment what went wrong? something seems very fishy on wet-lab level!

ADD REPLYlink written 3.6 years ago by H.Hasani910
0
gravatar for benjamin.saintpierre
3.6 years ago by
benjamin.saintpierre20 wrote:

Thnaks for all your answers! It's the strand specific experiment information that i didn't have in the first time.

Thanks for your time!

ADD COMMENTlink written 3.6 years ago by benjamin.saintpierre20
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: 776 users visited in the last hour