SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Orientation of paired-end reads marnilush Illumina/Solexa 6 09-13-2012 08:19 AM
Orientation of 454 paired end reads split by linker skblazer 454 Pyrosequencing 8 04-26-2012 05:30 PM
should I reverse paired end reads before mapping? supermario Bioinformatics 2 03-06-2012 08:59 PM
454 Paired End orientation problems pr0t3us Bioinformatics 5 06-03-2010 12:25 AM

Reply
 
Thread Tools
Old 03-03-2012, 07:46 PM   #1
darthsequencer
Member
 
Location: california

Join Date: Feb 2012
Posts: 35
Default Paired end reverse orientation?

Hi,
I was given a lot of Illumina paired end reads. ~300bp inserts were used.

If I wanted to combine these into one long read where it would be:

Read 1 ----- Read 2

Where --- = the bases between the ends.

Do I just need to reverse Read 2? or reverse complement?

I think reverse and not complement.

Thanks
darthsequencer is offline   Reply With Quote
Old 03-04-2012, 09:36 AM   #2
mastal
Senior Member
 
Location: uk

Join Date: Mar 2009
Posts: 667
Default Paired end

I think you need to reverse complement.

My understanding is that the second read comes from the opposite strand to the first read.

Maria
mastal is offline   Reply With Quote
Old 03-04-2012, 11:50 AM   #3
darthsequencer
Member
 
Location: california

Join Date: Feb 2012
Posts: 35
Default

Ok thanks - if anyone else can confirm that'd be really helpful.

Thanks
darthsequencer is offline   Reply With Quote
Old 03-04-2012, 12:17 PM   #4
BAMseek
Senior Member
 
Location: St. Louis, MO, USA

Join Date: Apr 2011
Posts: 124
Default

Yep, the second read of an Illumina paired-end read would need to be reversed and complemented. Don't forget to reverse the quality string, too. You might want to check out FLASH, which will detect paired end reads that overlap each other and merge them into a single read.

Justin
BAMseek is offline   Reply With Quote
Old 03-04-2012, 12:28 PM   #5
darthsequencer
Member
 
Location: california

Join Date: Feb 2012
Posts: 35
Default

Thanks for the help guys.
darthsequencer is offline   Reply With Quote
Old 05-25-2012, 02:21 PM   #6
modi2020
Member
 
Location: New York

Join Date: May 2012
Posts: 22
Default

Hi guys...
I am a bit confused about the way we deal with paired end reads...
Say we have a read that is AAACGCGCTTTCCC
does paired end mean we only sequence AAA and CCC?
how about the sequence in between?
modi2020 is offline   Reply With Quote
Old 05-25-2012, 07:54 PM   #7
BAMseek
Senior Member
 
Location: St. Louis, MO, USA

Join Date: Apr 2011
Posts: 124
Default

Quote:
Originally Posted by modi2020 View Post
Hi guys...
Say we have a read that is AAACGCGCTTTCCC
does paired end mean we only sequence AAA and CCC?
how about the sequence in between?
In this case, AAACGCGCTTTCCC would be the sequence fragment. Ideally, you would just read the entire fragment, but due to the current limitations of the technology, that may not be possible. Paired-end sequencing allows you to read some bases from the beginning and some bases from the end of the fragment. If the fragment is short enough, the first read might overlap the second read. But in cases where the fragment length is longer than the sum of the paired read lengths, you would have a gap in the middle.

Let's say your read length was 4 base pairs. Then the first read would be AAAC and the second read would be GGGA (because the second read is obtained by reading the reverse complement).

That's at least my understanding of the process. Hopefully that is of some help.

Justin
BAMseek is offline   Reply With Quote
Old 05-25-2012, 09:42 PM   #8
modi2020
Member
 
Location: New York

Join Date: May 2012
Posts: 22
Default

Thank you BAMseek for helping me with this.
I have an inquiry though. Assuming the reads are 200 bp. suppose we read the ends only. How would we know the sequence in between?
Also, while sequencing (by synthesis) wouldn't we get a TTTG read for the AAAC sequence ?
I thought that in sequencing by synthesis we add the compliment of the base in the genomic sequence. Thank you
Quote:
Originally Posted by BAMseek View Post
In this case, AAACGCGCTTTCCC would be the sequence fragment. Ideally, you would just read the entire fragment, but due to the current limitations of the technology, that may not be possible. Paired-end sequencing allows you to read some bases from the beginning and some bases from the end of the fragment. If the fragment is short enough, the first read might overlap the second read. But in cases where the fragment length is longer than the sum of the paired read lengths, you would have a gap in the middle.

Let's say your read length was 4 base pairs. Then the first read would be AAAC and the second read would be GGGA (because the second read is obtained by reading the reverse complement).

That's at least my understanding of the process. Hopefully that is of some help.

Justin
modi2020 is offline   Reply With Quote
Old 05-27-2012, 08:29 PM   #9
BAMseek
Senior Member
 
Location: St. Louis, MO, USA

Join Date: Apr 2011
Posts: 124
Default

Quote:
Originally Posted by modi2020 View Post
Assuming the reads are 200 bp. suppose we read the ends only. How would we know the sequence in between?
You would have to rely on other reads that align nearby to get the bases in between. Your fragmentation/sonification size selection would determine how big of a gap you expect between read pairs.

Quote:
Also, while sequencing (by synthesis) wouldn't we get a TTTG read for the AAAC sequence ?
I thought that in sequencing by synthesis we add the compliment of the base in the genomic sequence.
Yes, that's correct. But also keep in mind that you might be reading the cDNA of the original sequence, and the sequence might have gone though a few rounds of PCR. If the protocol is not strand specific, then the strand information is lost anyway. I think the important thing to note that if it is Illumina sequencing, you would expect the strands of a read pair to be opposite of each other (forward-reverse or reverse-forward) in a "proper" pair.
BAMseek is offline   Reply With Quote
Old 03-13-2013, 06:16 AM   #10
tahamasoodi
Success
 
Location: India

Join Date: May 2012
Posts: 130
Default

Hi BAMseek,
I also have a doubt in this. If we read the ends (say 100 bp) only for a fragment of 500 bp, there will be around 300 bp in between which would not be read. You say we have rely on other reads that align nearby, which are these reads if we have only forward and reverse reads for a fragment?
__________________
Thanks,
tahamasoodi is offline   Reply With Quote
Old 03-13-2013, 06:34 AM   #11
kwaraska
Senior Member
 
Location: Boston,MA

Join Date: Nov 2008
Posts: 122
Default

The basis of this method of sequencing is that you have random shearing.

While one fragment will read bases 1-300, the next one might read 200-500, and a third may be 100-400. In this way, every base gets read, but on different fragments.
kwaraska is offline   Reply With Quote
Old 03-15-2013, 10:45 AM   #12
tahamasoodi
Success
 
Location: India

Join Date: May 2012
Posts: 130
Default

Thanks Kwaraska, it was really helpful!
__________________
Thanks,
tahamasoodi is offline   Reply With Quote
Reply

Tags
end, hiseq, illumina, paired

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:56 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO