Benchmark of thirteen bioinformatic pipelines for metagenomic virus diagnostics using datasets from clinical samples

Jutte J.C. de Vries*, Julianne R. Brown, Nicole Fischer, Igor A. Sidorov, Sofia Morfopoulou, Jiabin Huang, Bas B.Oude Munnink, Arzu Sayiner, Alihan Bulgurcu, Christophe Rodriguez, Guillaume Gricourt, Els Keyaerts, Leen Beller, Claudia Bachofen, Jakub Kubacki, Cordey Samuel, Laubscher Florian, Schmitz Dennis, Martin Beer, Dirk HoeperMichael Huber, Verena Kufner, Maryam Zaheri, Aitana Lebrand, Anna Papa, Sander van Boheemen, Aloys C.M. Kroes, Judith Breuer, F. Xavier Lopez-Labrador, Eric C.J. Claas

*Corresponding author for this work

Research output: Contribution to journalArticleAcademicpeer-review

12 Citations (Scopus)

Abstract

Introduction: Metagenomic sequencing is increasingly being used in clinical settings for difficult to diagnose cases. The performance of viral metagenomic protocols relies to a large extent on the bioinformatic analysis. In this study, the European Society for Clinical Virology (ESCV) Network on NGS (ENNGS) initiated a benchmark of metagenomic pipelines currently used in clinical virological laboratories. Methods: Metagenomic datasets from 13 clinical samples from patients with encephalitis or viral respiratory infections characterized by PCR were selected. The datasets were analyzed with 13 different pipelines currently used in virological diagnostic laboratories of participating ENNGS members. The pipelines and classification tools were: Centrifuge, DAMIAN, DIAMOND, DNASTAR, FEVIR, Genome Detective, Jovian, MetaMIC, MetaMix, One Codex, RIEMS, VirMet, and Taxonomer. Performance, characteristics, clinical use, and user-friendliness of these pipelines were analyzed. Results: Overall, viral pathogens with high loads were detected by all the evaluated metagenomic pipelines. In contrast, lower abundance pathogens and mixed infections were only detected by 3/13 pipelines, namely DNASTAR, FEVIR, and MetaMix. Overall sensitivity ranged from 80% (10/13) to 100% (13/13 datasets). Overall positive predictive value ranged from 71-100%. The majority of the pipelines classified sequences based on nucleotide similarity (8/13), only a minority used amino acid similarity, and 6 of the 13 pipelines assembled sequences de novo. No clear differences in performance were detected that correlated with these classification approaches. Read counts of target viruses varied between the pipelines over a range of 2-3 log, indicating differences in limit of detection. Conclusion: A wide variety of viral metagenomic pipelines is currently used in the participating clinical diagnostic laboratories. Detection of low abundant viral pathogens and mixed infections remains a challenge, implicating the need for standardization and validation of metagenomic analysis for clinical diagnostic use. Future studies should address the selective effects due to the choice of different reference viral databases.

Original languageEnglish
Article number104908
JournalJournal of Clinical Virology
Volume141
DOIs
Publication statusPublished - 1 Aug 2021

Bibliographical note

Funding Information:
MH was supported by the Clinical Research Priority Program 'Comprehensive Genomic Pathogen Detection' of the University of Zurich. FXL received funding from Instituto de Salud Carlos III, Spain (Grant numbers PI18/01824 and PI18/01759 and CIBEResp ).

Funding Information:
We thank the COMPARE study group (https://www.compare-europe.eu/) and the EMBL-EBI (https://www.ebi.ac.uk/) for the availability of the Data Hubs. MH was supported by the Clinical Research Priority Program 'Comprehensive Genomic Pathogen Detection' of the University of Zurich. FXL received funding from Instituto de Salud Carlos III, Spain (Grant numbers PI18/01824 and PI18/01759 and CIBEResp).

Publisher Copyright:
© 2021

Fingerprint

Dive into the research topics of 'Benchmark of thirteen bioinformatic pipelines for metagenomic virus diagnostics using datasets from clinical samples'. Together they form a unique fingerprint.

Cite this