Skip to content

Add subworkflow for gridss somatic structural variant calling - #12574

Open
imsarath wants to merge 7 commits into
nf-core:masterfrom
imsarath:subworkflow_gridss_somatic_svs_calling
Open

Add subworkflow for gridss somatic structural variant calling#12574
imsarath wants to merge 7 commits into
nf-core:masterfrom
imsarath:subworkflow_gridss_somatic_svs_calling

Conversation

@imsarath

@imsarath imsarath commented Aug 4, 2026

Copy link
Copy Markdown
Contributor

This PR adds bam_tumor_normal_somatic_structural_variant_calling_gridss, a subworkflow that calls somatic structural variants from a tumor/normal BAM pair. It chains the full GRIDSS workflow - optional targeted fragment extraction, per-sample preprocess, joint assembly and call, then gridss_somatic_filter. It emits the unfiltered joint VCF plus the all-somatic and high-confidence somatic call sets.

PR checklist

Closes #4498

  • This comment contains a description of changes (with reason).
  • If you've fixed a bug or added code that should be tested, add tests!
  • If you've added a new tool - have you followed the module conventions in the contribution docs
  • If necessary, include test data in your PR.
  • Remove all TODO statements.
  • Broadcast software version numbers to topic: versions - See version_topics
  • Follow the naming conventions.
  • Follow the parameters requirements.
  • Follow the input/output options guidelines.
  • Add a resource label
  • Use BioConda and BioContainers if possible to fulfil software requirements.
  • Ensure that the test works with either Docker / Singularity. Conda CI tests can be quite flaky:
    • For modules:
      • nf-core modules test <MODULE> --profile docker
      • nf-core modules test <MODULE> --profile singularity
      • nf-core modules test <MODULE> --profile conda
    • For subworkflows:
      • nf-core subworkflows test <SUBWORKFLOW> --profile docker
      • nf-core subworkflows test <SUBWORKFLOW> --profile singularity
      • nf-core subworkflows test <SUBWORKFLOW> --profile conda

@imsarath imsarath self-assigned this Aug 4, 2026
Comment thread subworkflows/nf-core/bam_somatic_sv_calling_gridss/tests/main.nf.test Outdated
when {
workflow {
"""
input[0] = channel.of([

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Shouldn't need the channel.of or channel.value at all.

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

I tried removing it completely, but input[0] and input[1] still require it.

ch_fasta_fai_bwaindex // channel: [ val(meta2), fasta, fai, bwa_index ]
ch_pondir // channel: [ val(meta3), [ pondir ] ]
ch_gridss_config // channel: [ val(meta4), [ gridss_config ] ]
val_target_bed // string: target BED for panel/exome data, null to run untargeted

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

This isn't a value?

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

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

Yes, it's a groovy path value for target_bed. I realized we could implement this in a cleaner way, so I've updated the code a bit. Let me know what you think of this approach!

@SPPearce SPPearce left a comment

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

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

Can you format the main.nf with each module argument on its own line too please

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

Projects

None yet

Development

Successfully merging this pull request may close these issues.

new subworkflow: distributed computing GRIDSS subworkflow

2 participants