nf-core/ampliseq
Amplicon sequencing analysis workflow using DADA2 and QIIME2
Version history
nf-core/ampliseq version 2.15.0 - 2025-10-06
Summary of changes
Many issues were fixed and more control was given over QC and taxonomic classification.
Detailed changes
Added
- #909,#914 - Parameter
--truncqallows read truncation by quality score, where each read is truncated at the first instance of a quality score less than or equal totruncq(default value is 2). This quality score-based truncation occurs before read length truncation. If--trunc_qminand--trunc_rminare used to automatically calculate the values fortrunclenfandtrunclenrused for read length truncation, those calculations are based on the read metrics before quality score-based truncation (usingtruncq) is performed.truncqis passed directly astruncQinto DADA2’s filterAndTrim method. - #913 - Parameter
--dada_min_bootallows the user to specify a minimum bootstrap confidence (out of 100 trials) for assigning a taxonomic level with DADA2’s assignTaxonomy method.
Changed
- #903 - Corrects a broken link for midori2-co1 database.
Fixed
- #893,#896,#905,#912 - Template update for nf-core/tools version 3.3.2 and reorganization of nf-test
- #897 - Allow multiple region analysis with one sample.
- #899 - With
--dada_ref_tax_customall chunks defined by--dada_assign_chunksizeare now taxonomically annotated. - #904 - Update
untarmodule. - #906 - Fix
--mergepairs_strategy "consensus"causing an error in dada2_stats when analysing a run composed by only one sample.
Removed
- #894 - Remove BOLD databases
plantae-bold-its1&plantae-bold-trnL, for plant samples with markers ITS1 and trnL, versions 20240510. Database included low quality sequences.
nf-core/ampliseq version 2.14.0 - 2025-05-23
Summary of changes
- Multi-region analysis is now easier with custom databases, parameters were changed or added
- Updated
sbdi-gtdb,gtdb,PR2databases and addedBOLDplant databases for classification with DADA2 - Follow nextflow’s strict syntax (only impacts development)
Detailed changes
Added
- #879 - Add
sbdi-gtdb=R10-RS226-1as parameter to--dada_ref_taxonomyand make this version the default forsbdi-gtdb. - #882 - Add
gtdb=R10-RS226as parameter to--dada_ref_taxonomyand make this version the default forgtdb. - #883 - Add BOLD databases
plantae-bold-its1&plantae-bold-trnL, for plant samples with markers ITS1 and trnL, versions 20240510, for DADA2 as parameter to--dada2_ref_taxonomy. - #885 - Add
pr2=5.1.0as parameter to--dada_ref_taxonomyand make this version the default forpr2.
Changed
- #871 - Multi-region analysis is now easier with custom databases. No alignment file of the reference sequences are needed.
| Parameter | Description |
|---|---|
| sidle_ref_tax_custom | Changed: Previously three files, now one file: path to reference taxonomy strings |
| sidle_ref_seq_custom | New: Path to reference taxonomy sequences in fasta format |
| sidle_ref_aln_custom | New: Path to multiple sequence alignment of reference taxonomy sequences in fasta format |
| sidle_ref_degenerates | New: Only effective with --sidle_ref_tax_custom, filter reference sequences, default: 5 |
| sidle_ref_cleaning | New: Arguments regarding ad-hoc cleaning, with --sidle_ref_tax_custom default is ‘—p-database silva’ |
Fixed
- #872 - Follow nextflow’s strict syntax
- #876 - Pulled the updated vsearch/clusters module from nf-core to fix a bug where a wildcard expansion trigging an arguments-too-long error.
- #878 - Downgraded nf-schema from 2.3.0 to 2.2.0 due to incompatibilities with nextflow stable versions 25.05 and newer
- #881 - Template update for nf-core/tools version 3.2.1
- #887 - Updated documentation to include all updates
nf-core/ampliseq version 2.13.0 - 2025-04-04
Summary of changes
- Zenodo changed from zenodo.org/record to zenodo.org/records and some taxonomic databases could not be downloaded any more
- Add most recent database versions for SILVA, Greengenes2, and SBDI-GTDB for DADA2 or QIIME2
- Export of a TreeSummarizedExperiment R object
- Taxonomic assignment using DADA2 is done in chunks (default 10,000 ASVs) to reduce memory and runtime
- Consensus strategy to merge or concatenate paired end reads based on overlap is introduced
- A bunch of bug fixes
- nf-core template update from 3.0.2 to 3.2.0
Detailed changes
Added
- #798 - Added SILVA version 138.2 of DADA2 taxonomy database:
silva=138.2orsilvaas parameter to--dada2_ref_taxonomy - #801 - Parameter
--quality_typeallows specifying the type of quality scores in raw read data, by defaultAuto(i.e. default behavior did not change) - #804 - Added version 10 of Unite as parameter for
--dada_ref_taxonomy(issue #768) - #803,#856 - New parameters introduced related to
--mergepairs_strategy. These parameters would only be effective if--mergepairs_strategy consensusis set. - #807 - Export of TreeSummarizedExperiment R object by default, can be omitted with
--skip_tse, also added ability to skip phyloseq R object generation with--skip_phyloseq - #833 - Add paths to updated version of SBDI-GTDB database.
- #849 - Added Greengenes version 2024.09 of DADA2 taxonomy database:
greengenes2=2024.09orgreengenes2as parameter to--dada2_ref_taxonomy - #853 - Added Greengenes2 version 2024.09 for taxonomic classification with QIIME2:
greengenes2=2024.09orgreengenes2as parameter to--qiime_ref_taxonomy - #864 - Taxonomic assignment using DADA2 is done in chunks of 10,000 ASVs by default, the chunk size can be adjusted with
--dada_assign_chunksize
| New parameter | Description | Default Value |
|---|---|---|
| quality_type | Set the input data quality type for DADA2 | ”Auto” |
| mergepairs_strategy | Choose the method of pair end merging strategy | ”merge” |
| mergepairs_consensus_match | The score assigned for each matching base pair during sequence alignment. | 1 |
| mergepairs_consensus_mismatch | The penalty score assigned for each mismatched base pair during sequence alignment. | -2 |
| mergepairs_consensus_gap | The penalty score assigned for each gap introduced during sequence alignment. | -4 |
| mergepairs_consensus_minoverlap | The minimum number of overlapping base pairs required to merge forward and reverse reads. | 12 |
| mergepairs_consensus_maxmismatch | The maximum number of mismatches allowed within the overlapping region for merging reads. | 0 |
| mergepairs_consensus_percentile_cutoff | The percentile cutoff determining the minimum observed overlap in the dataset. | 0.001 |
| skip_tse | Skip generation of a TreeSummarizedExperiment R object | false |
| skip_phyloseq | Skip generation of a Phyloseq R object | false |
Changed
- #803 - Changed DADA2_DENOISING :
--concatenate_readsrenaming to--mergepairs_strategy; support new method named “consensus” by setting--mergepairs_strategy consensus; changed options of--mergepairs_strategyfrom TRUE/FALSE (boolean) to [“merge”, “concatenate”, “consensus”]. - #818 - Provide users the ability to not bump stack size in vsearch clustering.
Fixed
- #837 - Replaced zenodo.org/record to zenodo.org/records in ref_databases.config to reduce DB download failures
- #800 - Fixed SH files for UNITE9.0, they were missing some entries due to a bug caused by API update in PlutoF
- #808 - Add missing library declaration in R script.
- #847 - update the
usage.mdabout sample id format (no dashes allowed) - #832,#848 - Add pattern to input fields
- #845 - Template update for nf-core/tools version 3.2.0
- #846 - Validation of sample sheet doesnt allow intergers for sampleID any more
- #854,#858 - QIIME2 and Sidle are now also available with conda/mamba
- #862 - ITSx serves now as filter of the abundance table
Dependencies
| software | previously | now |
|---|---|---|
| QIIME2 | 2023.7 | 2024.10 |
| MultiQC | 1.25.1 | 1.27 |
Removed
- #860 - Removed unused test profile
test_its_dada_taxonomy.
nf-core/ampliseq version 2.12.0 - 2024-11-15
Summary of changes
- updated some reference databases
- several bug fixes
- updated to nf-core template 3.0.2
- dependency updates
Detailed changes
Added
- #786 - Added version 20240904 of Kraken2 standard taxonomy database:
standard=20240904orstandardas parameter to--kraken2_ref_taxonomy
Changed
- #785 - replaced https://data.ace.uq.edu.au/public/gtdb/data by https://data.gtdb.ecogenomic.org for GTDB versions via
--dada_ref_taxonomy
Fixed
- #785,#795 - Template update for nf-core/tools version 3.0.2
- #786 - Creation of the summary report was failing with
--dada_ref_taxonomy silva=132 - #792 - Remove white spaces in process tags to avoid errors on some systems
Dependencies
| software | previously | now |
|---|---|---|
| MultiQC | 1.21 | 1.25.1 |
| Kraken2 | 2.1.2 | 2.1.3 |
| Picrust2 | 2.5.2 | 2.5.3 |
nf-core/ampliseq version 2.11.0 - 2024-08-06
Summary of changes
This update is mostly about reference taxonomy databases:
- added reference databases
sbdi-gtdb=R09-RS220-1for DADA2 - added reference databases
unite-fungi=10.0andunite-alleuk=10.0for SINTAX - disabled UNITE reference databases for QIIME2
- fixed broken urls for UNITE databases
Detailed changes
Added
- #765 - Added version R09-RS220 of curated GTDB 16S taxonomy:
sbdi-gtdb=R09-RS220-1orsbdi-gtdbas parameter to--dada_ref_taxonomy - #766 - Added version 10 of Unite as parameter for
--sintax_ref_taxonomy:unite-fungi=10.0andunite-alleuk=10.0
Changed
- #762 - Improved output documentation section “Optional ASV filtering” and parameter documentation
- #766 - Modified warning filenames from
QIIME2_ANCOMto avoid collisions - #766,#769 - Disabled Unite databases from the
--qiime_ref_taxonomybecause of divergent results compared to the other classifiers
Fixed
nf-core/ampliseq version 2.10.0 - 2024-06-27
Summary of changes
- added ANCOM-BC for differential abundance analysis
- updated reference databases for
sbdi-gtdb>dbwithR08-RS214-1andR09-RS220, respectively - template update 2.13.1 to 2.14.1
- several bug fixes
Detailed changes
Added
- #751 - Added version R08-RS214 of curated GTDB 16S taxonomy:
sbdi-gtdb=R08-RS214-1orsbdi-gtdbas parameter to--dada_ref_taxonomy - #752 - Added version R09-RS220 of GTDB 16S taxonomy:
gtdb=R09-RS220orgtdbas parameter to--dada_ref_taxonomy - #753, #756, #757 - ANCOM-BC via QIIME2 can be used with
--ancombc,--ancombc_formula, and--ancombc_formula_reflvl, plotting can be modified with thresholds--ancombc_effect_sizeand--ancombc_significance
Changed
- #749 - Create barplot also when no metadata is given
- #753 - ANCOM via QIIME2 is not run anymore by default but on request whith
--ancom, therefore--skip_ancomwas removed
Fixed
- #747 - Template update for nf-core/tools version 2.14.1
- #748 - Updating misleading error message and documentation
- #750 - Numbers in
overall_summary.tsvwere fixed (sometimes misleading in 2.9.0 for columns “denoised[F/R]”, “merged”, and “nochim”)
Removed
- #753 -
--skip_ancomwas removed
nf-core/ampliseq version 2.9.0 - 2024-04-03
Summary of changes
- added multiple regions analysis (including 5R / SMURF / q2-sidle), integrating pipesidle and its software container with 1 & 2.
- template update 2.11 to 2.13.1
- several bug fixes
Detailed changes
Added
- #700 - Optional
--save_intermediatesto publish QIIME2 data objects (.qza) and visualisation objects (.qzv) - #702,#723,#728,#729 - Add multiple regions analysis (including 5R / SMURF / q2-sidle)
Changed
- #719 - Versions of all (instead of selected) processes are now exported to
pipeline_info/software_versions.yml
Fixed
- #697,#699,#713 - Template update for nf-core/tools version 2.13.1
- #711 - From r207 and onwards Archaea sequences were omitted when parsing GTDB databases. (This did not affect
sbdi-gtdbdatabases, onlygtdb.) - #715 - Fix filtering vsearch clusters for high number of clusters
- #717 - Fix edge case for sorting file names by using radix method
- #718 - Require a minimum sequence length of 50bp for taxonomic classifcation after using ITSx
- #721 - Fix error
unknown recognition error type: groovyjarjarantlr4.v4.runtime.LexerNoViableAltExceptioncaused by a missing\in nf-core modulepigz/uncompress(which had no consequences but was confusing) - #722 - When barrnap detects several genes select the lowest e-value
- #726 - Add fallback to
download_pipeline.ymlbecause the pipeline does not support stub runs (#2846)
Dependencies
- #720 - Updated nf-core modules, DADA2, and Phyloseq
| software | previously | now |
|---|---|---|
| cutadapt | 3.4 | 4.6 |
| DADA2 | 1.28.0 | 1.30.0 |
| Phyloseq | 1.44.0 | 1.46.0 |
Removed
- #710 - Removed Phyloref from DADA2 reference option because it’s part of PR2 5.0.0
nf-core/ampliseq version 2.8.0 - 2024-01-16
Summary of changes
- Three additional reference databases available:
--qiime_ref_taxonomy greengenes2(prokaryotic 16S),--dada_ref_taxonomy phytoref(eukaryotic plastid 16S rRNA),--dada_ref_taxonomy zehr-nifh(Nitrogenase iron protein NifH) - New parameters
--qiime_ref_tax_custom,--skip_qiime_downstream,--dada_taxonomy_rc, and--dada_addspecies_allowmultiple - Bug fixes for using phyloseq, VSEARCH, and error messages
- Updated documentation & reporting
- Pipeline template update
Detailed changes
Added
- #666 - Added Greengenes2 database, version 2022.10, support for QIIME2 taxonomic classification.
- #667,#691 - Added
--qiime_ref_tax_customto permit custom reference database for QIIME2 taxonomic classification - #674 - Add PhytoRef database for DADA2 taxonomy assignment using
--dada_ref_taxonomy phytoref - #675 - Add the Zehr lab nifH database for DADA2 taxonomy assignment using
--dada_ref_taxonomy zehr-nifh - #681 - For DADA2, with
--dada_addspecies_allowmultiplemultiple exact species matches are reported and with--dada_taxonomy_rcreverse-complement matches are also considered in taxonomic classification
Changed
- #677 - Added cut_its information to SDBI export
Fixed
- #672,#688,#691 - Updated documentation
- #676 - Phyloseq sometimes only produced one of multiple output files
- #679 - Prevent masking low complexity regions by VSEARCH with lower case letters
- #680,#673 - Improved pipeline summary report & error messages
- #683 - Template update for nf-core/tools version 2.11
- #687 - Correct conda package for ASV SSU filtering
nf-core/ampliseq version 2.7.1 - 2023-11-13
Summary of changes
Downstream analysis with QIIME2 fails on some system with the most recent version of nextflow (23.10.0), this release fixes the issue.
Detailed changes
Changed
- #657 - Improved parameter descriptions and sequence
Fixed
nf-core/ampliseq version 2.7.0 - 2023-10-20
Summary of changes
- Added Pipeline RMarkdown summary report
- Added optional ASV to OTU post-clustering with VSEARCH
- Added Kraken2 as another optional taxonomic classification tool
- Splitted
--inputinto 3 parameters,--inputnow points to a sampleheet, alternatives are--input_fasta&--input_folderto accommodate nf-validation’sfromSamplesheet - Update to nf-core template to 2.10
Detailed changes
Added
- #558,#619,#625,#632,#644 - Pipeline summary report
- #615 - Phyloseq R object creation
- #622 - ASV post-clustering with Vsearch
- #637 - Taxonomic classification with Kraken2, parameter
--kraken2_ref_taxonomy,--kraken2_ref_tax_custom,--kraken2_assign_taxlevels,--kraken2_confidence - #639 - GTDB release 214.1 for taxonomic classification with DADA2, using
--dada_ref_taxonomy gtdbor--dada_ref_taxonomy gtdb=R08-RS214 - #641 - Continue analysis even when individual files fail the filtering threshold, added parameter
--ignore_failed_filtering
Changed
- #616 - When using a sample sheet with
--inputcontaining forward and reverse reads, specifying--single_endwill only extract forward reads and treat the data as single ended instead of extracting forward and reverse reads. - #616 -
--inputwas split into three params: (1)--inputfor samplesheet, (2)--input_fastafor ASV/OTU fasta input, (3)--input_folderdirect FASTQ input
| Param updated | Param old | Accepts |
|---|---|---|
| input | input | samplesheet, .tsv/.csv/.yml/.yaml |
| input_fasta | input | ASV/OTU sequences, .fasta |
| input_folder | input | Folder containing compressed fastq files |
- #639 -
--dada_ref_taxonomy gtdbpoints towards GTDB release 214.1 instead of GTDB release 207 for taxonomic classification with DADA2 - #645 - Updated documentation, including workflow figure
Fixed
- #605 - Make
--sbdiexportcompatible with PR2 version 5.0.0 - #614,#620,#642 - Template update for nf-core/tools version 2.10
- #617 - Fix database compatibility check for
--sbdiexport - #628 - Fix edge case for sample sheet input when using specific combinations of sampleID and forwardReads or reverseReads that will forward one file too much to cutadapt
- #630 - ASV rRNA (barrnap), length, and codon filter now work with ASV fasta file input
- #633 - UNIFRAC in QIIME2_DIVERSITY_CORE is now prevented from using a GPU to avoid errors
- #643 - Fix using
--skip_dada_addspecieswithout--dada_ref_tax_custom_spwhich was broken in 2.6.0 & 2.6.1 - #647 - Update of credits
Dependencies
- #646 - Updated dependencies, see below:
| software | previously | now |
|---|---|---|
| FASTQC | 0.11.9 | 0.12.1 |
| DADA2 | 1.22.0 | 1.28.0 |
| PICRUSt2 | 2.5.0 | 2.5.2 |
| QIIME2 | 2022.11 | 2023.7 |
nf-core/ampliseq version 2.6.1 - 2023-06-28
Summary of changes
This release fixes a bug with container registry.
Detailed changes
Fixed
- #603 - Fix all containers registry
nf-core/ampliseq version 2.6.0 - 2023-06-27
Summary of changes
- Updated taxonomic databases (UNITE 9.0, PR2 5.0.0)
- Added more choice for taxonomic classification: phylogenetic placement, SINTAX
- Codon usage based filtering allows for additional QC of coding regions
- Add NF-TEST pipeline end-to-end tests for existing CI tests
- nf-core template update from 2.7.2 to 2.8
Detailed changes
Added
- #564,#567,#582 - Added phylogenetic placement
- #577 - Added SINTAX for taxonomic classification
- #575, #586 - Added filtering step for stop codons for ASVs that are of coding regions.
- #580 - Add NF-TEST pipeline end-to-end tests for existing CI tests
- #591 - New version of the Unite taxonomy databases: 9.0
- #596 - New version of the PR2 taxonomy database: 5.0.0, only available with DADA2 (
--dada_ref_taxonomy) - #597 - Samples with less reads than specified with
--min_read_counts(default: 1) stop the pipeline, previously the threshold was 1KB in size.
Changed
- #563 - Renamed DADA2 taxonomic classification files to include the chosen reference taxonomy abbreviation.
- #567 - Renamed
--dada_tax_agglom_minand--qiime_tax_agglom_minto--tax_agglom_minand--dada_tax_agglom_maxand--qiime_tax_agglom_maxto--tax_agglom_max - #580 - GitHub Actions CI - pull_request to
devtests with NXF_VERlatest-everything& pull_request tomastertests with NXF_VER22.10.1&latest-everything - #598 - Updated Workflow figure with SINTAX and phylogenetic placement
- #599 - For exact species assignment (DADA2’s addSpecies) PR2 taxonomy database (e.g.
--dada_ref_taxonomy pr2) now excludes any taxa that end with ” sp.”.
Fixed
- #553 - Handle empty barrnap results files
- #554 - Accept taxonomy strings that contain
#,' - #569 - Make header of overall_summary.tsv consistent between input data types
- #573 - Avoid parser error for single-end data when an empty read file is detected
- #578 - Template update for nf-core/tools version 2.8, including changing
System.exit(1)toNextflow.error() - #594 - Update metadata documentation
- #595 - Closing gaps in rarefaction depth for diversity calculations (
mindepthin QIIME2_DIVERSITY_CORE)
nf-core/ampliseq version 2.5.0 - 2023-03-02
Summary of changes
- Added COIDB DADA2 taxonomic database
- Small adjustments to taxonomic classification, ASV filtering by barrnap, fasta input
- sbdi-export works for more taxonomic databases now
- nf-core template update from 2.6 to 2.7.2
Detailed changes
Added
- #518,#534 - Add COIDB DADA2 reference taxonomy database
- #521 - Export svg in addition to pdf files for quality plots from DADA2
- #538 - Parameter
--diversity_rarefaction_depthcontrols the minimum rarefaction depth for diversity analysis, this allows increasing the rarefaction depth at the cost of excluding low count samples. Parameter--ancom_sample_min_countsets the minimum sample counts to retain a sample for ANCOM analysis.
Changed
- #537 - Update output generated with option sbdi-export
- #541 - Remove adjustments of taxonomic levels for RDP & SILVA & GTDB & UNITE database for DADA2 taxonomic classification, reduced default of
--dada_tax_agglom_maxfrom 7 to 6 - #548 -
--filter_ssuaccepted any barrnap hit to a kingdom (domain) (any occurence in resulting gff) to choose an ASV, now only ASVs with the kingdom (domain) that has the lowest evalue are accepted.
Fixed
- #513 - Template update for nf-core/tools version 2.7.2
- #519 - Adding the pipeline reference to the MultiQC report
- #520,#530 - Fix conda packages
- #531,#546 - Update documentation
- #535 - Make sure barrnap runs with fasta input
- #544 - Adding module to fix header in fasta input if needed
Dependencies
- #528 - Updated QIIME2
| Tool | Previous version | New version |
|---|---|---|
| QIIME2 | 2022.8 | 2022.11 |
Removed
- #513 - Removed parameter
--enable_conda.
nf-core/ampliseq version 2.4.1 - 2022-12-07
Summary of changes
Most importantly this release fixes a bug in 2.4.0 that prevented the use of the SILVA v138 database with QIIME2.
Detailed changes
Added
- #494 -
--metadata_category_barplotaccepts a comma separated list of metadata categories and plots for each barplots with average relative abundance.
Changed
- #492 -
--qiime_adonis_formulaaccepts a comma separated list of formulas.
Fixed
- #486 - Fixed typo in error message stating
--skip_classiferinstead of--classifier. - #487,#488 - Update stale links in usage documentation.
- #489 - Reduce linting warnings for nf-core tools version 2.5.1.
- #491 - Make output from —addSH match UNITE format by replacing spaces with underscores.
- #495 - Template update for nf-core/tools version 2.6
- #501 - Check for empty fields in samplesheet column “run” and raise an appropriate error.
- #503 - Changed environment for formatting databases.
- #504 - Fixed warnings with nextflow 22.10 (and later) about processes that are defined more than once.
Dependencies
| Tool | Previous version | New version |
|---|---|---|
| QIIME2 | 2021.8 | 2022.8 |
nf-core/ampliseq version 2.4.0 - 2022-09-07
Summary of changes
- Custom DADA2 reference taxonomy databases are now supported.
- Illumina NovaSeq data can be correctly processed.
- ASV length filter is available for stricter QC.
- Species classification is now reported in columns “Species” using DADA2’s assignTaxonomy function or “Species_exact” using DADA2’s addSpecies function, the latter only assigns exact sequence matches.
Detailed changes
Added
- #456 - An optional ASV length filter can be activated using
--min_len_asv <int>and/or--max_len_asv <int>. - #458 - Samplesheet, ASV fasta file, and/or metadata sheet is now exported into
<results>/input/ - #459 - MIDORI2 CO1 database with keys
midori2-co1=gb250andmidori2-co1for--dada_ref_taxonomy - #460 - Taxonomic ranks for DADA2 taxonomic classification can be now adjusted using
--dada_assign_taxlevels <comma separated string>. - #461 - A custom DADA2 reference taxonomy database can now be used with
--dada_ref_tax_customand--dada_ref_tax_custom_sp, typically accompanied by--dada_assign_taxlevels. - #446,#467 - Binned quality scores from Illumina NovaSeq data can be now corrected with
--illumina_novaseq. - #477 - QC plots of DADA2’s plotQualityProfile are now also produced after preprocessing.
- #478 - Added GTDB R07-RS207 DADA2 taxonomy reference databases
Changed
- #444,#457,#463,#465,#466,#469 - Updated the documentation.
- #445 - The minimum number of total bases to use for error rate learning by default is 1e8 (DADA2, learnErrors, nbases). Previously, samples were read in the provided order until enough reads were obtained (DADA2, learnErrors, randomize=FALSE). Now, samples are picked at random from those provided (DADA2, learnError, randomize=TRUE) and a seed is set.
- #453 - Export a few more basic QIIME2 fragments (zipped files) that can be easily imported into the correct QIIME2 version for custom analysis.
- #464 - Reported taxonomic classifications on species level based on DADA2’s assignTaxonomy (approximations) is now listed in column “Species” while exact matches based on DADA2’s addSpecies are now reported in column “Species_exact”.
Fixed
- #448 - Updated SBDI export scripts to include Unite species hypothesis information if available.
- #451 - Pairwise statistics will be now performed on a subset of metadata columns specified with
--metadata_categoryinstead of ignoring that setting. - #451 - Replace busybox with Ubuntu base image for GCP support.
- #455 - Stop with descriptive error when only one of
--trunclenfand--trunclenris given, earlier it was silently ignored. - #474 - Template update for nf-core/tools version 2.5.1
- #475 - Report software versions for DADA2_TAXONOMY
Dependencies
- #479 - Updated software
| Tool | Previous version | New version |
|---|---|---|
| PICRUSt2 | 2.4.2 | 2.5.0 |
| MultiQC | 1.12 | 1.13a |
nf-core/ampliseq version 2.3.2 - 2022-05-27
Summary of changes
This release fixes bugs and adds more options to manipulate cutadapt parameters.
Detailed changes
Added
- #429 -
--cutadapt_min_overlapsets cutadapt’s global minimum overlap (-O) and--cutadapt_max_error_ratesets cutadapt’s global maximum error rate (-e) for trimming primer sequences. - #431 -
--skip_dada_qualityallows to skip quality check with DADA2. This is only allowed when--trunclenfand--trunclenrare set.
Changed
- #432 - The number of records to sample from a fastq file was decreased from 5e+06 to 5e+04 for plotQualityProfile (DADA2_QUALITY), therefore a smaller subset of reads is sampled for determining
--trunlenfand--trunclenr. This should make the process more robust also from larger data sets.
Fixed
nf-core/ampliseq version 2.3.1 - 2022-04-05
Summary of changes
This release fixes a bug.
Detailed changes
Fixed
- #415 - ADONIS test was running by default in version 2.2.0 and 2.3.0 on metadata columns that did not always have the required format, now specifying
--qiime_adonis_formulais required to run that step.
nf-core/ampliseq version 2.3.0 - 2022-04-04
Summary of changes
- Added Barrnap that annotates whether ASV sequences are SSUs and filter on demand
- Primer sequences are now not required any more and cutadapt can be skipped
- Updated to nf-core template 2.3.2
- Fixed some bugs
Detailed changes
Added
- #385 -
--skip_cutadaptallows to skip primer trimmimg. - #390, #408 - Add values to option
--cut_its, specifying which ITS part to use. Also, add option--its_partial <x>to allow partial ITS sequences longer than given cutoff. - #395 -
--seedspecifies the random seed. - #396 - Barrnap annotates ASV sequences for SSU’s, it can be skipped with
--skip_barrnap.--filter_ssutakes a comma separated list of “bac,arc,mito,euk” and enables SSU filtering depending on Barrnap (default: off). - #397 - Complement README.md with links to the nf-core bytesize 25 (nf-core/ampliseq).
Changed
- #385 -
--FW_primerand--RV_primerare not obligatory any more, however primer sequences are still required with cutadapt (i.e. without--skip_cutadapt),--qiime_ref_taxonomy, and--cut_dada_ref_taxonomy(cuts reference with primer sequences).
Fixed
- #384 - For QIIME2 beta diversity, make directory before execution.
- #394 - Prevent simultaneous usage of
--qiime_ref_taxonomyand--classifier. - #402, #410 - Template update for nf-core/tools version 2.3.2
- #403 - Limit number of files for DADA2_QUALITY (plotQualityProfile) by read numbers.
- #405 - Fix omitting taxonomic filtering (QIIME2_FILTERTAXA).
- #404 - Rephrase error message for empty input files and empty files after trimming with cutadapt.
- #407 - Fix reformatting script for Unite.
Dependencies
| Tool | Previous version | New version |
|---|---|---|
| MultiQC | 1.11 | 1.12 |
nf-core/ampliseq version 2.2.0 - 2022-01-31
Summary of changes
- Added features (e.g. SBDI-GTDB db update, Adonis test)
- Added convenience updates (e.g. properly reporting and on demand ignoring low read samples)
- Updated to nf-core template 2.2
- Fixed some bugs
Detailed changes
Added
- #352, #372 -
--skip_dada_addspeciesallows to skip species level classification to reduce memory requirements, incompatible with--sbdiexportthat expect species annotation. - #364 - Adonis in QIIME2 for testing feature importance in beta diversity distances,
--qiime_adonis_formulacan be set to provide a custom formula. - #366 - New version of the SBDI-GTDB taxonomy database: v. 3. (Fixes problem with
Reverse_added to some domain strings.)
Changed
- #354 - Input files and files after primer trimming with cutadapt are required to be >1KB (i.e. not empty) and either the pipeline will stop if at least one sample file fails or the failing samples will be ignored when using
--ignore_empty_input_filesor--ignore_failed_trimming, respectively. - #376 - Forbid sampleIDs starting with a number when also
--metadatais used, because such strings are unintentionally modified and the metadata will not match any more.
Fixed
- #377- An error message will occur when
--sbdiexportis used with--skip_taxonomyor--skip_dada_addspecies - #375- Updated documentation regarding not using curly brackets in
--extensionwith--single_end - #362- Template update for nf-core/tools version 2.2, now requires nextflow version
>= 21.10.3 - #374- Cutadapt results can be now also viwed in the MultiQC report
Dependencies
| Tool | Previous version | New version |
|---|---|---|
| Cutadapt | 3.2 | 3.4 |
| DADA2 | 1.20.0 | 1.22.0 |
| QIIME2 | 2021.2 | 2021.8 |
| PICRUSt2 | 2.4.1 | 2.4.2 |
| MultiQC | 1.10.1 | 1.11 |
Removed
- #350 - Remove redundant derepFastq step (has no impact on results)
nf-core/ampliseq version 2.1.1 - 2021-10-28
Summary of changes
This release fixes bugs.
Detailed changes
Added
- #336 - Taxa agglomeration levels with
--dada_tax_agglom_min,--dada_tax_agglom_max,--qiime_tax_agglom_min,--qiime_tax_agglom_max, with defaults that go to genus level for abundance tables and ANCOM analysis
Changed
- 338 - Write empty space instead of
NAfor missing values in output files. - 342 - Added PICRUSt2 to summary figure.
Fixed
nf-core/ampliseq version 2.1.0 “Gray Steel Boa” - 2021-09-14
Summary of changes
- Updated documentation
- Added export to the Swedish Biodiversity Infrastructure (SBDI) on demand (
--sbdiexport) - Added newer versions of reference databases
- Predict community functions on demand (
--picrust) for 16S rRNA amplicons with PICRUSt2 - Fixed some bugs
Detailed changes
Added
- #322 - Export to the Swedish Biodiversity Infrastructure (SBDI)
- #294 - New version of the PR2 taxonomy database: 4.14.0, contains also some Bacteria, see PR2 release 4.14.0 notes
- #294 - New version of the Unite taxonomy databases: 8.3
- #302 - Pipeline workflow figure in README.md
- #307 - Functional predictions with PICRUSt2, on demand with
--picrust - #310 - Workflow figure in usage.md when using
--multiple_sequencing_runs - #312 - Added curated GTDB 16S taxonomy:
sbdi-gtdbas parameter to--dada_ref_taxonomy - #318 - Output information about the used reference taxonomy in a separate file in results folder
dada2/orqiime2/taxonomy
Changed
- #313 - Relative abundance tables in
qiime2/rel_abundance_tables/on ASV level were renamed and with either DADA2 (rel-table-ASV_with-DADA2-tax.tsv) or QIIME2 classifications (rel-table-ASV_with-QIIME2-tax.tsv), if available.
Fixed
- #306 - Sample names can now be identical to basenames of read files
- #299, #301- Template update for nf-core/tools version 2.1
- #303 - Reverse primer of PacBio and IonTorrent reads should be now given to
--RV_primerin usual direction (before: reverse complement) - #305 -
--max_lennow accepts integers as expected - #314, #315 - ASV fasta input via —input fixed (was broken in 2.0.0) and a test profile was added.
Dependencies
nf-core/ampliseq version 2.0.0 “Blue Copper Kangaroo” - 2021-06-29
Summary of changes
- Re-wrote whole pipeline in nextflow DSL2 instead of DSL1
- Three types of input is now supported: data folder, sample sheet, ASV fasta file
- ASV generation is now performed by DADA2 in R instead of DADA2 in QIIME2
- More reference taxonomy databases are now supported out of the box (SILVA, UNITE, PR2, GTDB, RDP)
- More data types can be analysed (Illumina SE & PE, PacBio, IonTorrent)
Detailed changes
Added
- #229 -
--single_endfor single-ended Illumina data - #229, #245, #267 - Taxonomic classification with DADA2
- #229 -
--dada_ref_taxonomyfor taxonomic classification with DADA2’s assignTaxonomy and addSpecies functions - #278 -
--qiime_ref_taxonomyfor taxonomic classification with QIIME2 - #239 - Support of RDP database for DADA2 classification
- #237 - Support of UNITE database for DADA2 classification
- #229 -
--inputmay point (1) at a fasta file ending with.fasta/.fna/.fathat will be taxonomically classified, (2) at a samples sheet ending with.tsvthat allows analysis of multiple sequencing runs by reading the optional columnrun, or (3) at a folder input - #229 -
--sample_inference,--concatenate_reads,--illumina_pe_its; please check the documentation for their function - #275 - Read count summary
- #274 -
--skip_qiimeto prevent any steps that are executed with QIIME2 - #272 -
--cut_itsto cut ASV sequence to ITS region before performing taxonomic classification with DADA2 - #280 - Added support for IonTorrent data
- #283 -
--cut_dada_ref_taxonomyallows extracting expected amplicons from DADA2 reference taxonomy database
Changed
- #254 - Updated CamelCase parameters to be lower_case_snake_case:
multipleSequencingRunstomultiple_sequencing_runsminLentomin_lenmaxLentomax_lenmaxEEtomax_ee
- #277 - Requires nextflow version
>= 21.04.0
Fixed
- #273 - Template update for nf-core/tools version 1.14
Dependencies
- #272 - New dependency ITSx v1.1.3
- #229 - Updated from cutadapt v2.8 to v3.2
- #229 - Updated DADA2 from v1.10 to v1.18.0, now not using QIIME2 for ASV generation any more
- #229 - Updated QIIME2 to v2021.2
Removed
- #229 -
--manifestis superseeded by--inputthat can now also handle a sample sheet file input (required extension:.tsv) - #229 -
--Q2importedanduntilQ2importare removed because pausing at that point is not neccessary - #229 -
--splitis no longer supported, therefore all sample IDs have to be unique - #229 -
--classifier_removeHashand--qiime_timezonebecame unnecessary - #229 -
--onlyDenoisingis deprecated in favour of--skip_taxonomy(which does the exact same thing) --taxon_referencebecame unnecessary- #229 -
--reference_databaseand--dereplicationare not supported any more.--qiime_ref_taxonomyallows now choosing a taxonomic reference
nf-core/ampliseq version 1.2.0 “Teal Bronze Lion” - 2021
Added
- #106 - Added support for PacBio data
- Added
--taxon_referenceto be able to support both ‘silva’ and ‘unite’ - #157 - Added possibility to run double cutadapt steps,
--double_primer - #211 - Added quality filter option
--maxEE
Fixed
- #182 - Fix input in case there are no underscores in sample IDs
- #186 - Update github actions
- #187 - Sample ids are in incorrect order in feature-table from PacBio data
- #201 - Template update for nf-core/tools version 1.12.1
- #147 - Split
make_classifierin two different processes that can be allocated different resources - #183 - Don’t fetch taxonomy/create classifier when run with
--skip_taxonomy - #180 - MultiQC, cutadapt and fastQC now work with
--multipleSequencingRuns
Dependencies
- Updated from cutadapt v2.6 to v2.8
nf-core/ampliseq version 1.1.3 - 2020
Added
- #170 - Cite paper for initial release
- #111 - Added parameter for user specified manifest file
- #118 - Added social preview images
- #135 - Added
--trunc_rminto make sure that auto trunc cutoff retaines a certain fraction of reads
Fixed
- #172 - Template update for nf-core/tools v1.11
- #163 - Template update for nf-core/tools v1.10.2
- #136 - Pipeline fails with remote working directory
- #152 - Don’t fetch taxonomy/create classifier when run with
--onlyDenoising
Dependencies
- Updated from MultiQC v1.6 to v1.9
Deprecated
--readsis replaced by--inputdue to nf-core/tools v1.10.2
nf-core/ampliseq version 1.1.2 - 2019
- No further changes, except a bugfix for the timezone issue found by @marchoeppner
- Specification of ‘—qiime_timezone’ might be required to run the analysis appropriately
nf-core/ampliseq version 1.1.1 - 2019
Pipeline Updates
- Update from QIIME2 v2018.6 to v2019.10, including DADA2 v1.6 to DADA2 v1.10
Bugfixes
- #78 - All sequenced classified to the same species
nf-core/ampliseq version 1.1.0 “Silver Lime Bee” - 2019
Pipeline updates
- #40 - Added support for data originating from multiple sequencing runs
- #53 - DADA2 report is always exported
- #49 - Allowed more filtering options
- #5 - Introduced check for existence of input files
- Extended parameter sanity check, including #15
- #61 - Improved documentation
- #62 - Utilize nf-core/configs centrally for this pipeline
- #63 - QIIME imports files by using a manifest, giving more freedom with input file names
- #84 - Add proper nf-core logo
Bug fixes
Initial release of nfcore/ampliseq!
The pipeline is a bioinformatics analysis pipeline used for 16S rRNA amplicon sequencing data.
The workflow processes raw data from FastQ inputs (FastQC), trims primer sequences from the reads (Cutadapt), imports data into QIIME2, generates amplicon sequencing variants (ASV, DADA2), classifies features against the SILVA v132 database, excludes unwanted taxa, produces absolute and relative feature/taxa count tables and plots, plots alpha rarefaction curves, computes alpha and beta diversity indices and plots thereof, and finally calls differentially abundant taxa (ANCOM). See the output documentation for more details of the results.