Merge remote-tracking branch 'upstream/dev'

CHANGELOG.md

@@ -3,9 +3,12 @@
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/)
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
-## Unreleased
+## Dev vTBD - TBD - TBD
 
-- Fix [[#488]](https://github.com/nf-core/smrnaseq/issues/488) - Fix runtime error that can occur when params.genome is set.
+- [[#481]](https://github.com/nf-core/smrnaseq/pull/481) - Fix [MIRTOP_STATS IndexError](https://github.com/nf-core/smrnaseq/issues/477) - Fix mirtop process execution when mirgenedb is used.
+- [[#482]](https://github.com/nf-core/smrnaseq/pull/482) - Update documentation regarding MirgeneDB input files.
+- [[#486]](https://github.com/nf-core/smrnaseq/pull/486) - Replace `CSVTK_JOIN` to improve processing in large amount of files.
+- [[#493]](https://github.com/nf-core/smrnaseq/pull/493) - Fix [[#488]](https://github.com/nf-core/smrnaseq/issues/488) - Fix runtime error that can occur when params.genome is set.
 
 ## v2.4.0 - 2024-10-14 - Navy Iron Boxer
 

bin/pivot_longer.R

@@ -0,0 +1,29 @@
+#!/usr/bin/env Rscript
+
+library(optparse)
+library(tidyr)
+library(vroom)
+
+option_list <- list(
+    make_option(c("--input"), type = "character", help = "Input TSV file", metavar = "character"),
+    make_option(c("--output"), type = "character", help = "Output CSV file", metavar = "character")
+)
+
+opt_parser <- OptionParser(option_list = option_list)
+opt <- parse_args(opt_parser)
+
+# Read CSV with vroom
+data <- vroom::vroom(opt$input, delim = "\t", col_types = c(.default = "c"))
+
+last_col <- names(data)[ncol(data)]
+
+# Convert from wide to long format
+long_data <- data %>%
+    pivot_longer(
+        cols = last_col,
+        names_to = "Sample_ID",
+        values_to = "Counts"
+    )
+
+vroom_write(long_data, opt$output, delim = ",")
+

bin/pivot_wider.R

@@ -0,0 +1,32 @@
+#!/usr/bin/env Rscript
+
+library(optparse)
+library(tidyr)
+library(vroom)
+library(dplyr)
+
+option_list <- list(
+    make_option(c("--input"), type = "character", help = "Input CSV file in long format", metavar = "character"),
+    make_option(c("--output"), type = "character", help = "Output CSV file in wide format", metavar = "character")
+)
+
+opt_parser <- OptionParser(option_list = option_list)
+opt <- parse_args(opt_parser)
+
+# Read CSV with vroom
+long_data <- vroom::vroom(opt$input, delim = ",",
+        col_types = c(
+        Counts = "d",
+        .default = "c"
+    ))
+
+# Transform to wide format
+wide_data <- long_data %>%
+    pivot_wider(
+        names_from = Sample_ID,
+        values_from = Counts,
+        values_fill = 0
+    )
+
+# Export wide format
+vroom_write(wide_data, opt$output, delim = "\t")

conf/modules.config

@@ -468,9 +468,11 @@ process {
         ]
     }
 
-    withName: 'NFCORE_SMRNASEQ:MIRNA_QUANT:CSVTK_JOIN' {
-        ext.args = "--fields 'UID,Read,miRNA,Variant,iso_5p,iso_3p,iso_add3p,iso_snp,iso_5p_nt,iso_3p_nt,iso_add3p_nt,iso_snp_nt' --tabs --outer-join --na \"0\" --out-delimiter \"\t\""
-        ext.prefix = "joined_samples_mirtop"
+    withName: 'NFCORE_SMRNASEQ:MIRNA_QUANT:PIVOT_LONGER' {
+        publishDir = [ enabled: false ]
+    }
+
+    withName: 'NFCORE_SMRNASEQ:MIRNA_QUANT:PIVOT_WIDER' {
         publishDir = [
             path: { "${params.outdir}/mirna_quant/mirtop" },
             mode: params.publish_dir_mode,

docs/usage.md

@@ -45,6 +45,10 @@ If MirGeneDB should be used instead it needs to be specified using `--mirgenedb`
 - `mirgenedb_mature`: This parameter should point to the FASTA file containing mature miRNA sequences. The file can be manually downloaded from [MirGeneDB](https://mirgenedb.org/download).
 - `mirgenedb_hairpin`: This parameter should point to the FASTA file containing precursor miRNA sequences. Note that MirGeneDB does not offer a dedicated hairpin file, but the precursor sequences can be downloaded from [MirGeneDB](https://mirgenedb.org/download) and used instead.
 
+See examples in: [the test-datasets repository of nf-core](https://github.com/nf-core/test-datasets/tree/smrnaseq/MirGeneDB).
+
+> [!NOTE] > `mirtop` is hard-coded to use the `pre` sequences, which originate from the hairpin FASTA, rather than the `pri` sequences, which come from the mature FASTA. Users must provide `pre` files from the start to ensure consistency between the FASTA and GFF files, as the coordinates in the GFF file are referenced to `pre` sequences. This also ensures that names in the BAM file will match those in the GFF.
+
 ### Genome
 
 - `fasta`: the reference genome FASTA file

modules/local/pivot/longer/environment.yml

@@ -0,0 +1,5 @@
+channels:
+  - conda-forge
+dependencies:
+  - conda-forge::r-optparse
+  - conda-forge::r-tidyverse

modules/local/pivot/longer/main.nf

@@ -0,0 +1,29 @@
+process PIVOT_LONGER {
+    tag"$meta.id"
+    label 'process_single'
+
+    conda "${moduleDir}/environment.yml"
+    container "community.wave.seqera.io/library/r-optparse_r-tidyverse_r-vroom:3cbb224fea84a0e1"
+
+    input:
+    tuple val(meta), path(tsv)
+
+    output:
+    tuple val(meta), path("*_long.csv") , emit: csv
+    path "versions.yml"                 , emit: versions
+
+    script:
+    """
+    pivot_longer.R \\
+        --input ${tsv} \\
+        --output ${meta.id}_long.csv
+
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        r-base: \$(echo \$(R --version 2>&1) | sed 's/^.*R version //; s/ .*\$//')
+        tidyr: \$(Rscript -e "library(tidyr); cat(as.character(packageVersion('tidyr')))")
+        optparse: \$(Rscript -e "library(optparse); cat(as.character(packageVersion('optparse')))")
+    END_VERSIONS
+    """
+
+}

modules/local/pivot/wider/environment.yml

@@ -0,0 +1,5 @@
+channels:
+  - conda-forge
+dependencies:
+  - conda-forge::r-optparse
+  - conda-forge::r-tidyverse

modules/local/pivot/wider/main.nf

@@ -0,0 +1,35 @@
+process PIVOT_WIDER {
+    tag"$meta.id"
+    label 'process_high'
+
+    conda "${moduleDir}/environment.yml"
+    container "community.wave.seqera.io/library/r-optparse_r-tidyverse_r-vroom:3cbb224fea84a0e1"
+
+    input:
+    tuple val(meta), path(csvs)
+
+    output:
+    tuple val(meta), path("*joined_samples_mirtop.tsv") , emit: csv
+    path "versions.yml"                                 , emit: versions
+
+    script:
+    """
+    awk 'NR == 1 || FNR > 1' ${csvs.join(' ')} > final_long_results_temp.csv
+
+    pivot_wider.R \\
+        --input final_long_results_temp.csv \\
+        --output ${meta.id}_concatenated_temp.csv
+
+    sort -t\$'\t' -k1,1 ${meta.id}_concatenated_temp.csv > joined_samples_mirtop.tsv
+
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        r-base: \$(echo \$(R --version 2>&1) | sed 's/^.*R version //; s/ .*\$//')
+        tidyr: \$(Rscript -e "library(tidyr); cat(as.character(packageVersion('tidyr')))")
+        dplyr: \$(Rscript -e "library(dplyr); cat(as.character(packageVersion('dplyr')))")
+        optparse: \$(Rscript -e "library(optparse); cat(as.character(packageVersion('optparse')))")
+        vroom: \$(Rscript -e "library(vroom); cat(as.character(packageVersion('vroom')))")
+    END_VERSIONS
+    """
+
+}

nextflow_schema.json

@@ -183,6 +183,13 @@
                     "help_text": "Saving generated references means that you can use them for future pipeline runs, reducing processing times.",
                     "fa_icon": "fas fa-save"
                 },
+                "igenomes_base": {
+                    "type": "string",
+                    "format": "directory-path",
+                    "description": "Directory / URL base for iGenomes references.",
+                    "default": "s3://ngi-igenomes/igenomes/",
+                    "fa_icon": "fas fa-cloud-download-alt"
+                },
                 "igenomes_ignore": {
                     "type": "boolean",
                     "description": "Do not load the iGenomes reference config.",
@@ -499,11 +506,5 @@
         {
             "$ref": "#/$defs/generic_options"
         }
-    ],
-    "properties": {
-        "igenomes_base": {
-            "type": "string",
-            "default": "s3://ngi-igenomes/igenomes/"
-        }
-    }
+    ]
 }

subworkflows/local/mirna_quant.nf

@@ -24,6 +24,9 @@ include { EDGER_QC               } from '../../modules/local/edger_qc/main'
 include { BAM_STATS_MIRNA_MIRTOP } from '../../subworkflows/nf-core/bam_stats_mirna_mirtop/main'
 include { CSVTK_JOIN             } from '../../modules/nf-core/csvtk/join/main'
 
+include { PIVOT_LONGER           } from '../../modules/local/pivot/longer/main'
+include { PIVOT_WIDER            } from '../../modules/local/pivot/wider/main'
+
 workflow MIRNA_QUANT {
     take:
     ch_reference_mature  // channel: [ val(meta), fasta file]
@@ -94,12 +97,7 @@ workflow MIRNA_QUANT {
 
     ch_mirtop_logs = Channel.empty()
 
-    // nf-core/mirtop
-
-    ch_mirna_gtf_species = ch_mirna_gtf.map{ meta,gtf -> gtf }
-            .combine(ch_mirtrace_species)
-            .map{ gtf, species -> [ [id:species.toString()], gtf, species ] }
-            .collect()
+    ch_mirna_gtf_species = ch_mirna_gtf.map{ meta, gtf-> [ meta, gtf, meta.species ] }.collect()
 
     BAM_STATS_MIRNA_MIRTOP(BOWTIE_MAP_SEQCLUSTER.out.bam, FORMAT_HAIRPIN.out.formatted_fasta, ch_mirna_gtf_species )
 
@@ -110,10 +108,20 @@ workflow MIRNA_QUANT {
         .collect{it[1]}
         .map{it -> return [[id:"TSVs"], it]}
 
-    CSVTK_JOIN ( ch_tsvs )
-    ch_versions = ch_versions.mix(CSVTK_JOIN.out.versions)
+    PIVOT_LONGER( BAM_STATS_MIRNA_MIRTOP.out.counts )
+    ch_versions = ch_versions.mix(PIVOT_LONGER.out.versions)
+
+    ch_long_files = PIVOT_LONGER.out.csv
+        .map { meta, file -> file }
+        .collect()
+        .map { files ->
+            return [[id: "Long_Files"], files]
+        }
+
+    PIVOT_WIDER( ch_long_files )
+    ch_versions = ch_versions.mix(PIVOT_WIDER.out.versions)
 
-    DATATABLE_MERGE ( CSVTK_JOIN.out.csv )
+    DATATABLE_MERGE ( PIVOT_WIDER.out.csv )
     ch_versions = ch_versions.mix(DATATABLE_MERGE.out.versions)
 
     ch_reads_genome = BOWTIE_MAP_HAIRPIN.out.fastq

subworkflows/local/prepare_genome/main.nf

@@ -62,6 +62,15 @@ workflow PREPARE_GENOME {
     ch_mirtrace_species       = val_mirtrace_species          ? Channel.value(val_mirtrace_species) : Channel.empty()
     mirna_gtf_from_species    = val_mirtrace_species          ? (val_mirtrace_species == 'hsa' ? "https://raw.githubusercontent.com/nf-core/test-datasets/smrnaseq/reference/hsa.gff3" : "https://mirbase.org/download/CURRENT/genomes/${val_mirtrace_species}.gff3") : false
     ch_mirna_gtf              = val_mirna_gtf                 ? Channel.fromPath(val_mirna_gtf, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect() : ( mirna_gtf_from_species ? Channel.fromPath(mirna_gtf_from_species, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect() :  Channel.empty() )
+
+    // Add species of the gtf in the meta
+    ch_mirna_gtf = ch_mirna_gtf
+        .combine(ch_mirtrace_species.ifEmpty('unknown'))
+        .map { meta, gtf, species ->
+            def new_meta = meta.clone() + [species: species]
+            [new_meta, gtf]
+        }
+
     ch_mirna_adapters         = params.with_umi               ? [] : Channel.fromPath(val_fastp_known_mirna_adapters, checkIfExists: true).collect()
 
     ch_rrna                   = val_rrna                      ? Channel.fromPath(val_rrna, checkIfExists: true).map{ it -> [ [id:'rRNA'], it ] }.collect()                 : Channel.empty()
@@ -118,15 +127,29 @@ workflow PREPARE_GENOME {
 
     // Genome options
     if (!params.mirgenedb) {
-        ch_reference_mature  = params.mature  ? Channel.fromPath(params.mature, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect()  : { exit 1, "Mature miRNA fasta file not found: ${params.mature}" }
-        ch_reference_hairpin = params.hairpin ? Channel.fromPath(params.hairpin, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect() : { exit 1, "Hairpin miRNA fasta file not found: ${params.hairpin}" }
+        ch_reference_mature  = params.mature  ? Channel.fromPath(params.mature, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect()  : { exit 1, "Specify a Mature miRNA fasta file via '--mature'" }
+        ch_reference_hairpin = params.hairpin ? Channel.fromPath(params.hairpin, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect() : { exit 1, "Specify a Hairpin miRNA fasta file via '--hairpin'" }
     } else {
         if (!params.mirgenedb_species) {
             exit 1, "MirGeneDB species not set, please specify via the --mirgenedb_species parameter"
         }
-        ch_reference_mature  = params.mirgenedb_mature  ? Channel.fromPath(params.mirgenedb_mature, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect()  : { exit 1, "Mature miRNA fasta file not found via --mirgenedb_mature: ${params.mirgenedb_mature}" }
-        ch_reference_hairpin = params.mirgenedb_hairpin ? Channel.fromPath(params.mirgenedb_hairpin, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect() : { exit 1, "Hairpin miRNA fasta file not found via --mirgenedb_hairpin: ${params.mirgenedb_hairpin}" }
-        ch_mirna_gtf         = params.mirgenedb_gff     ? Channel.fromPath(params.mirgenedb_gff, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect()     : { exit 1, "MirGeneDB gff file not found via --mirgenedb_gff: ${params.mirgenedb_gff}"}
+        ch_reference_mature  = params.mirgenedb_mature  ? Channel.fromPath(params.mirgenedb_mature, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect()  : { exit 1, "Specify a mirgenedb Mature miRNA fasta file via '--mirgenedb_mature'" }
+        ch_reference_hairpin = params.mirgenedb_hairpin ? Channel.fromPath(params.mirgenedb_hairpin, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect() : { exit 1, "Specify a mirgenedb Hairpin miRNA fasta file via '--mirgenedb_hairpin'" }
+        ch_mirna_gtf         = params.mirgenedb_gff     ? Channel.fromPath(params.mirgenedb_gff, checkIfExists: true).map{ it -> [ [id:it.baseName], it ] }.collect()     : { exit 1, "Specify a MirGeneDB gff file via '--mirgenedb_gff'"}
+
+        // Create a channel for mirgenedb_species
+        ch_mirgenedb_species = Channel.value(params.mirgenedb_species)
+
+        // Add species of the gtf
+        // When mirgenedb workflow is not indicated, species defaults to val_mirtrace_species.
+        // If mirgenedb workflow parameters are indicated, the params.mirgenedb_species is used instead.
+        // If both mirgenedb workflow parameters and mirtrace_species (or mirna_gtf) are provided, params.mirgenedb_species is used as species value
+        ch_mirna_gtf = ch_mirna_gtf
+            .combine(ch_mirgenedb_species)
+            .map { meta, gtf, species ->
+                def new_meta = meta.clone() + [species: species]
+                [new_meta, gtf]
+            }
     }
 
     emit:

tests/test_contamination_tech_reps.nf.test

@@ -20,7 +20,7 @@ nextflow_pipeline {
             assertAll(
                 { assert workflow.success },
                 { assert snapshot(UTILS.removeNextflowVersion("$outputDir")).match("software_versions") },
-                { assert workflow.trace.succeeded().size() == 100 },
+                { assert workflow.trace.succeeded().size() == 103 },
 
                 { assert snapshot(
                     path("$outputDir/contaminant_filter/filter/Clone1_N1_trimmed.contamination_mqc.yaml").exists(), //TODO see if we can make these deterministic or why they are non-deterministic

tests/test_contamination_tech_reps.nf.test.snap

@@ -34,13 +34,13 @@
     },
     "software_versions": {
         "content": [
-            "{BLAT_CDNA={blat=36}, BLAT_NCRNA={blat=36}, BOWTIE2_ALIGN_CDNA={bowtie2=2.5.2, samtools=1.18, pigz=2.6}, BOWTIE2_ALIGN_NCRNA={bowtie2=2.5.2, samtools=1.18, pigz=2.6}, BOWTIE2_ALIGN_TRNA={bowtie2=2.5.2, samtools=1.18, pigz=2.6}, BOWTIE_MAP_HAIRPIN={bowtie=1.3.0, samtools=1.16.1}, BOWTIE_MAP_MATURE={bowtie=1.3.0, samtools=1.16.1}, BOWTIE_MAP_SEQCLUSTER={bowtie=1.3.0, samtools=1.16.1}, CAT_FASTQ={cat=8.3}, CSVTK_JOIN={csvtk=0.30.0}, DATATABLE_MERGE={r-base=3.6.2}, FASTP={fastp=0.23.4}, FILTER_STATS={BusyBox=1.32.1}, FORMAT_HAIRPIN={fastx_toolkit=0.0.14}, FORMAT_MATURE={fastx_toolkit=0.0.14}, GAWK_CDNA={gawk=5.3.0}, GAWK_NCRNA={gawk=5.3.0}, INDEX_CDNA={bowtie2=2.5.2}, INDEX_HAIRPIN={bowtie=1.3.0}, INDEX_MATURE={bowtie=1.3.0}, INDEX_NCRNA={bowtie2=2.5.2}, INDEX_TRNA={bowtie2=2.5.2}, MIRTOP_COUNTS={mirtop=0.4.28}, MIRTOP_EXPORT={mirtop=0.4.28}, MIRTOP_GFF={mirtop=0.4.28}, MIRTOP_STATS={mirtop=0.4.28}, MIRTRACE_QC={mirtrace=1.0.1}, PARSE_HAIRPIN={seqkit=2.6.1}, PARSE_MATURE={seqkit=2.6.1}, SAMTOOLS_FLAGSTAT={samtools=1.21}, SAMTOOLS_IDXSTATS={samtools=1.21}, SAMTOOLS_INDEX={samtools=1.21}, SAMTOOLS_SORT={samtools=1.21}, SAMTOOLS_STATS={samtools=1.21}, SEQCLUSTER_COLLAPSE={seqcluster=1.2.9}, SEQKIT_GREP_CDNA={seqkit=2.8.0}, SEQKIT_GREP_NCRNA={seqkit=2.8.0}, STATS_GAWK_CDNA={gawk=5.3.0}, STATS_GAWK_NCRNA={gawk=5.3.0}, STATS_GAWK_TRNA={gawk=5.3.0}, Workflow={nf-core/smrnaseq=v2.4.0}}"
+            "{BLAT_CDNA={blat=36}, BLAT_NCRNA={blat=36}, BOWTIE2_ALIGN_CDNA={bowtie2=2.5.2, samtools=1.18, pigz=2.6}, BOWTIE2_ALIGN_NCRNA={bowtie2=2.5.2, samtools=1.18, pigz=2.6}, BOWTIE2_ALIGN_TRNA={bowtie2=2.5.2, samtools=1.18, pigz=2.6}, BOWTIE_MAP_HAIRPIN={bowtie=1.3.0, samtools=1.16.1}, BOWTIE_MAP_MATURE={bowtie=1.3.0, samtools=1.16.1}, BOWTIE_MAP_SEQCLUSTER={bowtie=1.3.0, samtools=1.16.1}, CAT_FASTQ={cat=8.3}, DATATABLE_MERGE={r-base=3.6.2}, FASTP={fastp=0.23.4}, FILTER_STATS={BusyBox=1.32.1}, FORMAT_HAIRPIN={fastx_toolkit=0.0.14}, FORMAT_MATURE={fastx_toolkit=0.0.14}, GAWK_CDNA={gawk=5.3.0}, GAWK_NCRNA={gawk=5.3.0}, INDEX_CDNA={bowtie2=2.5.2}, INDEX_HAIRPIN={bowtie=1.3.0}, INDEX_MATURE={bowtie=1.3.0}, INDEX_NCRNA={bowtie2=2.5.2}, INDEX_TRNA={bowtie2=2.5.2}, MIRTOP_COUNTS={mirtop=0.4.28}, MIRTOP_EXPORT={mirtop=0.4.28}, MIRTOP_GFF={mirtop=0.4.28}, MIRTOP_STATS={mirtop=0.4.28}, MIRTRACE_QC={mirtrace=1.0.1}, PARSE_HAIRPIN={seqkit=2.6.1}, PARSE_MATURE={seqkit=2.6.1}, PIVOT_LONGER={r-base=4.4.2, tidyr=1.3.1, optparse=1.7.5}, PIVOT_WIDER={r-base=4.4.2, tidyr=1.3.1, dplyr=1.1.4, optparse=1.7.5, vroom=1.6.5}, SAMTOOLS_FLAGSTAT={samtools=1.21}, SAMTOOLS_IDXSTATS={samtools=1.21}, SAMTOOLS_INDEX={samtools=1.21}, SAMTOOLS_SORT={samtools=1.21}, SAMTOOLS_STATS={samtools=1.21}, SEQCLUSTER_COLLAPSE={seqcluster=1.2.9}, SEQKIT_GREP_CDNA={seqkit=2.8.0}, SEQKIT_GREP_NCRNA={seqkit=2.8.0}, STATS_GAWK_CDNA={gawk=5.3.0}, STATS_GAWK_NCRNA={gawk=5.3.0}, STATS_GAWK_TRNA={gawk=5.3.0}, Workflow={nf-core/smrnaseq=v2.4.0}}"
         ],
         "meta": {
             "nf-test": "0.9.0",
-            "nextflow": "24.04.4"
+            "nextflow": "24.10.2"
         },
-        "timestamp": "2024-10-08T23:16:26.853242481"
+        "timestamp": "2024-12-10T00:29:32.052341"
     },
     "mirna_quant_bam": {
         "content": [
@@ -65,9 +65,9 @@
         ],
         "meta": {
             "nf-test": "0.9.0",
-            "nextflow": "24.04.4"
+            "nextflow": "24.10.2"
         },
-        "timestamp": "2024-10-01T20:06:04.974546479"
+        "timestamp": "2024-12-10T00:29:32.116301175"
     },
     "mirna_quant_edger_qc": {
         "content": [
@@ -90,9 +90,9 @@
         ],
         "meta": {
             "nf-test": "0.9.0",
-            "nextflow": "24.04.4"
+            "nextflow": "24.10.2"
         },
-        "timestamp": "2024-10-01T20:06:05.025175321"
+        "timestamp": "2024-12-10T00:29:32.164075991"
     },
     "contaminant_filter_filter": {
         "content": [
@@ -113,8 +113,8 @@
         ],
         "meta": {
             "nf-test": "0.9.0",
-            "nextflow": "24.04.4"
+            "nextflow": "24.10.2"
         },
-        "timestamp": "2024-10-01T20:06:05.070939602"
+        "timestamp": "2024-12-10T00:29:32.208602197"
     }
-}
+}

tests/test_mirgenedb.nf.test

@@ -19,7 +19,7 @@ nextflow_pipeline {
             assertAll(
                 { assert workflow.success },
                 { assert snapshot(UTILS.removeNextflowVersion("$outputDir")).match("software_versions") },
-                { assert workflow.trace.succeeded().size() == 90 },
+                { assert workflow.trace.succeeded().size() == 107 },
                 { assert workflow.trace.failed().size() == 1 },
 
                 { assert snapshot(