Drop-seq cell barcode whitelist
1
2
Entering edit mode
2.1 years ago

G'morning,

I have been searching for a cell barcode (CB) whitelist for Drop-seq, but it seems difficult to find a standard list. On the other hand, it is relatively easy to find whitelists for 10x protocols.

If no whitelist is available, what would be the best procedure here? I have used STARsolo to align the reads so far, and it accepts a parameter configuration of --soloCBwhitelist None if no whitelist is available. For all intents and purposes, I presume that I can then just assume that the cell barcodes that STARsolo identifies from one sample is a whitelist that I can use across all samples?

I presume that other programs, e.g., UMI-tools and Alevin, can also take a FASTQ pair and be used here?

Happy Easter / 复活节快乐

Kevin

Drop-seq scrna-seq CB cell barcode • 1.9k views
ADD COMMENT
2
Entering edit mode

Not really answering your question but as an alternative you could use Alevin for the lowlevel processing. It has a --dropseq flag and (from what I understand) uses a ML approach to decide which CBs are present and trustworthy, so no need for an external list. That is then probably similar to what UMI-tools does since Alevin uses that internally from what I know. Rob and i.sudbery will know for sure.

ADD REPLY
1
Entering edit mode
ADD REPLY
0
Entering edit mode

Thank you, genomax

ADD REPLY
0
Entering edit mode

Dear Kevin,

Would you kindly tell me if you successfully apply STARsolo with Drop-seq? I also plan to learn how to use STARsolo with Drop-seq analysis. It would be very kind of you if you can provide me an example for such application.

Best regards, Kaj

ADD REPLY
1
Entering edit mode

Hey Kaj, if I recall correctly, I eventually used --soloCBwhitelist None, which instructs STAR to find cell barcodes automatically. There was a thread on GitHub: https://github.com/alexdobin/STAR/issues/584

ADD REPLY
0
Entering edit mode

Thank you very much!!

ADD REPLY
0
Entering edit mode

Dear Kevin,

Would you please enlighten me a little bit more? After STARsolo alignment, how could you acquire the expression matrix of all cells of a file? I have tried --quantMode GeneCounts, however it rendered me the overall gene counts of the sam file not gene count per cell.

Best regards Kaj

ADD REPLY
1
Entering edit mode

Hey Kaj, after STARsolo alignment, we then use featureCounts

ADD REPLY
1
Entering edit mode

Thank you very much!!

ADD REPLY
1
Entering edit mode
2.1 years ago

Drop Seq doesn't have a white list of barcodes. I think the process by which the barcodes are made doesn't restrict them to a white list.

I use umi_tools whitelist, and umi_tools will also help you pick where the cutoff is between good cell barcodes and empty ones.

ADD COMMENT
0
Entering edit mode

Thanks swbarnes2

ADD REPLY

Login before adding your answer.

Traffic: 2555 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6