Skip to content

Commit

Permalink
continue documentation
Browse files Browse the repository at this point in the history
  • Loading branch information
fmalmeida committed May 10, 2024
1 parent f842cba commit adfda0f
Show file tree
Hide file tree
Showing 4 changed files with 14 additions and 16 deletions.
22 changes: 13 additions & 9 deletions docs/usage.md
Original file line number Diff line number Diff line change
Expand Up @@ -260,16 +260,20 @@ PBMC_10K_CMV,https://raw.githubusercontent.com/nf-core/test-datasets/modules/dat

#### Additional supporting files

Cellranger-arc needs a reference index directory that you can provide with `--cellranger_index`. Be aware, you can use
for cellranger-arc the same index you use for cellranger ([see](https://kb.10xgenomics.com/hc/en-us/articles/4408281606797-Are-the-references-interchangeable-between-pipelines)).
Yet, a cellranger-arc index might include additional data (e.g., TF binding motifs). Therefore, please first check if
you have to create a new cellranger-arc index ([see here](https://support.10xgenomics.com/single-cell-multiome-atac-gex/software/pipelines/latest/advanced/references) for
more information)
Cellranger multi needs a reference for GEX and VDJ analysis. They are calculated on the fly given the reference files (`--fasta` and `--gtf`) provided, but users can also provide their own with: `--cellranger_index` and `--cellranger_vdj_index`, for GEX and VDJ, respectively.

If you decide to create a cellranger-arc index, then you need to create a config file to generate the index. The pipeline
can do this autmatically for you if you provide a `--fasta`, `--gtf`, and an optional `--motif` file. However, you can
also decide to provide your own config file with `--cellrangerarc_config`, then you also have to specify with `--cellrangerarc_reference`
the reference genome name that you have used and stated as _genome:_ in your config file.
> When running cellranger multi, without any VDJ data, users can also skip VDJ automated ref building with: `--skip_cellrangermulti_vdjref`.
Cellranger multi also requires some additional data/information when utilizing some specific Feature Barcode Technologies:

TODO: talk about these:
gex_frna_probe_set = null
gex_target_panel = null
gex_cmo_set = null
fb_reference = null
vdj_inner_enrichment_primers = null
gex_barcode_sample_assignment = null
cellranger_multi_barcodes = null

## Running the pipeline

Expand Down
1 change: 0 additions & 1 deletion nextflow.config
Original file line number Diff line number Diff line change
Expand Up @@ -55,7 +55,6 @@ params {
skip_emptydrops = false

// CellrangerMulti parameters
cellranger_gex_index = null
cellranger_vdj_index = null
skip_cellrangermulti_vdjref = false // skip mkvdjref if not using VDJ data in cellranger/multi samplesheet
gex_frna_probe_set = null
Expand Down
4 changes: 0 additions & 4 deletions nextflow_schema.json
Original file line number Diff line number Diff line change
Expand Up @@ -313,10 +313,6 @@
"description": "",
"default": "",
"properties": {
"cellranger_gex_index": {
"type": "string",
"description": "Specify a pre-built cellranger index for GEX analysis"
},
"cellranger_vdj_index": {
"type": "string",
"description": "Specify a pre-built cellranger index for VDJ analysis"
Expand Down
3 changes: 1 addition & 2 deletions workflows/scrnaseq.nf
Original file line number Diff line number Diff line change
Expand Up @@ -77,7 +77,6 @@ workflow SCRNASEQ {
ch_universc_index = params.universc_index ? file(params.universc_index) : []

//cellrangermulti params
cellranger_gex_index = params.cellranger_gex_index ? file(params.cellranger_gex_index, checkIfExists: true) : []
cellranger_vdj_index = params.cellranger_vdj_index ? file(params.cellranger_vdj_index, checkIfExists: true) : []
ch_multi_samplesheet = params.cellranger_multi_barcodes ? file(params.cellranger_multi_barcodes, checkIfExists: true) : []
empty_file = file("$projectDir/assets/EMPTY", checkIfExists: true)
Expand Down Expand Up @@ -273,7 +272,7 @@ workflow SCRNASEQ {
ch_genome_fasta,
ch_filter_gtf,
ch_cellrangermulti_collected_channel,
cellranger_gex_index,
ch_cellranger_index,
cellranger_vdj_index,
ch_multi_samplesheet
)
Expand Down

0 comments on commit adfda0f

Please sign in to comment.