Skip to content

Releases: epi2me-labs/wf-somatic-variation

v1.4.1

16 Dec 10:40
Compare
Choose a tag to compare

Changed

  • Reconcile workflow with wf-template v5.3.3

v1.4.0

16 Sep 16:10
Compare
Choose a tag to compare

Added

  • BED files for the VNTR regions in hg19/hg38 from Severus.
    • The appropriate file will be automatically selected for the appropriate genome, unless a user provides a custom bed with --tr_bed.
  • IGV configuration supporting VCF files from --snv and --sv.
  • Aligned BAM to the output, when re-alignment is required.

Changed

  • Tweaked parameters for Severus to refine SV calling.
  • More informative log when a BAM called with an invalid basecaller model is provided.
  • Failures of processes involved in differentially modified loci and regions detection will not cause workflow to fail.
  • Updated modkit to v0.3.3.
  • Reconciled workflow _ingress.nf from wf-human-variation v2.4.0 and wf-template v5.2.6.

Fixed

  • -resume failing for some snv processes.
  • Excessive memory usage for sample_probs process when using --mod leading to exit code 137.
  • Erroneous handling of inputs for the joint report.
  • Missing re-aligned XAM files in the output directory.

v1.3.1

14 Aug 10:07
Compare
Choose a tag to compare

Added

  • --override_basecaller_cfg parameter allows users to provide a basecall configuration name in cases where automatic basecall model detection fails.
  • --diff_mod option to allow users to turn off differential modified loci (DMR) and regions (DMR) analysis with DSS by setting --diff_mod false.

Changed

  • Updated modkit to v0.3.0.
  • Reconciled workflow _ingress.nf from wf-human-variation v2.3.1.

Fixed

  • --snv crashing with --include_all_ctgs true.
  • --snv always referring to model in --override_basecaller_cfg when deciding whether to call Indels.
  • Automated basecaller detection not finding a basecaller model.

Removed

  • --basecaller_cfg as the workflow now automatically detects the basecaller model from the input data.

v1.3.0

11 Jul 08:26
Compare
Choose a tag to compare

Added

  • Tumor-only mode for the base workflow, small variant calling with ClairS-TO, modified base aggregation and somatic SV calling.

Changed

  • If available basecaller_cfg will be inferred from the basecall_model DS key of input read groups.
    • Providing --basecaller_cfg will not be required if basecall_model is present in the DS tag of the read groups of the input BAM.
    • basecaller_cfg will be ignored if a basecall_model is found in the input BAM.
    • The workflow will fail if the tumor and normal BAM files have not been called with the same basecall_model.
  • Updated to Severus v1.1.

Fixed

  • Workflow crashing when the input BED file has overlapping intervals.
  • Returning error in annotate_sv when END position smaller than POS.
  • Workflow not starting in nextflow v24.04.

v1.2.2

17 Jun 16:34
Compare
Choose a tag to compare

Changed

  • Update to ClairS v0.2.0.
    • The workflow now uses normal heterozygote sites to haplotag reads.
    • This behaviour can be changed using --use_normal_hets_for_phasing and --use_tumor_hets_for_phasing.
    • Moreover, it uses the indels for the phasing as default; this can be changed with --use_het_indels_for_phasing.
    • The workflow now uses longphase to haplotag reads; this can be changed with --use_longphase_haplotag.
    • The workflow now accepts a --liquid_tumor option, enabling presets and, where available, models specific for liquid tumors.
  • Update to Clair3 v1.0.8.
    • Added --clair3_base_err and --clair3_gq_bin_size options.
  • modkit now runs by contig.
  • modkit bedMethyl are now in the top level output directory.

v1.2.1

05 Jun 10:22
Compare
Choose a tag to compare

Added

  • A report with name [sample name].wf-somatic-variation-report.html, linking the individual detailed reports.

Changed

  • Use ezcharts SeqCompare to in QC report.
  • Memory usage of alignment report reduced by using histograms.
  • Retry process when clairs.py predict crashes with error 134.

v1.2.0

21 May 14:07
Compare
Choose a tag to compare

Added

  • Support for input folders of BAM files for --bam_tumor and --bam_normal (instead of only allowing single BAM files).

Changed

  • ClinVar version in SnpEff container updated to version 20240307
  • Update to Clair3 v1.0.7.
  • Update to modkit v0.2.6.
  • Improved modkit runtime by increasing the default interval size.
  • Increased minimum CPU requirement for the workflow to 16.
  • bedMethyl output files now follow the pattern {{ alias }}.wf-somatic-mods.{{ type }}.bedmethyl.gz.
  • Structural variant (SV) calling is now performed with Severus (v0.1.2).
  • ARM-compatible base workflow and modified base calling.
  • minimap2 alignments will be in BAM format when --sv is set.

Fixed

  • Force minimap2 to clean up memory more aggressively. Empirically this reduces peak-memory use over the course of execution.
  • Workflow occasionally repeating QC analyses when resuming, even if successful.
  • Alignment report script using too much memory.
  • HAC models not recognised as valid.
  • Spurious bamstats crashes with implausible alignment information.

Removed

  • CRAM as supported input format.
  • Reference genome and its indexes from the output directory.
  • Insert classification and support for mismatching panel of control.
  • Options --qv, --classify_insert, --min_ref_support, --genotype_sv and --control_panel, as these are no longer used by the workflow.

v1.1.0

05 Mar 18:24
Compare
Choose a tag to compare

Changed

  • Updated ClairS to v0.1.7, with the new dorado 4KHz/5KHz HAC models.
  • Several performance improvements which should noticeably reduce the running time of the workflow
  • makeQCreport allows for one retry to prevent the workflow failing on report generation

v1.0.0

07 Feb 22:04
Compare
Choose a tag to compare

Added

  • Tumor-only mode for the base workflow and modified base aggregation.
  • --control_panel option to provide a non-matching control panel produced with nanomonsv merge_control.
  • Memory directives for every process.

Changed

  • Run ClairS haplotype_filter by contig.
  • --sv will no longer emit genotypes by default, and only save the same fields from nanomonsv.
    • Users can still request a genotype using --genotype_sv.
    • --min_ref_support defines the minimum number of REF-supporting reads to call a heterozygote site.

Fixed

  • --snv calling genome-wide variants when --bed is specified.
  • snv:makeReport crashing when --annotation false.

Removed

  • --annotation_threads removed as snpEff v5 uses only one thread.

v0.5.2

06 Dec 20:57
Compare
Choose a tag to compare

Changed

  • New documentation
  • Updated ClairS to v0.1.6
  • ClinVar annotation of SVs has been temporarily removed due to not being correctly incorporated. SnpEff annotations are still produced as part of the final SV VCF.

Fixed

  • rVersion retried also upon success

Removed

  • Default local executor CPU and RAM limits