Skip to content

12801 bug in gsea module related with multiple gene sets processing - #12821

Open
brovolia wants to merge 4 commits into
masterfrom
12801-bug-in-gsea-module-related-with-multiple-gene-sets-processing
Open

12801 bug in gsea module related with multiple gene sets processing#12821
brovolia wants to merge 4 commits into
masterfrom
12801-bug-in-gsea-module-related-with-multiple-gene-sets-processing

Conversation

@brovolia

Copy link
Copy Markdown
Contributor

I narrowed the removal pattern and used quotation marks, so the white space will not break the processes

@brovolia brovolia linked an issue Aug 26, 2026 that may be closed by this pull request
2 tasks
Comment thread modules/nf-core/gsea/gsea/main.nf Outdated
Comment thread modules/nf-core/gsea/gsea/main.nf Outdated
Comment thread modules/nf-core/gsea/gsea/main.nf Outdated

@pinin4fjords pinin4fjords left a comment

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Ran this for real before going further — pulled the gsea-cli 4.3.2 container and ran the actual script against real GCT/CLS/CHIP data plus two split Hallmark GMT files.

  • Comma-joined -gmx file1.gmt,file2.gmt: succeeds, combines both files (Got gsets: 50), output dir is plain <rpt_label>.Gsea.<timestamp> — no bracketed name list.
  • Space-separated (what this module currently produces): also exits 0, but gsea-cli silently keeps only the first file and drops the rest (to parse>hallmark_part1.gmt< got: [hallmark_part1.gmt]). No crash.

So the real bug is silent data loss, not a crash, and join(',') fixes it. Couldn't reproduce the crash or the bracketed-name naming issue — details inline.

Comment thread modules/nf-core/gsea/gsea/main.nf Outdated
Comment thread modules/nf-core/gsea/gsea/main.nf Outdated
Comment thread modules/nf-core/gsea/gsea/main.nf Outdated
Comment thread modules/nf-core/gsea/gsea/main.nf Outdated
@brovolia
brovolia force-pushed the 12801-bug-in-gsea-module-related-with-multiple-gene-sets-processing branch from 8c1d769 to 2f890ad Compare August 31, 2026 14:02
@github-actions github-actions Bot added size/m and removed size/s labels Aug 31, 2026
@brovolia

Copy link
Copy Markdown
Contributor Author

I added some changes that we've discussed. But still I believe that the pipeline-level changes should also take place, eg in conf/modules.config, where it takes gene_sets.baseName, a single file, in assets/differentialabundance_report.qmd that renders one report tab per individual GMT file

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Bug in gsea module related with multiple gene sets processing

2 participants