Abstract:Abstract Background Environmental DNA and metabarcoding allow the identification of a mixture of species and launch a new era in bio- and eco-assessment. Many steps are required to obtain taxonomically assigned matrices from raw data. For most of these, a plethora of tools are available; each tool's execution parameters need to be tailored to reflect each experiment's idiosyncrasy. Adding to this complexity, the computation capacity of high-performance computing systems is frequently required for such analyses. To address the difficulties, bioinformatic pipelines need to combine state-of-the art technologies and algorithms with an easy to get-set-use framework, allowing researchers to tune each study. Software containerization technologies ease the sharing and running of software packages across operating systems; thus, they strongly facilitate pipeline development and usage. Likewise programming languages specialized for big data pipelines incorporate features like roll-back checkpoints and on-demand partial pipeline execution. Findings PEMA is a containerized assembly of key metabarcoding analysis tools that requires low effort in setting up, running, and customizing to researchers’ needs. Based on third-party tools, PEMA performs read pre-processing, (molecular) operational taxonomic unit clustering, amplicon sequence variant inference, and taxonomy assignment for 16S and 18S ribosomal RNA, as well as ITS and COI marker gene data. Owing to its simplified parameterization and checkpoint support, PEMA allows users to explore alternative algorithms for specific steps of the pipeline without the need of a complete re-execution. PEMA was evaluated against both mock communities and previously published datasets and achieved results of comparable quality. Conclusions A high-performance computing–based approach was used to develop PEMA; however, it can be used in personal computers as well. PEMA's time-efficient performance and good results will allow it to be used for accurate environmental DNA metabarcoding analysis, thus enhancing the applicability of next-generation biodiversity assessment studies.

pipemake: A pipeline creation tool using Snakemake for reproducible analysis of biological datasets

Panpipes: a Pipeline for Multiomic Single-Cell and Spatial Transcriptomic Data Analysis

Sustainable data analysis with Snakemake

Biopipe: A Lightweight System Enabling Comparison of Bioinformatics Tools and Workflows

Pipeasm: a tool for automated large chromosome-scale genome assembly and evaluation

Dadasnake, a Snakemake implementation of DADA2 to process amplicon sequencing data for microbial ecology

VIPER: Visualization Pipeline for RNA-seq, a Snakemake workflow for efficient and complete RNA-seq analysis

hppRNA-a Snakemake-based handy parameter-free pipeline for RNA-Seq analysis of numerous samples

SimpleMetaPipeline: Breaking the bioinformatics bottleneck in metabarcoding

Snaq: A Dynamic Snakemake Pipeline for Microbiome Data Analysis With QIIME2

Spacemake: processing and analysis of large-scale spatial transcriptomics data

A snakemake toolkit for the batch assembly, annotation, and phylogenetic analysis of mitochondrial genomes and ribosomal genes from genome skims of museum collections

Scalepopgen: bioinformatics workflow resources implemented in Nextflow for comprehensive population genomic analyses

ExplorePipolin: reconstruction and annotation of piPolB-encoding bacterial mobile elements from draft genomes

PEMA: a flexible Pipeline for Environmental DNA Metabarcoding Analysis of the 16S/18S ribosomal RNA, ITS, and COI marker genes

How do users design scientific workflows? The Case of Snakemake

Biopipe: A Flexible Framework for Protocol-Based Bioinformatics Analysis

MetaWRAP—a flexible pipeline for genome-resolved metagenomic data analysis

PAPipe: a pipeline for comprehensive population genetic analysis

Accelerating Single-Cell Sequencing Data Analysis with SciDAP: A User-Friendly Approach

A Snakemake Toolkit for the Batch Assembly, Annotation and Phylogenetic Analysis of Mitochondrial Genomes and Ribosomal Genes From Genome Skims of Museum Collections