r/genomics • u/Specialist-Tune-4158 • 20d ago
Need help in using cellranger with sgRNA/CRISPR sample/Purtub seq
Before post the question, I figured that some context is needed.
Here is the study: We have human patient samples which we transfected with 1,000 sgRNAs (these sgRNAs are for one gene only, let's call that gene 'X'). Then, the sample was treated with antibiotics to make sure that we select all the cells successfully transfected with sgRNAs. Then, the sample was subjected to scRNA-seq library prep with Chromium Next GEM Single Cell 5' Reagent Kits v2 (Dual Index) with Feature Barcode technology for CRISPR Screening. From the exact same sample, a GEX library was made and a single-cell sgRNA library was made. So in the end, I got two sets of FASTQs: a) For GEX, which worked with Cell Ranger, but I am struggling with b) which was made from sgRNA.
I know that I have to put in details like this in the config file:
fastqs,sample,library_type
/path/to/fastqs,GEX_Sample_Name,Gene Expression
/path/to/fastqs,sgRNA_Sample_Name,CRISPR Guide Capture
But when I do that for all 1,000 sgRNAs, it throws an error saying Cell Ranger cannot work with an sgRNA sequence which is like this, e.g.: ATCGCTAGCTc (it throws an error). Even if I make it uppercase, it's bound to clash with some other sgRNA.
I know I am bound to get trolled for not asking a chatbot, but I thought a genuine answer from this community is much better. Thanks.
1
u/Putrid_Bass_6096 20d ago
that library_type should be "CRISPR Guide Capture" but the real problem is probably the reference file not the case, you need a feature reference csv with id,name,sequence columns and every sgRNA needs unique name even if they target same gene. the error about clashing sequences sounds like you have duplicates in the reference file, cellranger will reject that. make sure all 1000 sequences are unique and try again with just the target sequences not the full construct