nf-core / scrnaseq

A single-cell RNAseq pipeline for 10X genomics data
https://nf-co.re/scrnaseq
MIT License
214 stars 172 forks source link

smartseq 1/2 support #258

Closed matq007 closed 1 year ago

matq007 commented 1 year ago

Description of feature

Hi, I would like to add smartseq and smartseq2 support to the pipeline. StarSolo supports it which means it's simple to adjust it. From quick night session it requires only 2 files to be changed.

and by adding custom configuration

process {
    withName: STAR_ALIGN {
        ext.args = "--readFilesCommand zcat --soloUMIdedup Exact --soloStrand Unstranded"
    }
}

From my understanding this should be enough to make this work. I can implement and submit a PR.

grst commented 1 year ago

How would the samplesheet look like in this case? Since with smartseq every cell is sequenced independently, wouldn't this lead to one sample per cell in the end? We likely would need a different logic for merging biological samples in that case.

In the past we have recommended to use the rnaseq pipeline for smart-seq2 data. Is there any reason why this doesn't work for you?

nasjr08 commented 4 months ago

Is this now an option for this nfcore or is the recommendation still to use the bulk rnaseq nfcore?

grst commented 4 months ago

rnaseq is still the way to go for SmartSeq1/2.