--clip3pAdapterMMp has to contain 2 values to match the number of mates

[SuhasSrinivasan]

Description of the bug

Dear Developers,

Thank you for developing and maintaining the pipeline.

When following the instructions for 3′ digital gene expression assays there is a consistent issue when specifying the custom STAR arguments.

Pipeline completed with errors-
ERROR ~ Error executing process > 'NFCORE_RNASEQ:RNASEQ:ALIGN_STAR:STAR_ALIGN (D0_2)'

Caused by:
  Process `NFCORE_RNASEQ:RNASEQ:ALIGN_STAR:STAR_ALIGN (D0_2)` terminated with an error exit status (102)


Command executed:

  STAR \
      --genomeDir star \
      --readFilesIn input1/D0_2_trimmed_1_val_1.fq.gz input2/D0_2_trimmed_2_val_2.fq.gz \
      --runThreadN 12 \
      --outFileNamePrefix D0_2. \
       \
      --sjdbGTFfile GRCh38.primary_assembly.genome.filtered.gtf \
      --outSAMattrRGline 'ID:D0_2'  'SM:D0_2'  \
      --quantMode TranscriptomeSAM --twopassMode Basic --outSAMtype BAM Unsorted --readFilesCommand zcat --runRNGseed 0 --outFilterMultimapNmax 20 --alignSJDBoverhangMin 1 --outSAMattributes NH HI AS NM MD --outSAMstrandField intronMotif --alignIntronMax 1000000 --alignIntronMin 20 --alignMatesGapMax 1000000 --alignSJoverhangMin 8 --outFilterMismatchNmax 999 --outFilterType BySJout --outFilterMismatchNoverLmax 0.1 --clip3pAdapterSeq AAAAAAAA AAAAAAAA --quantTranscriptomeSAMoutput BanSingleEnd



  if [ -f D0_2.Unmapped.out.mate1 ]; then
      mv D0_2.Unmapped.out.mate1 D0_2.unmapped_1.fastq
      gzip D0_2.unmapped_1.fastq
  fi
  if [ -f D0_2.Unmapped.out.mate2 ]; then
      mv D0_2.Unmapped.out.mate2 D0_2.unmapped_2.fastq
      gzip D0_2.unmapped_2.fastq
  fi

  cat <<-END_VERSIONS > versions.yml
  "NFCORE_RNASEQ:RNASEQ:ALIGN_STAR:STAR_ALIGN":
      star: $(STAR --version | sed -e "s/STAR_//g")
      samtools: $(echo $(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*$//')
      gawk: $(echo $(gawk --version 2>&1) | sed 's/^.*GNU Awk //; s/, .*$//')
  END_VERSIONS

Command exit status:
  102

Command output:
  (empty)

Command error:

  EXITING because of fatal PARAMETER error: --clip3pAdapterMMp has to contain 2 values to match the number of mates.
  SOLUTION: specify 2values in --clip3pAdapterMMp
  Apr 21 02:26:45 ...... FATAL ERROR, exiting

Workaround

--extra_star_args "--clip3pAdapterMMp 0 0"

Command used and terminal output

nextflow run nf-core/rnaseq 
-profile singularity 
--save_reference \
--input /mnt/storage/projects/r2/rna_initial/fastqs/samples.csv \
--outdir  /mnt/storage/projects/r2/rna_initial/ \
--gtf /mnt/storage/projects/reference/rna/gencode.v47.primary_assembly.annotation.gtf \
--fasta /mnt/storage/projects/reference/rna/GRCh38.primary_assembly.genome.fa \
--extra_star_align_args "--alignIntronMax 1000000 --alignIntronMin 20 --alignMatesGapMax 1000000 --alignSJoverhangMin 8 --outFilterMismatchNmax 999 --outFilterMultimapNmax 20 --outFilterType BySJout --outFilterMismatchNoverLmax 0.1 --clip3pAdapterSeq AAAAAAAA AAAAAAAA" \
--extra_salmon_quant_args="--noLengthCorrection"

Relevant files

No response

System information

Version of nf-core/rnaseq: 3.18.0
Nextflow version: v24.10.5
Container engine: singularity-ce version 4.1.2-jammy
OS: Ubuntu 22.04.5 LTS x86_64
Executor: Local
Hardware: Desktop

[SuhasSrinivasan]

Hi @pinin4fjords
Thank you for reviewing this and the previous issues.

When using paired-end sequencing with Lexogen QuantSeq 3‘ mRNA-Seq FWD, the Custom STAR parameters in the docs need an update.

--extra_star_args "--clip3pAdapterMMp 0 0"

Possible values for clip3pAdapterMMp are 0, 1 and 2.

[pinin4fjords]

I don't see clip3pAdapterMMp in your initial command at all- did you specify it, or did STAR produce this error without that?

The documentation in STAR on this is confusing- as others have noted alexdobin/STAR#1260

[pinin4fjords]

Actually, I suspect the issue here (and in the other issue) is that the documentation was written for the typical 3' DGE use case, where data are expected to be single-end. So I'm disinclined to change the documentation for a paired-end use case.

Could you confirm which assay you're using please, and why your data are paired end?

[SuhasSrinivasan]

Hi,

Please see details in this comment: #1543 (comment)

I don't see clip3pAdapterMMp in your initial command at all- did you specify it, or did STAR produce this error without that?

Yes, STAR produced the error documented above.

[pinin4fjords]

As discussed above, docs are correct for the typical single-endedness for this tech. PRs welcome to describe the settings needed for paired end data, but for now things are correct as designed.

[SuhasSrinivasan]

Thank you!