Exception: Event ENSG00000004961 not found in pickled directory index. Are you sure this is the right directory for the event?

[mdozmorov]

Description of the bug

This issue appears in various forms in other issues, but I could not find a workable solution. This error occurs when running the pipeline with --genome GRCh38 or with manually provided --fasta and --gtf (I use those from 10x Genomics).

-[nf-core/rnasplice] Pipeline completed with errors-
ERROR ~ Error executing process > 'NFCORE_RNASPLICE:RNASPLICE:VISUALISE_MISO:MISO_SASHIMI (1)'

Caused by:
  Process `NFCORE_RNASPLICE:RNASPLICE:VISUALISE_MISO:MISO_SASHIMI (1)` terminated with an error exit status (1)

Command executed:

  sashimi_plot --plot-event ENSG00000004961 index miso_settings.txt --output-dir sashimi

  cat <<-END_VERSIONS > versions.yml
  "NFCORE_RNASPLICE:RNASPLICE:VISUALISE_MISO:MISO_SASHIMI":
      python: $(python --version | sed "s/Python //g")
      misopy: $(python -c "import pkg_resources; print(pkg_resources.get_distribution('misopy').version)")
  END_VERSIONS

Command exit status:
  1

Command output:
  (empty)

Command error:
  INFO:    Converting SIF file to temporary sandbox...
  /usr/local/lib/python2.7/site-packages/matplotlib/cbook/deprecation.py:107: MatplotlibDeprecationWarning: The mpl_toolkits.axes_grid module was deprecated in version 2.1. Use mpl_toolkits.axes_grid1 and mpl_toolkits.axisartist provies the same functionality instead.
    warnings.warn(message, mplDeprecation, stacklevel=1)
  Traceback (most recent call last):
    File "/usr/local/bin/sashimi_plot", line 11, in <module>
      sys.exit(main())
    File "/usr/local/lib/python2.7/site-packages/misopy/sashimi_plot/sashimi_plot.py", line 276, in main
      plot_label=plot_label)
    File "/usr/local/lib/python2.7/site-packages/misopy/sashimi_plot/sashimi_plot.py", line 142, in plot_event
      %(event_name, pickle_dir)
  Exception: Event ENSG00000004961 not found in pickled directory index. Are you sure this is the right directory for the event?
  INFO:    Cleaning up image...

Work dir:
  /lustre/home/mdozmorov/data/WorkData/Tony/Krista/RPE1_GFP_MYCN_DAC_72hrs/work/fe/1da2e608fc7aac0f8f63541f62c08b

A workaround is to set --sashimi_plot false. But is there a better solution?

Command used and terminal output

export NXF_OPTS='-Xms8g -Xmx24g'
export NXF_SINGULARITY_CACHEDIR=/lustre/home/mdozmorov/singularity_cache
DIRIN=/lustre/home/mdozmorov/data/WorkData/
INPUT=${DIRIN}/samplesheet.csv
DIROUT=${DIRIN}/OUT_full
GENOME=/lustre/home/mdozmorov/data/ExtData/10x/refdata-gex-GRCh38-2024-A/fasta/genome.fa
GTF=/lustre/home/mdozmorov/data/ExtData/10x/refdata-gex-GRCh38-2024-A/genes/genes.gtf

nextflow run nf-core/rnasplice \
   --input ${INPUT} \
   --contrasts contrastsheet.csv \
   --fasta ${GENOME} \
   --gtf ${GTF} \
   --outdir ${DIROUT} \
   -profile singularity \
   --sashimi_plot false \
   -resume

Relevant files

No response

System information

nextflow version 23.10.0.5891, run on HPC in local mode, Singularity, Rocky Linux 9.4, nf-core/rnasplice: 1.0.4