SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Pacific Biosciences



Similar Threads
Thread Thread Starter Forum Replies Last Post
Pacbio sequel for single cell human WGS? WhatsOEver Pacific Biosciences 5 06-27-2017 12:02 PM
Sequel Sequencing Costs - Microbial Genomes docbio Pacific Biosciences 9 06-22-2017 12:17 PM
Enter our Microbial SMRT Grant Program - Win SMRT Sequencing on the Sequel System pacbio Pacific Biosciences 0 06-06-2017 09:35 AM
PacBio Sequel demultiplexing Ali May Pacific Biosciences 4 12-06-2016 10:43 PM
Sequel SMRT Cell chip SEQanalyst532 Pacific Biosciences 2 11-11-2016 12:15 PM

Reply
 
Thread Tools
Old 06-21-2017, 10:27 PM   #1
weigrc
Member
 
Location: Hong Kong

Join Date: Oct 2011
Posts: 43
Default Problem of demultiplexing 4 barcoded E.coli genomes in one Sequel SMRT cell

Dear all,

As captioned, this sequencing included 4 barcoded libraries prepared by the following protocol, and barcode lbc01, lbc02, lbc03 and lbc04 were chosen, respectively. The library size is around 7.5kb (the longest subread is 5259 bp), and the SMRT Link we use is Version 4.

http://www.pacb.com/wp-content/uploa...Sequencing.pdf


Refer to attached, demultiplexing the data with full set of 96 barcodes, the problem is that most of the sequencing reads are sorted as other 92 no-existing barcodes, but the 4 barcodes we used!!!!

I think the size as 5,259 bp is too long for PacBio demultiplex.


Anyone has experience on sequencing barcoded libraries?? suggestions?


Thanks,
Wei
Attached Files
File Type: pdf 4xBarcodedEcoli.pdf (220.5 KB, 18 views)

Last edited by weigrc; 06-22-2017 at 01:49 AM.
weigrc is offline   Reply With Quote
Old 06-22-2017, 05:05 AM   #2
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,550
Default

Have you contacted PacBio tech support directly about this?

While it is common to see some sequences getting assigned to indexes you did not use, you should obviously get the correct results, if the experiment worked as designed.
GenoMax is offline   Reply With Quote
Old 06-22-2017, 11:19 AM   #3
jharting
Junior Member
 
Location: CA Bay Area

Join Date: Jun 2017
Posts: 4
Default barcode orientation?

What barcode set did you use for barcode scoring? The v4.0 software is orientation specific. If you used barcodes in the wrong orientation, they will not be binned correctly.
jharting is offline   Reply With Quote
Old 06-22-2017, 06:51 PM   #4
weigrc
Member
 
Location: Hong Kong

Join Date: Oct 2011
Posts: 43
Default

We were using PacBio SMRTbell Barcoded Adapter Complete Prep Kit-96.

And we did contact our tech support to follow up this issue since this Tuesday with the full information they needed, but the answer I got so far is only PacBio does not recommend to do barcoded WGS.

I understand some reads are possibly mis-assigned to undersigned barcodes, but I worry about the possibility is much bigger than my expectation, I was told 1% or more may be normal for PacBio sequencing. Is it true?

Anyone has experience on barcoded WGS using Sequel system?


Thanks,
Wei

Last edited by weigrc; 06-23-2017 at 12:26 AM.
weigrc is offline   Reply With Quote
Old 06-23-2017, 09:17 AM   #5
jharting
Junior Member
 
Location: CA Bay Area

Join Date: Jun 2017
Posts: 4
Default

We have been doing multiplexed microbial assembly on Sequel with success. Barcoded adapters are for general use -- basically a smrtbell adapter with extra sequence, so there's no reason why they cannot be used for your purpose.

You should be using the barcodes labeled something like 'pacbio_barcode_adapters_96' for scoring them. That set is in the correct orientation. If you have used barcoded adapters previously with smrt analysis version 2.3 or earlier, you may have used a barcode file with the sequences in the other orientation without issue. That software was not orientation-specific, and so it did not matter. The barcoding tool was re-written for v3+, and it is now orientation-specific.
jharting is offline   Reply With Quote
Old 06-23-2017, 05:00 PM   #6
weigrc
Member
 
Location: Hong Kong

Join Date: Oct 2011
Posts: 43
Default

Thank you, jharting.

Do you know what algorithm SMRT Link is using for demultiplexing? I mean it requests CCS read or full-pass subread? Is it possible to access such setting?

The barcoded WGS library I've tested has 7.5 kb peak size with smear size range 3k-10k, and the polymerase read length shown in the attached. Although the official protocol I followed is for 10 kb library, I guess the 7.5 kb library is still no short enough for a good demultiplexing result. May I know what library size you have been targeting for barcoded microbial WGS?

Regards,
Wei
weigrc is offline   Reply With Quote
Old 06-26-2017, 10:36 AM   #7
jharting
Junior Member
 
Location: CA Bay Area

Join Date: Jun 2017
Posts: 4
Default

Hi Wei,

We typically target libraries of ~10-12kb for multiplexed assembly, although there is work to optimize this, including a size selection step to remove short reads and pre-extension runs where the laser is turned off for ~120min to let the polymerase get near the end of the insert before sequencing begins. The trade-off in multiplexed assembly on PB platforms is between length and barcode yield. You want long reads to generate nice, closed assemblies, but too long and the polymerase is less likely to reach a smrtbell adapter/barcode before dying, in which case you may get a nice long read but no barcode for binning!

The underlying tool in smrtlink for barcoding is 'bam2bam'. You can find some information about the tool in the smrt tools reference guide:
http://www.pacb.com/wp-content/uploa...ide-v4.0.0.pdf

This tools works only on the raw subreads. To label a zmw (sequencing hole) with a barcode, at least one smrtbell adapter must be included in the hqregion of the polymerase read for the zmw.

Here also is a link to a "third-party" demultiplexing tool that you can use, if for nothing else than to look at the demultiplexing report, which can be useful for diagnosing your dataset.

https://github.com/armintoepfer/uhu

Although this is technically "third-party" because is exists outside of pacbio an is not supported, its written by the same developer who wrote bam2bam. I find it useful and fast.
jharting is offline   Reply With Quote
Old 06-26-2017, 05:28 PM   #8
weigrc
Member
 
Location: Hong Kong

Join Date: Oct 2011
Posts: 43
Default

Hi Jharting,

Thank you for these information, and I am reading the documents, they seem useful. thanks!
weigrc is offline   Reply With Quote
Old 09-19-2017, 01:35 AM   #9
CarambaKaracho
Junior Member
 
Location: Basel

Join Date: Nov 2013
Posts: 4
Default

Hi weigrc and jharting, just wanted to thank you for the concise information regarding pacbio multiplexing. Your excellent questions and answers helped me a lot to understand what I'm getting from a sequencing provider.
CarambaKaracho is offline   Reply With Quote
Old 09-19-2017, 08:58 AM   #10
jharting
Junior Member
 
Location: CA Bay Area

Join Date: Jun 2017
Posts: 4
Default

Here is a beta version of the next generation barcoding tool for pacbio if you are interested in trying. It has a lot of useful options and a nicer user interface!

https://github.com/PacificBiosciences/barcoding
jharting is offline   Reply With Quote
Old 11-11-2017, 08:55 AM   #11
dho
Junior Member
 
Location: USA

Join Date: Sep 2017
Posts: 5
Default

I know this is a while after the last reply, but I struggled for much of the summer with a similar issue. Basically, if the barcode is not immediately adjacent to the SMRTBell adapter, the demultiplexer in SMRTLink behaves inconsistently. There are two solutions:

1) Make CCS files and demultiplex manually by looking for the sequences of your barcodes in the CCS files.

2) You can use the new PacBio Lima demultiplexer (the one mentioned by jharting) with the -W flag to set a larger window for barcode detection. For example, in my data, there is about 30bp of sequence between the adapter and 15bp barcode. Setting -W 50 allows barcodes to be detected correctly.

--dave
dho is offline   Reply With Quote
Reply

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 11:20 PM.


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