Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Issue with 5' technologies with single-ends #16

Closed
TomKellyGenetics opened this issue Sep 6, 2023 · 2 comments
Closed

Issue with 5' technologies with single-ends #16

TomKellyGenetics opened this issue Sep 6, 2023 · 2 comments
Assignees
Labels
bug Something isn't working

Comments

@TomKellyGenetics
Copy link
Collaborator

TomKellyGenetics commented Sep 6, 2023

As mentioned by @Davidwei7 (#12) there are known issues with STRT-Seq. An ex-colleague from RIKEN has also contacted us about supporting this technique. I've suggested using the (undocumented) "SC5P-R1" chemistry to run Cell Ranger. However, the latest version of UniverSC has automated checks that require both R1 and R2 files to exist, assuming that all technologies are paired-end.

As a work-around, I suggest copying R1 to and R2 FASTQ file, hard-trimming the barcodes from R2 and running 3' paired-end chemistry in Cell Ranger with a custom barcode whitelist of STRT-seq presets.

The stable solution is to refactor the input checks for UniverSC to conditionally allow single-end FASTQ inputs (missing R2 files) for 5' scRNA-Seq data and the associate technologies (i.e., STRT-seq, SCRB-seq, ICELL8). I won't have time to do this for a while but I'll note here publicly that this is a known issue and is being addressed.

@TomKellyGenetics TomKellyGenetics self-assigned this Sep 6, 2023
@TomKellyGenetics TomKellyGenetics added the bug Something isn't working label Sep 6, 2023
@TomKellyGenetics TomKellyGenetics changed the title Issue with 5' Issue with 5' technologies with single-ends Sep 6, 2023
@TomKellyGenetics
Copy link
Collaborator Author

TomKellyGenetics commented Sep 19, 2023

Technology "strt-seq" has been updated to default to single-end 5' chemistry "SC5P-R1" and allows --read1 to give single-end FASTQ input without a ---read2 file input.

I've tested these parameters on data from:
GSE29087 (SRR212477) published by Islam et al. (2011): https://genome.cshlp.org/content/21/7/1160.full

In principle manual input of SC5P-R1 chemistry for any technology should allow the same parameters.

@TomKellyGenetics
Copy link
Collaborator Author

Release 1.2.7 resolves this.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

1 participant