Reproducibly sampling SARS-CoV-2 genomes across time, geography, and viral diversity
Open Access
- 28 October 2020
- journal article
- research article
- Published by F1000 Research Ltd in F1000Research
- Vol. 9, 657
Abstract
The COVID-19 pandemic has led to a rapid accumulation of SARS-CoV-2 genomes, enabling genomic epidemiology on local and global scales. Collections of genomes from resources such as GISAID must be subsampled to enable computationally feasible phylogenetic and other analyses. We present genome-sampler, a software package that supports sampling collections of viral genomes across multiple axes including time of genome isolation, location of genome isolation, and viral diversity. The software is modular in design so that these or future sampling approaches can be applied independently and combined (or replaced with a random sampling approach) to facilitate custom workflows and benchmarking. genome-sampler is written as a QIIME 2 plugin, ensuring that its application is fully reproducible through QIIME 2’s unique retrospective data provenance tracking system. genome-sampler can be installed in a conda environment on macOS or Linux systems. A complete default pipeline is available through a Snakemake workflow, so subsampling can be achieved using a single command. genome-sampler is open source, free for all to use, and available at https://caporasolab.us/genome-sampler. We hope that this will facilitate SARS-CoV-2 research and support evaluation of viral genome sampling approaches for genomic epidemiology.Keywords
Funding Information
- Arizona Board of Regents
- Australian National University
- National Cancer Institute (U54CA143925)
- Alfred P Sloan Foundation (G-2019-12432)
- National Cancer Institute (U54CA143924)
- National Science Foundation (1565100)
- Chan-Zuckerberg Initiative
This publication has 11 references indexed in Scilit:
- Defining the Pandemic at the State Level: Sequence-Based Epidemiology of the SARS-CoV-2 virus by the Arizona COVID-19 Genomics Union (ACGU)medRxiv, 2020
- Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2Nature Biotechnology, 2019
- Nextstrain: real-time tracking of pathogen evolutionBioinformatics, 2018
- The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 updateNucleic Acids Research, 2018
- Data, disease and diplomacy: GISAID's innovative contribution to global healthGlobal Challenges, 2017
- VSEARCH: a versatile open source tool for metagenomicsPeerJ, 2016
- MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and UsabilityMolecular Biology and Evolution, 2013
- GenBankNucleic Acids Research, 2012
- Snakemake—a scalable bioinformatics workflow engineBioinformatics, 2012
- Nomenclature for incompletely specified bases in nucleic acid sequences: rcommendations 1984Nucleic Acids Research, 1985