Abstract:Abstract Background In the last few years, 16S rRNA gene sequencing (16S rDNA-seq) has seen a surprisingly rapid increase in election rate as a methodology to perform microbial community studies. Despite the considerable popularity of this technique, an exiguous number of specific tools are currently available for proper 16S rDNA-seq count data preprocessing and simulation. Indeed, the great majority of tools have been developed adapting methodologies previously used for bulk RNA-seq data, with poor assessment of their applicability in the metagenomics field. For such tools and the few ones specifically developed for 16S rDNA-seq data, performance assessment is challenging, mainly due to the complex nature of the data and the lack of realistic simulation models. In fact, to the best of our knowledge, no software thought for data simulation are available to directly obtain synthetic 16S rDNA-seq count tables that properly model heavy sparsity and compositionality typical of these data. Results In this paper we present metaSPARSim, a sparse count matrix simulator intended for usage in development of 16S rDNA-seq metagenomic data processing pipelines. metaSPARSim implements a new generative process that models the sequencing process with a Multivariate Hypergeometric distribution in order to realistically simulate 16S rDNA-seq count table, resembling real experimental data compositionality and sparsity. It provides ready-to-use count matrices and comes with the possibility to reproduce different pre-coded scenarios and to estimate simulation parameters from real experimental data. The tool is made available at http://sysbiobig.dei.unipd.it/?q=Software#metaSPARSim and https://gitlab.com/sysbiobig/metasparsim . Conclusion metaSPARSim is able to generate count matrices resembling real 16S rDNA-seq data. The availability of count data simulators is extremely valuable both for methods developers, for which a ground truth for tools validation is needed, and for users who want to assess state of the art analysis tools for choosing the most accurate one. Thus, we believe that metaSPARSim is a valuable tool for researchers involved in developing, testing and using robust and reliable data analysis methods in the context of 16S rRNA gene sequencing.

Readsynth: short-read simulation for consideration of composition-biases in reduced metagenome sequencing approaches

Reduced metagenome sequencing for strain-resolution taxonomic profiles

A Comprehensive Evaluation of Alignment Software for Reduced Representation Bisulfite Sequencing Data

SimReadUntil for Benchmarking Selective Sequencing Algorithms on ONT Devices

scReadSim: a single-cell RNA-seq and ATAC-seq read simulator

metaSPARSim: a 16S rRNA gene sequencing count data simulator

Semisynthetic Simulation for Microbiome Data Analysis

Characterization and simulation of metagenomic nanopore sequencing data with Meta-NanoSim

SimpleMetaPipeline: Breaking the bioinformatics bottleneck in metabarcoding

BEERS2: RNA-Seq simulation through high fidelity in silico modeling

MetaSMC: a Coalescent-Based Shotgun Sequence Simulator for Evolving Microbial Populations

AsaruSim: a single-cell and spatial RNA-Seq Nanopore long-reads simulation workflow

Systematic Review of Next-Generation Sequencing Simulators: Computational Tools, Features and Perspectives

Assessment of Metagenomic Assemblers Based on Hybrid Reads of Real and Simulated Metagenomic Sequences

Comprehensive Simulation of Metagenomic Sequencing Data with Non-Uniform Sampling Distribution

TKSM: highly modular, user-customizable, and scalable transcriptomic sequencing long-read simulator

Simple, reference-independent analyses help optimize hybrid assembly of microbial community metagenomes

MCSS: microbial community simulator based on structure

SQANTI-reads: a tool for the quality assessment of long read data in multi-sample lrRNA-seq experiments

METASEED: a novel approach to full-length 16S rRNA gene reconstruction from short read data

SEQ2MGS: an effective tool for generating realistic artificial metagenomes from the existing sequencing data