The ViReflow Pipeline Enables User Friendly Large Scale Viral Consensus Genome Reconstruction

Overview

Journal Sci Rep

Specialty Science

Date 2022 Mar 25

PMID 35332213

Authors

Niema Moshiri

Kathleen M Fisch

Amanda Birmingham

Peter DeHoff

Gene W Yeo

Kristen Jepsen

Louise C Laurent

Rob Knight

Affiliations

Soon will be listed here.

Abstract

Throughout the COVID-19 pandemic, massive sequencing and data sharing efforts enabled the real-time surveillance of novel SARS-CoV-2 strains throughout the world, the results of which provided public health officials with actionable information to prevent the spread of the virus. However, with great sequencing comes great computation, and while cloud computing platforms bring high-performance computing directly into the hands of all who seek it, optimal design and configuration of a cloud compute cluster requires significant system administration expertise. We developed ViReflow, a user-friendly viral consensus sequence reconstruction pipeline enabling rapid analysis of viral sequence datasets leveraging Amazon Web Services (AWS) cloud compute resources and the Reflow system. ViReflow was developed specifically in response to the COVID-19 pandemic, but it is general to any viral pathogen. Importantly, when utilized with sufficient compute resources, ViReflow can trim, map, call variants, and call consensus sequences from amplicon sequence data from 1000 SARS-CoV-2 samples at 1000X depth in < 10 min, with no user intervention. ViReflow's simplicity, flexibility, and scalability make it an ideal tool for viral molecular epidemiological efforts.

Citing Articles

ViralFlow v1.0-a computational workflow for streamlining viral genomic surveillance.

da Silva A, da Silva Neto A, Aksenen C, Jeronimo P, Dezordi F, Almeida S NAR Genom Bioinform. 2024; 6(2):lqae056.

PMID: 38800829 PMC: 11127631. DOI: 10.1093/nargab/lqae056.

A hepatitis B virus (HBV) sequence variation graph improves alignment and sample-specific consensus sequence construction.

Duchen D, Clipman S, Vergara C, Thio C, Thomas D, Duggal P PLoS One. 2024; 19(4):e0301069.

PMID: 38669259 PMC: 11051683. DOI: 10.1371/journal.pone.0301069.

ViralWasm: a client-side user-friendly web application suite for viral genomics.

Ji D, Aboukhalil R, Moshiri N Bioinformatics. 2024; 40(1).

PMID: 38200583 PMC: 10809900. DOI: 10.1093/bioinformatics/btae018.

COWID: an efficient cloud-based genomics workflow for scalable identification of SARS-COV-2.

Lim H, Fann Y, Lee Y Brief Bioinform. 2023; 24(5).

PMID: 37738400 PMC: 10516370. DOI: 10.1093/bib/bbad280.

A hepatitis B virus (HBV) sequence variation graph improves sequence alignment and sample-specific consensus sequence construction for genetic analysis of HBV.

Duchen D, Clipman S, Vergara C, Thio C, Thomas D, Duggal P bioRxiv. 2023; .

PMID: 36711598 PMC: 9882026. DOI: 10.1101/2023.01.11.523611.

References

Li D, Liu C, Luo R, Sadakane K, Lam T . MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015; 31(10):1674-6. DOI: 10.1093/bioinformatics/btv033. View

Topfer A, Marschall T, Bull R, Luciani F, Schonhuth A, Beerenwinkel N . Viral quasispecies assembly via maximal clique enumeration. PLoS Comput Biol. 2014; 10(3):e1003515. PMC: 3967922. DOI: 10.1371/journal.pcbi.1003515. View

Kim D, Paggi J, Park C, Bennett C, Salzberg S . Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019; 37(8):907-915. PMC: 7605509. DOI: 10.1038/s41587-019-0201-4. View

Chen S, Zhou Y, Chen Y, Gu J . fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. 2018; 34(17):i884-i890. PMC: 6129281. DOI: 10.1093/bioinformatics/bty560. View

Chikhi R, Rizk G . Space-efficient and exact de Bruijn graph representation based on a Bloom filter. Algorithms Mol Biol. 2013; 8(1):22. PMC: 3848682. DOI: 10.1186/1748-7188-8-22. View

Wick R, Judd L, Gorrie C, Holt K . Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol. 2017; 13(6):e1005595. PMC: 5481147. DOI: 10.1371/journal.pcbi.1005595. View

Shu Y, McCauley J . GISAID: Global initiative on sharing all influenza data - from vision to reality. Euro Surveill. 2017; 22(13). PMC: 5388101. DOI: 10.2807/1560-7917.ES.2017.22.13.30494. View

Li H . A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011; 27(21):2987-93. PMC: 3198575. DOI: 10.1093/bioinformatics/btr509. View

Langmead B, Salzberg S . Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012; 9(4):357-9. PMC: 3322381. DOI: 10.1038/nmeth.1923. View

10.

Li H . Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018; 34(18):3094-3100. PMC: 6137996. DOI: 10.1093/bioinformatics/bty191. View

11.

Rambaut A, Holmes E, OToole A, Hill V, McCrone J, Ruis C . A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. Nat Microbiol. 2020; 5(11):1403-1407. PMC: 7610519. DOI: 10.1038/s41564-020-0770-5. View

12.

Dezordi F, da Silva Neto A, Campos T, Jeronimo P, Aksenen C, Almeida S . ViralFlow: A Versatile Automated Workflow for SARS-CoV-2 Genome Assembly, Lineage Assignment, Mutations and Intrahost Variant Detection. Viruses. 2022; 14(2). PMC: 8877152. DOI: 10.3390/v14020217. View

13.

Moshiri N, Smith D, Mirarab S . HIV Care Prioritization Using Phylogenetic Branch Length. J Acquir Immune Defic Syndr. 2021; 86(5):626-637. PMC: 7933099. DOI: 10.1097/QAI.0000000000002612. View

14.

Bhoyar R, Jain A, Sehgal P, Divakar M, Sharma D, Imran M . High throughput detection and genetic epidemiology of SARS-CoV-2 using COVIDSeq next-generation sequencing. PLoS One. 2021; 16(2):e0247115. PMC: 7888613. DOI: 10.1371/journal.pone.0247115. View

15.

Li H, Durbin R . Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009; 25(14):1754-60. PMC: 2705234. DOI: 10.1093/bioinformatics/btp324. View

16.

Bolger A, Lohse M, Usadel B . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30(15):2114-20. PMC: 4103590. DOI: 10.1093/bioinformatics/btu170. View

17.

Zagordi O, Bhattacharya A, Eriksson N, Beerenwinkel N . ShoRAH: estimating the genetic diversity of a mixed sample from next-generation sequencing data. BMC Bioinformatics. 2011; 12:119. PMC: 3113935. DOI: 10.1186/1471-2105-12-119. View

18.

OToole A, Scher E, Underwood A, Jackson B, Hill V, McCrone J . Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool. Virus Evol. 2021; 7(2):veab064. PMC: 8344591. DOI: 10.1093/ve/veab064. View

19.

Karthikeyan S, Nguyen A, McDonald D, Zong Y, Ronquillo N, Ren J . Rapid, Large-Scale Wastewater Surveillance and Automated Reporting System Enable Early Detection of Nearly 85% of COVID-19 Cases on a University Campus. mSystems. 2021; 6(4):e0079321. PMC: 8409724. DOI: 10.1128/mSystems.00793-21. View

20.

Baaijens J, El Aabidine A, Rivals E, Schonhuth A . De novo assembly of viral quasispecies using overlap graphs. Genome Res. 2017; 27(5):835-848. PMC: 5411778. DOI: 10.1101/gr.215038.116. View