Back

Laemple: A Benchmarking Framework for Virus Lineage Deconvolution Tools for SARS-CoV-2 from Wastewater

Schedl, A.; Bergthaler, A.; Amman, F.

2026-09-07 public and global health
10.64898/2026.09.02.26362081 medRxiv
Show abstract

Background. Correct and accurate deconvolution of SARS-CoV-2 lineages from wastewater sequencing data is a challenging task, given the intricacy of wastewater amplicon-sequencing data and the ever-growing complexity of the lineage classification. Existing benchmarking studies made use of artificial spike-in compositions, thereby falling short of reflecting the prevailing complexity of wastewater samples. Results. We present a modular, expandable, simulation-based benchmarking framework as a reproducible Snakemake workflow, named Laemple, to evaluate the performance of virus lineage deconvolution tools. Using in silico simulated data sets of varying complexity and sequencing quality, we demonstrate its utility by evaluating seven publicly available tools, based on their precision, sensitivity, and reproducibility. Freyja showed robust sensitivity and consistent performance across diverse data set complexities, alongside user-friendly installation and documentation, while VaQuERo demonstrated the highest precision. Conclusions. Our results reveal substantial variation in tool performance across conditions, emphasizing the need to benchmark with diverse and complex scenarios. This framework enables informed tool selection for researchers and public health agencies and allows developers to stress-test their tools during development and maintenance, i.e., updating the lineage definition for newly emerging virus lineages. To this end, Laemple was designed in a modular fashion for customization and future expansion to support ongoing software development across all stages of the application life-cycle management.

Matching journals

The top 7 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.