Hecatomb: an integrated software platform for viral metagenomics

Michael J Roach,Sarah J Beecroft,Kathie A Mihindukulasuriya,Leran Wang,Anne Paredes,Luis Alberto Chica Cárdenas,Kara Henry-Cocks,Lais Farias Oliveira Lima,Elizabeth A Dinsdale,Robert A Edwards,Scott A Handley
DOI: https://doi.org/10.1093/gigascience/giae020
IF: 7.658
2024-06-06
GigaScience
Abstract:Modern sequencing technologies offer extraordinary opportunities for virus discovery and virome analysis. Annotation of viral sequences from metagenomic data requires a complex series of steps to ensure accurate annotation of individual reads and assembled contigs. In addition, varying study designs will require project-specific statistical analyses. Here we introduce Hecatomb, a bioinformatic platform coordinating commonly used tasks required for virome analysis. Hecatomb means "a great sacrifice." In this setting, Hecatomb is "sacrificing" false-positive viral annotations using extensive quality control and tiered-database searches. Hecatomb processes metagenomic data obtained from both short- and long-read sequencing technologies, providing annotations to individual sequences and assembled contigs. Results are provided in commonly used data formats useful for downstream analysis. Here we demonstrate the functionality of Hecatomb through the reanalysis of a primate enteric and a novel coral reef virome. Hecatomb provides an integrated platform to manage many commonly used steps for virome characterization, including rigorous quality control, host removal, and both read- and contig-based analysis. Each step is managed using the Snakemake workflow manager with dependency management using Conda. Hecatomb outputs several tables properly formatted for immediate use within popular data analysis and visualization tools, enabling effective data interpretation for a variety of study designs. Hecatomb is hosted on GitHub (github.com/shandley/hecatomb) and is available for installation from Bioconda and PyPI.
multidisciplinary sciences
What problem does this paper attempt to address?