# Metavisitor, a Suite of Galaxy Tools for Simple and Rapid Detection and Discovery of Viruses in Deep Sequence Data

Guillaume D Carissimo, Marius D van den Beek, Kenneth D Vernick, Christophe D Antoniewski

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Guillaume D Carissimo, Marius D van den Beek, Kenneth D Vernick, Christophe D Antoniewski. Metavisitor, a Suite of Galaxy Tools for Simple and Rapid Detection and Discovery of Viruses in Deep Sequence Data. PLoS ONE, 2017, 12 (1), pp.e0168397. 10.1371/journal.pone.0168397. hal-01465804

HAL Id: hal-01465804

https://hal.sorbonne-universite.fr/hal-01465804v1

Submitted on 13 Feb 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
RESEARCH ARTICLE

# Metavisitor, a Suite of Galaxy Tools for Simple and Rapid Detection and Discovery of Viruses in Deep Sequence Data

Guillaume Carissimo<sup>1,2,3*,</sup> Marius van den Beek<sup>4,5*</sup>, Kenneth D. Vernick<sup>1,2,6</sup>, Christophe Antoniewski<sup>4,5*</sup>

1 Institut Pasteur, Unit of Insect Vector Genetics and Genomics, Department of Parasites and Insect Vectors, Paris, FRANCE, 2 CNRS, Unit of Hosts, Vectors and Pathogens (URA3012), Paris, FRANCE, 3 Laboratory of Microbial Immunity, Singapore Immunology Network, A\*STAR, 8A Biomedical Grove, Biopolis, Singapore, Singapore, 4 Sorbonne Universités, Université Pierre et Marie Curie (UPMC), CNRS, Institut de Biologie Paris Seine (IBPS), Developmental Biology Department, Paris, France, 5 Sorbonne Universités, Université Pierre et Marie Curie (UPMC), CNRS, Institut de Biologie Paris Seine (IBPS), ARTbio Bioinformatics Analysis Facility, Paris, France, 6 Department of Microbiology, University of Minnesota, Minneapolis, MN, United States of America

© These authors contributed equally to this work.

* christophe.antoniewski@upmc.fr

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

## OPEN ACCESS

**Citation:** Carissimo G, van den Beek M, Vernick KD, Antoniewski C (2017) Metavisitor, a Suite of Galaxy Tools for Simple and Rapid Detection and Discovery of Viruses in Deep Sequence Data. PLoS ONE 12(1): e0168397. doi:10.1371/journal. pone.0168397

**Editor:** Sebastien Pfeffer, Institut de Biologie Moleculaire et Cellulaire, FRANCE

**Received:** May 20, 2016

**Accepted:** November 29, 2016

**Published:** January 3, 2017

**Copyright:** © 2017 Carissimo et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All data, codes and workflows are currently fully accessible in public repositories: https://toolshed.g2.bx.psu.edu/, https://github.com/ARTbio/tools-artbio, https://artbio.github.io/Metavisitor-manual/, https://mississippi.snv.jussieu.fr/, https://dx.doi.org/10.6084/m9 figshare.3179026.v1, http://www.myexperiment.org/users/65016/workflows.

**Funding:** This work received financial support to CA from the Agence Nationale de la Recherche grant ANR-13-BSV2-0007 « Plastisipi » and from

## Abstract

Metavisitor is a software package that allows biologists and clinicians without specialized bioinformatics expertise to detect and assemble viral genomes from deep sequence datasets. The package is composed of a set of modular bioinformatic tools and workflows that are implemented in the Galaxy framework. Using the graphical Galaxy workflow editor, users with minimal computational skills can use existing Metavisitor workflows or adapt them to suit specific needs by adding or modifying analysis modules. Metavisitor works with DNA, RNA or small RNA sequencing data over a range of read lengths and can use a combination of *de novo* and guided approaches to assemble genomes from sequencing reads. We show that the software has the potential for quick diagnosis as well as discovery of viruses from a vast array of organisms. Importantly, we provide here executable Metavisitor use cases, which increase the accessibility and transparency of the software, ultimately enabling biologists or clinicians to focus on biological or medical questions.

## Introduction

Viruses infect cells and manipulate the host machinery for their replication and transmission. Genomes of viruses show high diversity and can consist of single- or double-stranded RNA or DNA. Many types of viral replication cycles exist which may involve various cellular components, various DNA or RNA replication intermediates, and diverse strategies for viral RNA transcription and viral protein translation. Thus, deep-sequencing has become a powerful approach for virologists in their quest to detect and identify viruses in biological samples, even when they are present at low levels. Plants and invertebrates use RNA interference as an antiviral mechanism [1,2]. Antiviral RNAi activity results in accumulation of viral interfering small
the Institut Universitaire d'Ingénierie en Santé IUIS grant "ARTIMED" and to KDV from the European Research Council (http://erc.europa.eu/funding-and-grants). Support for frontier research, Advanced Grant #323173 AnoPath; European Commission (http://ec.europa.eu/research/fp7), FP7 Infrastructures #228421 Infravec; and French Laboratoire d'Excellence (http://www.enseignementsup-recherche.gouv.fr/cid55551/investissements-d-avenir-projets-laboratoires-d-excellence-par-region-et-domaine.html), "Integrative Biology of Emerging Infectious Diseases" #ANR-10-LABX-62-IBEID. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** CA is an academic editor at PLoS ONE.

RNAs (viRNAs), whose extent depends on several factors such as the ability of a virus to replicate in the host and to evade the host RNAi machinery. Moreover, viRNAs derived from a variety of viruses can be detected in host organisms, regardless if these viruses have positive single strand, negative single strand or double-stranded RNA genomes, or DNA genomes [2]. Together, these features make small RNA deep sequencing a potent approach to detect viruses regardless of their genomic specificities, and different bioinformatic tools have been developed for detection or *de novo* assembly of viral genomes.

Accordingly, viRNAs produced by the insect model *Drosophila melanogaster* in response to viral infections were sufficient to reconstruct and improve the genomic consensus sequence of the Nora virus [3] using the Paparazzi software [4] which is based on the SSAKE assembler [5]. In that study, Paparazzi improved the consensus sequence and the coverage of the Nora virus genome by ~20%, as compared to the previous Nora virus reference genome. SearchS-mallRNA, a standalone tool with a graphical interface, used a similar approach to reconstruct viral genomes [6]. Importantly, both programs require known, closely related viral references for proper guidance of genome reconstructions from viRNAs, precluding the identification of more distant viral species or discovery of novel or unexpected viruses.

To circumvent the need for viral reference sequences, Velvet [7] *de novo* assembled contigs from plant [8], fruit fly and mosquito [9] have been aligned to NCBI sequence databases, allowing the identification of partial or complete viral genomes. Several studies improved this strategy by combining two *de novo* assemblers [10–13], or scaffolding virus-aligned contigs using an additional translation-guided assembly step [14].

Collectively, these studies allowed important progress in virus assembly and identification from deep sequencing data. However, the existing computational workflows require specialist skills for installation, execution and adaptation to specific research, making them poorly accessible to a broad user base of biologists. In some cases, the tools lack documentation or their source codes are not available.

In this context, we developed Metavisor as an open source set of tools and preset workflows [15,16] which allow effective implementation of the computational strategies in the Galaxy framework, with short read as well as long read sequence datasets. In addition, Metavisor workflows can be easily adapted to suit specific needs, by adding analysis steps or replacing/ modifying existing ones with the numerous tools available in the Galaxy tool sheds. Here, we report a series of use cases of Metavisor and we show that it provides biologists and medical practitioners with an easy-to-use and adaptable software for the detection or identification of viruses from high-throughput sequence datasets.

## Experimental Procedures

Metavisor consists of a set of Galaxy tools (Fig 1) that can be combined to (i) retrieve up-to-date nucleotide as well as protein sequences of viral genomes deposited in Genbank [17] and index these sequences for subsequent alignments; (ii) extract sequencing reads that do not align to the host genomes, known symbionts or parasites; (iii) perform *de novo* assembly of these reads using assembly tools available in Galaxy, align the *de novo* contigs against the viral nucleotide or protein blast databases using blastn or blastx, respectively, and generate reports from blast outputs to help in known viruses diagnosis or in candidate virus discovery; (iv) use CAP3 (optional, see Use Case 3–3), blast and viral scaffolds for selected viruses to generate guided final viral sequence assemblies of blast sequence hits. Below, we group analysis steps in functional tasks i to iv and provide details on the Metavisor tools. These tasks are linked together to build full workflows adapted to the analysis of the use cases described in the result section.
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Global view of the Metavisitor workflow. The workflow is organised in sub workflows (dashed line) corresponding to functional tasks as described in the manuscript. All Galaxy Tools (square boxes) are available in the main Galaxy tool shed (https://toolshed.g2.bx.psu.edu/).</figcaption></figure>

doi:10.1371/journal.pone.0168397.g001
### (i) Get reference viral sequences

The “Get reference viral sequences” task is performed using the “Retrieve FASTA from NCBI” tool that sends a query string to the Genbank database [17] and retrieves the corresponding nucleotide or protein sequences. For the viral nucleotide and protein sequences referred to as “vir1”, we used the tool to query Genbank (oct 2015) and retrieve viruses sequences filtered out from cellular organisms and bacteriophage sequences (see S1 Fig). However, users can change the tool settings by entering query strings that fit their specific needs. As retrieving vir1 from NCBI takes several hours, we allow users to skip the step by directly accessing the nucleotides or protein vir1 datasets on the Mississippi server (http://mississippi.fr) or to download them from figshare (https://dx.doi.org/10.6084/m9 figshare.3179026). For convenience, nucleotide and protein blast indexes of vir1 are also available in the public library of the Mississippi server, but can also be generated using the “NCBI BLAST+ makeblastdb” Galaxy tool [18]. Bowtie [19] as well as bowtie2 [20] indexes of the vir1 nucleotide sequences have been generated in the Mississippi Galaxy instance using the corresponding “data manager” Galaxy tools.

Finally, users can upload their own viral nucleotide and protein sequences using ftp and transfer them to a Galaxy history (Fig 1), where they can use the Galaxy data manager tools to produce the blast and bowtie indexes necessary for Metavisor.

### (ii) Prepare data

The “Prepare data” task (Fig 1) processes Illumina sequencing datasets in order to optimize the subsequent *de novo* assembly of viral sequencing reads. Fastq files of sequence reads are first clipped from library adapters and converted to fasta format using our tool “Clip adapter” tool (S1 Table). The clipped reads may be further converted to a fasta file of unique sequences headed by a character string that contains a unique identifier and the number of times that the sequences were found in the dataset, thus reducing the size of the dataset without loss of information. This optional treatment removes sequence duplicates and drastically reduces the workload of the next steps as well as the coverage variations after *de novo* assembly (see Use Cases 1–1 to 1–3). Clipped reads are then depleted from non-viral sequences by sequential alignments to the host genome, to other genomes from known or potential symbionts and parasites, as well as to PhiX174 genome sequences which are commonly used as internal controls in Illumina sequencing and may contaminate the datasets (Fig 1). The sequence reads that did not match the reference genomes are retained and returned as a fasta file that can be used subsequently by a *de novo* assembly tool. Note that these subtraction steps can be skipped when the host genome is not known or if the aim of the user is to discover endogenous viral elements [21].

### (iii) Assembly, Blast and Parsing

**De novo assembly.** In the task “Assemble, Blast and Parse” (iii), retained RNA sequences are subjected to *de novo* assembly. For short reads (<50 nt), we tested several rounds of *de novo* assembly by Velvet [7] using the Oases software package [22] (Fig 1) and k-mer lengths ranging from 15 to 35 (S1 Table). For reads between 50 nt and 100 nt, we also used the Oases with k-mer lengths ranging from 13 to 69. Finally in Use Case 3–3, we used the Trinity assembly software which is available as a Galaxy tool and was reported to perform well with long reads [23]. Trinity as well as SPAdes [24] assembly softwares were also tested as alternate option to Oases in the Use Case 2–2 (S1 Table), giving similar outputs. It is noteworthy that users can adapt a Metavisor workflow using any assembly software available in the Galaxy tool shed.
**Blast.** Next, *de novo* assembled contigs are aligned to both nucleotide and protein vir1 BLAST databases built from the viral reference sequences (Fig 1) using the blastn or blastx Galaxy tools [18]. These tools search nucleotide or protein databases using nucleotide or translated nucleotide queries, respectively [25]. The default parameters are adjusted in order to report only the 5 best alignments per contig (Maximum hits option is set to 5) and to generate a tabular blast output that includes the 12 standard columns plus a column containing the length of the aligned subject sequences (extended columns option, “slen” checked). Note that this additional column in the blast output is required for subsequent parsing of the blast output by the “Parse blast output and compile hits” tool.

**Parsing.** Tabular outputs generated by blastn and blastx alignments are processed by the “Parse blast output and compile hits” tool (S1 Table), which returns 4 files, namely “blast analysis, by subjects”, “hits”, “Blast aligned sequences” and “Blast unaligned sequences”.

In the “blast analysis, by subjects” file (S2 Fig), the subject sequences in the viral nucleotide or protein blast databases that produced significant blast alignments (hits) with *de novo* assembled contigs are listed, together with those contigs and hit information (% Identity, Alignment Length, start and end coordinates of hits relatively to the subject sequence, percentage of the contig length covered by the hit, E-value and Bit Score of the hit). In addition, for each subject sequence in the list, the length in nucleotide or amino-acid of the subject sequence (Subject Length), the summed coverage of the subject by all contig hits (Total Subject Coverage) as well as the fraction of the subject length that this coverage represents (Relative Subject Coverage), and the best (Best Bit Score) and mean (Mean Bit Score) bit scores produced by contig hits are computed and indicated. A simplified output can be generated without contigs and blast information by using the “compact” option for the reporting mode of the “Parse blast output and compile hits” tool. Note that the total and relative subject coverages indicate how much of the virus sequence is covered by the reconstructed contigs, whereas the Bit scores allow to estimate the distances between the reconstructed contigs and the subject sequence.

The “hits” file contains the sequences of contig portions that produced significant alignment in the BLAST step (i.e. query hit sequences), flanked by additional contig nucleotides 5’ and 3’ to the hit (the size of these margins is set to 5 by default and can be modified by the user). These margins allow to include sequences that might not have significant homology but could still be of viral origin.

Finally, the “Blast aligned sequences” file contains contigs that produced significant blast hits, whereas the “Blast unaligned sequences” file contains those that did not.

#### (iv) Blast-Guided Scaffolding

This last task allows to integrate hit sequences matching a candidate virus into a virus scaffold (Fig 1). First, blastn or blastx hits are retrieved from the “hits” file using the tool “Pick Fasta sequences” (S1 Table) and the appropriate query string (for instance, “Dengue” will retrieve hit sequences that significantly blast aligned with Dengue virus sequences). Next, these hit sequences can be further clustered in longer contigs using the “cap3 Sequence Assembly” Galaxy tool (S1 Table) adapted from CAP3 [26]. Finally, if there are still multiple unlinked contigs at this stage, they can be integrated (uppercase characters) in the matched viral sequence taken as a scaffold (lowercase characters). This scaffolding is achieved by (a) retrieving the viral sequence from the NCBI nucleotide database to be used as the backbone of the scaffold, generating a blast index from this sequence and aligning the contigs to this index with blastn or tblastx tools (b) running the “blast_to_scaffold” tools (S1 Table), taking as inputs the contigs, the viral guide sequence and the blastn or blastx output (Fig 1, bottom).
Availability of Metavisor

All Metavisor tools, workflows and use cases are available on the Galaxy server http://
mississippi.snv.jussieu.fr. Readers can import in their personal account the published Metavi-
sitor use case histories and their corresponding workflows to re-run the described analyses or
adapt them to their studies.

We made all tools and workflows that compose Metavisor available from the main Galaxy
tool shed (https://toolshed.g2.bx.psu.edu/), in the form of a tool suite (suite_metavisor_1_2)
which thus can be installed and used on any Galaxy server instance. The Metavisor workflows
are also available from the *myexperiment* repository (http://www.myexperiment.org/) They
can be freely modified or complemented with additional analysis steps within the Galaxy
environment.

The Metavisor tool codes are accessible in our public GitHub repository (https://github.com/ARTbio/tools-artbio/). We also provide a Docker image artbio/metavisor:1.2 as well as an ansible playbook that both allow to deploy a Galaxy server instance with preinstalled Metavisor tools and workflows in local infrastructures. Extensive documentation on how to install and use Metavisor is available at https://artbio.github.io/Metavisor-manual/.

Results

The strategy implemented by Metavisor (Fig 1) is to perform *de novo* assembly of sequencing
reads and to detect contigs of viral origin through blast alignments to a nucleotide or protein
sequence database of known viruses (vir1). These contig alignments can be further clustered to
reconstruct a viral genome.

Below, we report use cases to demonstrate the use of Metavisor in specific situations. For
each use case, we briefly present the purpose of the original study from which the datasets orig-
inate and we describe an adapted Metavisor workflow as well as its main outputs. Readers
can further examine the workflows (https://mississippi.snv.jussieu.fr/workflow/listpublished)
and use case analyses (https://mississippi.snv.jussieu.fr/history/listpublished) in every detail
at http://mississippi.fr. Indicative execution times of the workflows are given in S2 Table.

1. Detection of known viruses

**Use Cases 1–1, 1–2 and 1–3: detection and reconstruction of the Nora virus genome in small RNA sequencing datasets.** Using small RNA sequencing libraries SRP013822 (EBI ENA) and the Paparazzi software [4] we were previously able to propose a novel reference genome (NCBI JX220408) for the Nora virus strain infecting *Drosophila melanogaster* stocks in laboratories [3]. This so-called rNora genome differed by 3.2% nucleotides from the Nora virus reference NC_007919.3 and improved the alignment rate of viral siRNAs by ~121%. Thus, we first tested Metavisor on the small RNA sequencing datasets SRP013822 using the Oases *de novo* assembly tool which is well suited to assembly of short read [9].

Three Metavisor workflows were run on the merged SRP013822 small RNA sequence
reads and the NC_007919.3 genome as a guide for final scaffolding. The workflow for Use
Case 1–1 (S3 Fig) used raw reads collapsed to unique sequences (experimental procedures) to
reconstruct a Nora virus genome referred to as Nora_MV (S1 File). In a second workflow for
Use Case 1–2 (S4 Fig), we did not collapse the SRP013822 reads to unique sequences, which
allowed the reconstruction of a Nora_raw_reads genome (S2 File). Finally, the workflow for
Use Case 1–3 (S5 Fig) normalized the abundances of SRP013822 sequence reads using the Gal-
axy tool “Normalize by median” [27] and reconstructed a Nora_Median-Norm-reads genome
(S3 File).
All three reconstructed genomes as well as the Paparazzi-reconstructed JX220408 genome had a high sequence similarity (>96.6% nucleotide identity) with the NC_007919.3 guide genome (S4 File). The final *de novo* (capital letters) assemblies of both the Nora_raw_reads and Nora_Median-Norm-reads genomes entirely covered the JX220408 and NC_007919.3 genomes (both 12333 nt), whereas the *de novo* assembled part of the Nora_MV genome was marginally shorter (12298 nt, the 31 first 5' nucleotides are in lowercase to indicate that they were not *de novo* assembled but instead recovered from the guide genome). To evaluate the quality of assemblies, we remapped the SRP013822 reads to the 3 reconstituted Nora virus genomes as well as to the JX220408 guide genome using the "workflow for remapping in Use Cases 1–1,2,3" (S6 Fig). As can be seen in Fig 2, SRP013822 reads matched the genomes with almost identical profiles and had characteristic size distributions of viral siRNAs with a major peak at 21 nucleotides. Importantly, the numbers of reads re-matched to the Nora virus genomes were 1,578,704 (Nora_MV) > 1,578,135 (Paparazzi—JX220408) > 1,566,909 (Nora_raw_reads) > 1,558,000 (Nora_Median-Norm-reads) > 872,128 (NC_007919.3 reference genome guide).

Thus, Metavisor reconstructed a Nora virus genome Nora_MV whose sequence maximizes the number of vsiRNA read alignments which suggests it is the most accurate genome for the Nora virus present in the datasets. Of note, the Nora_MV genome differs from the JX220408 rNora genome generated by Paparazzi by only two mismatches at positions 367 and 10707, and four 2nt-deletions at positions 223, 365, 9059 and 12217 (see S4 File). These variations did not change the amino acid sequence of the 4 ORFs of the Nora virus. We conclude that Metavisor performs slightly better than Paparazzi for a known virus, using *de novo* assembly of small RNA reads followed by blast-guided scaffolding. We did not observe any benefits of using raw reads or normalized-by-median reads for *de novo* assembly with Oases, but rather a decrease in the accuracy of the reconstructed genome as measured by the number of reads re-mapped to the final genomes (Fig 2).

**Use Case 1–4: detection of multiple viruses in small RNA sequencing datasets.** In order to show the ability of Metavisor in detecting multiple known viruses in small RNA sequencing datasets, we built another workflow Case that performs blastn alignments of Oases contigs on the vir1 reference and reports for all significant alignments without filtering (S7 Fig). Applying this workflow to the SRP013822 sequence datasets produced a list of alignments which contains, as expected, the Nora virus. In addition, contigs were found to align with high significance (Mean BitScore > 200) to the Drosophila A virus and to the Drosophila C virus (S5 File and Table 1), strongly suggesting that the fly stocks analyzed in our previous work were also subject to persistent infection by these viruses [3].

## 2. Discovery of novel viruses

**Use Case 2–1: identification of new viruses in small RNA sequencing datasets.** Using Metavisor, we recently discovered two novel viruses infecting a laboratory colony of *Anopheles coluzzii* mosquitoes [28]. In this case, a workflow (S8 Fig) was used to process small RNA datasets from these mosquitoes (EBI SRA ERP012577) and to assemble a number of Oases contigs that show significant blastx hits with *Dicistroviridae* proteins, including *Drosophila* C virus (DCV) and Cricket paralysis virus (CrPV) proteins (S6 File).

The viral family of *Dicistroviridae* was named from the dicistronic organisation of their genome with a 5' open reading frame encoding a non-structural polyprotein and a second non-overlapping 3' open reading frame encoding the structural polyprotein. In order to construct a potential new *A. coluzzii* dicistrovirus genome, the “Pick Fasta Sequences” tool (S8 Fig) collected blastx hits showing significant alignment with both *Drosophila* C virus and
<figure><img src="image_6.png" /><figcaption>Fig 2. Realignments of small RNA sequence reads to reconstructed (Nora_MV, Nora_raw_reads and Nora_Median-Norm-reads) or published (JX220408.1 and NC_007919.3) Nora virus genomes. Plots (left)</figcaption></figure>
show the abundance of 18-30-nucleotide (nt) small RNA sequence reads matching the genome sequences and histograms (middle) show length distributions of these reads. Positive and negative values correspond to sense and antisense reads, respectively. Total read counts are indicated to the right hand side.

doi:10.1371/journal.pone.0168397.g002

Cricket paralysis viral polyproteins (S7 File) that were further clustered with the "cap3 Sequence Assembly" tool in 4 contigs of 1952, 341, 4688 and 320 nt, respectively (S8 File). These 4 contigs were further aligned to the DCV genome NC_001834.1 sequence with tblastx and integrated in this scaffold sequence with the "blast_to_scaffold" tool to produce a final assembly (S9 File). Re-mapping of the ERP012577 small RNA reads using the "Workflow for remapping in Use Cases 1-1,2,3" (S6 Fig) showed that they mostly align to *de novo* assembled regions (uppercase nucleotides) of this chimeric genome and have a typical size distribution of viral derived siRNA (S9 Fig), suggesting that the NC_001834.1 DCV sequences of the scaffold (lowercase nucleotides) are loosely related to the actual sequence of the novel *A. coluzzii* dicistrovirus. Nevertheless, the composite assembly allowed designing primers in the *de novo* assembled regions to PCR amplify and sequence the regions of the viral genome that could not be *de novo* assembled [28].

Several teams have used siRNA signature (a peak at 21 nt in the size distribution of re-aligned small RNA sequences) as an alternate approach to sequence similarity to identify contigs of potential viral origin [12,13]. In order to further illustrate the flexibility of Metavisor for implementing this strategy, we built a workflow (S10 Fig) to realign ERP012577 small RNA sequences to Oases contigs of length higher than 300 nt and to generate in batch read maps and read size distributions for these contigs using the "Generate readmap and histograms from alignment files" tool (S1 Table). We manually inspected these read maps and size distributions (S10 File) and collected all contigs with a clear siRNA signature (a pick at 21 nt for both forward and reverse strands of contig sequences), 2 sets of contigs with a modest excess of 21 nt reads from the forward strand only and 3 sets of contigs with no siRNA signature as negative controls (S3 Table). With the notable exceptions of loci 3 and 46 contigs, all contigs with a clear siRNA signature blastx aligned to vir1 viral sequences (S3 Table and see here the public Galaxy history for details). Loci 3 and 46 contigs did not align either to the non-redundant protein database of the NCBI and may therefore be of potential viral origin (S3 Table). All 5 negative control contigs with unclear or no RNA signature only aligned significantly to non-viral proteins (S3 Table). Together, these results illustrate the use of Metavisor to implement a sequence-independent strategy based on siRNAs for virus identification [12,13].

**Use Case 2-2: identification of new viruses in mRNA sequencing datasets.** In our study [28], we also used RNAseq libraries from the same *A. coluzzii* colony (EBI-SRA, ERS977505), demonstrating the use of a Metavisor workflow for long RNA sequencing read datasets (S11 Fig). Thus, 100 nt reads were aligned without adapter clipping to the *Anopheles gambiae* genome using bowtie2, and unmatched read were subjected to Oases assembly (kmer range 25 to 69, to take into account longer reads). Oases contigs were then filtered for a size > 5000 nt and aligned to the protein viral reference using blastx. Parsing of blastx alignments with the "blast analysis, by subjects" tool repeatedly pointed to a 8919 nt long Oases contig that matched to structural and non-structural polyproteins of DCV and CrPV (S11 File). This 8919 nt contig (S12 File) completely includes the contigs generated with the small RNA datasets (S8 File) and shows a dicistronic organization which is typical of Dicistroviridae and is referred to as a novel *Anopheles* C Virus [28]. The sequence of this *Anopheles* C Virus is deposited to the NCBI nucleotide database under accession number KU169878. As expected, when realigned to this genome (S12 Fig), the ERP012577 small RNA reads now show a typical alignment profile all along the AnCV genome sequence with a size distribution peaking at the 21 nt length of viral
<table><caption>Table 1. Report table generated by the “Parse blast output and compile hits” tool in Use Case 1–4 showing the presence of *Drosophila A* virus and *Drosophila C* virus in addition to the *Nora* virus in the small RNA sequencing of laboratory *Drosophila*. See Method section for a description of the columns.</caption><thead><tr><th>subject</th><th>subject length</th><th>Total Subject Coverage (nt)</th><th>Relative Subject Coverage</th><th>Best Bit Score</th><th>Mean Bit Score</th></tr></thead><tbody><tr><td>gi|157325505|gb|DQ321720.2|_Nora_virus_complete_genome</td><td>11908</td><td>10211</td><td>0.857</td><td>11840</td><td>4041</td></tr><tr><td>gi|822478532|gb|KP970099.1|_Nora_virus_isolate_RAKMEL13_gp1_(gp1)_gene,_partial_cds;_and_replication_p</td><td>11416</td><td>8736</td><td>0.765</td><td>11441</td><td>3673</td></tr><tr><td>gi|822478537|gb|KP970100.1|_Nora_virus_isolate_GEO58_gp1_(gp1)_gene,_partial_cds;_and_replication_p</td><td>11416</td><td>2463</td><td>0.216</td><td>4028</td><td>3607</td></tr><tr><td>gi|346421290|ref|NC_007919.3|_Nora_virus_complete_genome</td><td>12333</td><td>10530</td><td>0.854</td><td>11809</td><td>2653</td></tr><tr><td>gi|284022350|gb|GQ257737.1|_Nora_virus_isolate_Umea_2007,_complete_genome</td><td>12333</td><td>10530</td><td>0.854</td><td>11809</td><td>2573</td></tr><tr><td>gi|822478527|gb|KP970098.1|_Nora_virus_isolate_AM04_gp1_(gp1)_gene,_partial_cds;_and_replication_po</td><td>11413</td><td>7654</td><td>0.671</td><td>5745</td><td>2489</td></tr><tr><td>gi|822478512|gb|KP970095.1|_Nora_virus_isolate_RAK11_gp1_(gp1)_and_replication_polypeptide_(gp2)_ge</td><td>11416</td><td>6174</td><td>0.541</td><td>5368</td><td>2419</td></tr><tr><td>gi|822478141|gb|KP969947.1|_Drosophila_A_virus_isolate_ywiP_DrosophilaA_RNA-dependent_RNA_polymers</td><td>4516</td><td>4157</td><td>0.921</td><td>6980</td><td>2361</td></tr><tr><td>gi|402295620|gb|JX220408.1|_Nora_virus_isolate_FR1_complete_genome</td><td>12333</td><td>12302</td><td>0.997</td><td>12720</td><td>2324</td></tr><tr><td>gi|822478147|gb|KP969949.1|_Drosophila_A_virus_isolate_delta11_DrosophilaA_RNA-dependent_RNA_polyme</td><td>4481</td><td>4442</td><td>0.991</td><td>7081</td><td>2264</td></tr><tr><td>gi|822478417|gb|KP970078.1|_Nora_virus_isolate_D167_gp1_(gp1),_replication_polypeptide_(gp2),_gp3_</td><td>11895</td><td>7315</td><td>0.615</td><td>5503</td><td>2192</td></tr><tr><td>gi|822478517|gb|KP970096.1|_Nora_virus_isolate_K09_gp1_(gp1)_gene,_partial_cds;_and_replication_pol</td><td>11419</td><td>2027</td><td>0.178</td><td>3490</td><td>2050</td></tr><tr><td>gi|822478144|gb|KP969948.1|_Drosophila_A_virus_isolate_XIB_DrosophilaA_RNA-dependent_RNA_polymerase</td><td>4516</td><td>4507</td><td>0.998</td><td>7092</td><td>2009</td></tr><tr><td>gi|822478150|gb|KP969950.1|_Drosophila_A_virus_isolate_Qdelta_DrosophilaA_RNA-dependent_RNA_polymers</td><td>4476</td><td>4446</td><td>0.993</td><td>7092</td><td>1847</td></tr><tr><td>gi|822478440|gb|KP970082.1|_Nora_virus_isolate_RAKMEL12_gp1_(gp1),_replication_polypeptide_(gp2),_g</td><td>11968</td><td>7214</td><td>0.603</td><td>6396</td><td>1815</td></tr><tr><td>gi|822478497|gb|KP970092.1|_Nora_virus_isolate_delta11_gp1_(gp1)_gene,_partial_cds;_replication_pol</td><td>11157</td><td>2347</td><td>0.210</td><td>3314</td><td>1695</td></tr><tr><td>gi|822478522|gb|KP970097.1|_Nora_virus_isolate_JJ17_gp1_(gp1)_gene,_partial_cds;_and_replication_po</td><td>11420</td><td>993</td><td>0.087</td><td>1674</td><td>1659</td></tr><tr><td>gi|822478482|gb|KP970089.1|_Nora_virus_isolate_IM13_gp1_(gp1)_gene,_partial_cds;_replication_polypeptide</td><td>11103</td><td>1828</td><td>0.165</td><td>1977</td><td>1565</td></tr><tr><td>gi|225356593|gb|FJ150422.1|_Drosophila_A_virus_isolate_HD_complete_genome</td><td>4806</td><td>4753</td><td>0.989</td><td>6902</td><td>1419</td></tr><tr><td>gi|822478430|gb|KP970080.1|_Nora_virus_isolate_MONSIM03_gp1_(gp1),_replication_polypeptide_(gp2),_g</td><td>11968</td><td>2086</td><td>0.174</td><td>2993</td><td>1416</td></tr><tr><td>gi|822478445|gb|KP970083.1|_Nora_virus_isolate_SAF04_gp1_(gp1),_replication_polypeptide_(gp2),_gp3_</td><td>11142</td><td>649</td><td>0.058</td><td>1121</td><td>1121</td></tr><tr><td>gi|822478135|gb|KP969945.1|_Drosophila_A_virus_isolate_XID_DrosophilaA_RNA-dependent_RNA_polymerase</td><td>4516</td><td>4074</td><td>0.902</td><td>3103</td><td>1112</td></tr><tr><td>gi|822478403|gb|KP970076.1|_Nora_virus_isolate_ATH56_gp1_(gp1)_and_replication_polypeptide_(gp2)_ge</td><td>11965</td><td>2344</td><td>0.196</td><td>2812</td><td>1025</td></tr><tr><td>gi|822478412|gb|KP970077.1|_Nora_virus_isolate_IM09_gp1_(gp1),_replication_polypeptide_(gp2),_gp3_</td><td>11965</td><td>3689</td><td>0.308</td><td>2989</td><td>1023</td></tr></tbody></table>

(Continued)
<table><caption>**Table 1. (Continued)**</caption><thead><tr><th>subject</th><th>subject length</th><th>Total Subject Coverage (nt)</th><th>Relative Subject Coverage</th><th>Best Bit Score</th><th>Mean Bit Score</th></tr></thead><tbody><tr><td>gi|822478435|gb|KP970081.1|_Nora_virus_isolate_MON28_gp1_(gp1)_and_replication_polypeptide_(gp2)_gene</td><td>11967</td><td>5859</td><td>0.490</td><td>5445</td><td>1004</td></tr><tr><td>gi|822478507|gb|KP970094.1|_Nora_virus_isolate_K02_gp1_(gp1)_gene,_partial_cds;_replication_polypeptide</td><td>11160</td><td>1289</td><td>0.116</td><td>957</td><td>778</td></tr><tr><td>gi|822478132|gb|KP969944.1|_Drosophila_A_virus_isolate_wipe_Drosophila_A_RNA-dependent_RNA_polymers</td><td>4516</td><td>3070</td><td>0.680</td><td>3097</td><td>742</td></tr><tr><td>gi|822478477|gb|KP970088.1|_Nora_virus_isolate_IM12_gp1_(gp1)_gene,_partial_cds;_replication_polypeptide</td><td>11413</td><td>952</td><td>0.083</td><td>1153</td><td>732</td></tr><tr><td>gi|253761971|ref|NC_012958.1|_Drosophila_A_virus,_complete_genome</td><td>4806</td><td>607</td><td>0.126</td><td>1045</td><td>674</td></tr><tr><td>gi|822478542|gb|KP970101.1|_Nora_virus_isolate_SAFSIM01_gp1_(gp1)_gene,_partial_cds;_and_replication</td><td>11413</td><td>384</td><td>0.034</td><td>661</td><td>661</td></tr><tr><td>gi|822478138|gb|KP969946.1|_Drosophila_A_virus_isolate_LJ35_Drosophila_A_RNA-dependent_RNA_polymers</td><td>4468</td><td>959</td><td>0.215</td><td>848</td><td>502</td></tr><tr><td>gi|9629650|ref|NC_001834.1|_Drosophila_C_virus,_complete_genome</td><td>9264</td><td>6345</td><td>0.685</td><td>1276</td><td>445</td></tr><tr><td>gi|2388672|gb|AF014388.1|_Drosophila_C_virus_strain_EB,_complete_genome</td><td>9264</td><td>6587</td><td>0.711</td><td>1276</td><td>431</td></tr><tr><td>gi|300871949|gb|GU983882.2|_Drosophila_C_virus_isolate_ZW141_polypeptide_nGene,_partial_cds</td><td>500</td><td>272</td><td>0.544</td><td>482</td><td>395</td></tr><tr><td>gi|300871965|gb|GU983892.2|_Drosophila_C_virus_isolate_psjmg_polypeptide_nGene,_partial_cds</td><td>500</td><td>310</td><td>0.620</td><td>491</td><td>353</td></tr><tr><td>gi|300871979|gb|GU983900.2|_Drosophila_C_virus_isolate_AL7_polypeptide_nGene,_partial_cds</td><td>500</td><td>271</td><td>0.542</td><td>489</td><td>342</td></tr><tr><td>gi|300871957|gb|GU983888.2|_Drosophila_C_virus_isolate_Bam73_H_polypeptide_nGene,_partial_cds</td><td>500</td><td>453</td><td>0.906</td><td>491</td><td>323</td></tr><tr><td>gi|300871941|gb|GU983878.2|_Drosophila_C_virus_isolate_mel15_H_polypeptide_nGene,_partial_cds</td><td>500</td><td>453</td><td>0.906</td><td>489</td><td>321</td></tr><tr><td>gi|300871955|gb|GU983885.2|_Drosophila_C_virus_isolate_16a9_polypeptide_nGene,_partial_cds</td><td>490</td><td>151</td><td>0.308</td><td>273</td><td>262</td></tr><tr><td>gi|300871953|gb|GU983884.2|_Drosophila_C_virus_isolate_Tam15_polypeptide_nGene,_partial_cds</td><td>500</td><td>151</td><td>0.302</td><td>273</td><td>262</td></tr></tbody></table>

doi:10.1371/journal.pone.0168397.t001

derived siRNAs and no gap (Fig 3). Of note, we tested in Use Case 2–2 two alternate workflows substituting the Oases assembly tool with Trinity (S13 Fig) and SPAdes (S14 Fig), respectively. Both these workflows were equally able to assemble the genome KU169878 of the Anopheles C Virus (S13 File).

Taken together, the Metavisor Use Cases 2–1 and 2–2 illustrate that when short read datasets do not provide enough sequencing information, an adapted Metavisor workflow (S11 Fig) is able to exploit long reads of RNA sequencing datasets, if available, to assemble a complete viral genome [28].

### 3. Virus detection in human RNAseq libraries

Having illustrated that Metavisor is able to generate robust genome assemblies from known and novel viruses in *Drosophila* and *Anopheles* sequencing datasets, we tested whether it can be used to diagnose viruses in RNA sequencing datasets of human patients from three different studies [29–31].

**Use Case 3–1.** Innate lymphoid cells (ILCs) play a central role in response to viral infection by secreting cytokines crucial for immune regulation, tissue homeostasis, and repair.
Readmaps and size distributions

<figure><img src="image_7.png" /><figcaption>**Fig 3.** Alignments of small RNA sequence reads to the Anopheles C virus genome reconstructed in **Use Case 2–2**. Plot shows the abundance of 18–30-nucleotide (nt) small RNA sequence reads matching the genome sequence and histogram shows the length distribution of these reads. Positive and negative values correspond to sense and antisense reads, respectively.</figcaption></figure>

doi:10.1371/journal.pone.0168397.g003
Therefore, the pathogenic effect of HIV on these cells was recently analyzed in infected or uninfected patients using various approaches, including transcriptome profiling [30]. ILCs are unlikely to be infected *in vivo* by HIV as they lack expression of the CD4 co-receptor of HIV and they are refractory *in vitro* to HIV infection. However, we reasoned that ILCs samples could still be contaminated by infected cells. This might allow Metavisor to detect and assemble HIV genomes from patient's ILC sequencing data (EBI SRP068722).

We imported 40 ICL sequence datasets from the EBI SRP068722 archive and merged the datasets belonging to the same patients. As the data contained short 32 nt reads that in addition had to be 3' trimmed to 27 nt to retain acceptable sequence quality, we designed a workflow for Use Case 3–1 (S15 Fig) that is similar to the workflows used in cases 1–1 and 2–1 for small RNA sequencing data. Thus, the sequencing datasets were depleted from reads aligning to the human genome (hg19) and viral reads were selected by alignment to the NCBI viral sequences using the sRbowtie tool (S1 Table). These reads were further submitted to Oases assembly (kmers 11 to 27, to take into account short reads) and the resulting contigs were aligned to the Nucleotide Viral Blast Database using blastn. Alignments were parsed using the “Parse blast output and compile hits” tool, removing alignments to NCBI sequences related to patents to simplify the report (“Patent” term in the filter option of the “Parse blast output and compile hits” tool). A final report was generated by concatenating the reports produced by this tool for each patient (S14 File and Table 2). In summary, we were able to detect HIV RNAs in samples from 3 out of 4 infected patients whereas all samples from control uninfected patients remained negative for HIV. This Metavisor workflow was able to accurately detect HIV RNA, even in samples where the number of sequence reads was expected to be low.

**Use Case 3–2.** Yozwiak *et al*. searched the presence of viruses in RNA Illumina sequencing data from serums of children suffering from fevers of unknown origins [29]. In this study, paired-end sequencing datasets were depleted from reads aligning to the human genome and the human transcriptome using BLAT and BLASTn, respectively, and the remaining reads were aligned to the NCBI nucleotide database using BLASTn. A virus was considered identified when 10 reads or more aligned to a viral genome which was not tagged as a known lab contaminant.

For a significant number of Patient IDs reported in Table 1 of the article [29], we were not able to find the corresponding sequencing files in the deposited EBI SRP011425 archive. In addition, we did not find the same read counts for these datasets as those indicated by the authors. With these limitations in mind, we downloaded 86 sequencing datasets that could be further concatenated and assigned to 36 patients in Yozwiak *et al* [29]. As sequence reads in SRP011425 datasets are 97 nt long, we adapted a workflow for this Use Case 3–2 (S16 Fig) from the one used in the Use Case 3–1 with the following modifications: (i) sequences reads were depleted from human sequences and viral reads were selected by alignment to the NCBI viral sequences using the Galaxy bowtie2 tool (S1 Table) instead of the sRbowtie tool; (ii) viral reads were submitted to Oases assembly using kmer values ranging from 13 to 69 to take into account long reads; (iii) the SAM file with reads alignments to the vir1 bowtie2 index was parsed using the “join” and “sort” Galaxy tools in order to detect putative false negative datasets with viral reads that fails to produce significant Oases viral contigs.

This workflow generated a report file (S15 File) summarized in Table 3. The results show that Metavisor detected the same viruses as those reported by Yozwiak *et al*. in 17 patients. Although viral reads were detected in 16 other patients, they were not covering sufficient portions of viral genomes to produce significant viral assemblies. Finally, in the three remaining patients (patients 363, 330 and 345 in Table 3), we detected viruses (Dengue virus 2, Stealth virus 1 and Dengue virus 4, respectively) other than those identified by Yozwiak *et al*. These discrepancies are most likely due to misannotation of some of the deposited datasets, which precludes further detailed comparisons.
<table><caption>**Table 2.** HIV detection in ILC patient samples of Use Case 3-1. The table summarizes the report generated by Metavisor from a batch of 40 sequence datasets (S14 File). Metadata associated with each indicated sequence dataset as well as the ability of Metavisor to detect HIV in datasets and patients are indicated.</caption><thead><tr><th># GSM ID</th><th>ID-1</th><th>Patient</th><th>Treatment (SRR annotations)</th><th>HIV status</th><th>Days post HIV infection</th><th>Treat. status</th><th>SRR ID</th><th>Metavisor HIV detection by sample</th><th>number of raw reads</th><th>number of raw reads by patient</th><th>History for Use Case 3-1</th></tr></thead><tbody><tr><td>GSM2043730</td><td>110314</td><td>0450-318</td><td>ILC2</td><td>HIV+</td><td>1</td><td>untreated</td><td>SRR3111582</td><td>+</td><td>7 013 962</td><td>34 252 732</td><td>gil|45357423[gbl]<br>AY535449.1, gi|<br>45357419[gbl|AY535447.1</td></tr><tr><td>GSM2043731</td><td>110314</td><td>0450-318</td><td>ILC3</td><td>HIV+</td><td>1</td><td>untreated</td><td>SRR3111583</td><td>-</td><td>3 246 980</td><td></td><td></td></tr><tr><td>GSM2043732</td><td>180314</td><td>0450-318</td><td>ILC2</td><td>HIV+</td><td>7</td><td>untreated</td><td>SRR3111584</td><td>+</td><td>2 833 634</td><td></td><td></td></tr><tr><td>GSM2043733</td><td>180314</td><td>0450-318</td><td>ILC3</td><td>HIV+</td><td>7</td><td>untreated</td><td>SRR3111585</td><td>+</td><td>2 989 628</td><td></td><td></td></tr><tr><td>GSM2043734</td><td>170414</td><td>0450-318</td><td>ILC2</td><td>HIV+</td><td>38</td><td>untreated</td><td>SRR3111586</td><td>-</td><td>16 248 912</td><td></td><td></td></tr><tr><td>GSM2043735</td><td>170414</td><td>0450-318</td><td>ILC3</td><td>HIV+</td><td>38</td><td>untreated</td><td>SRR3111587</td><td>-</td><td>1 919 616</td><td></td><td></td></tr><tr><td>GSM2043736</td><td>110614</td><td>0387-272</td><td>ILC2</td><td>HIV+</td><td>1</td><td>untreated</td><td>SRR3111588</td><td>-</td><td>60 342 796</td><td>227 307 414</td><td>-</td></tr><tr><td>GSM2043737</td><td>110614</td><td>0387-272</td><td>ILC3</td><td>HIV+</td><td>1</td><td>untreated</td><td>SRR3111589</td><td>-</td><td>34 189 278</td><td></td><td></td></tr><tr><td>GSM2043738</td><td>170614</td><td>0387-272</td><td>ILC2</td><td>HIV+</td><td>7</td><td>untreated</td><td>SRR3111590</td><td>-</td><td>38 030 394</td><td></td><td></td></tr><tr><td>GSM2043739</td><td>170614</td><td>0387-272</td><td>ILC3</td><td>HIV+</td><td>7</td><td>untreated</td><td>SRR3111591</td><td>-</td><td>29 100 534</td><td></td><td></td></tr><tr><td>GSM2043740</td><td>290714</td><td>0387-272</td><td>ILC2</td><td>HIV+</td><td>49</td><td>untreated</td><td>SRR3111592</td><td>-</td><td>43 022 506</td><td></td><td></td></tr><tr><td>GSM2043741</td><td>290714</td><td>0387-272</td><td>ILC3</td><td>HIV+</td><td>49</td><td>untreated</td><td>SRR3111593</td><td>-</td><td>22 621 906</td><td></td><td></td></tr><tr><td>GSM2043742</td><td>41214</td><td>0629-453</td><td>Acute ART + ILC2</td><td>HIV+</td><td>1</td><td>ART</td><td>SRR3111594</td><td>+</td><td>5 061 920</td><td>54 052 098</td><td>gil|296033826[gbl]<br>GU474419.1, gi|<br>269294806[gbl]<br>DM461231.1, gi|<br>269294805[gbl]<br>DM461230.1, gi|<br>296556482[gbl]<br>AF324493.2gi|<br>296556485[gbl|M19921.2,<br>AY535447.1, gi|<br>45357423[gbl|AY535449.1</td></tr><tr><td>GSM2043743</td><td>41214</td><td>0629-453</td><td>Acute ART + ILC3</td><td>HIV+</td><td>1</td><td>ART</td><td>SRR3111595</td><td>-</td><td>8 455 026</td><td></td><td></td></tr><tr><td>GSM2043744</td><td>101214</td><td>0629-453</td><td>Acute ART + ILC2</td><td>HIV+</td><td>6</td><td>ART</td><td>SRR3111596</td><td>-</td><td>12 451 684</td><td></td><td></td></tr><tr><td>GSM2043745</td><td>101214</td><td>0629-453</td><td>Acute ART + ILC3</td><td>HIV+</td><td>6</td><td>ART</td><td>SRR3111597</td><td>+</td><td>6 419 868</td><td></td><td></td></tr><tr><td>GSM2043746</td><td>130115</td><td>0629-453</td><td>Acute ILC2</td><td>HIV+</td><td>40</td><td>ART</td><td>SRR3111598</td><td>+</td><td>6 837 584</td><td></td><td></td></tr><tr><td>GSM2043747</td><td>130115</td><td>0629-453</td><td>Acute ILC3</td><td>HIV+</td><td>40</td><td>ART</td><td>SRR3111599</td><td>+</td><td>14 826 016</td><td></td><td></td></tr><tr><td>GSM2043748</td><td>150714</td><td>0444-312</td><td>3dR10 ILC2</td><td>HIV+</td><td>2</td><td>ART</td><td>SRR3111600</td><td>+</td><td>15 618 282</td><td>39 610 902</td><td>gil|45357423[gbl]<br>AY535449.1, gi|<br>45357419[gbl|AY535447.1</td></tr><tr><td>GSM2043749</td><td>150714</td><td>0444-312</td><td>3dR10 ILC3</td><td>HIV+</td><td>2</td><td>ART</td><td>SRR3111601</td><td>+</td><td>13 491 804</td><td></td><td></td></tr><tr><td>GSM2043750</td><td>220814</td><td>0444-312</td><td>ILC2</td><td>HIV+</td><td>41</td><td>ART</td><td>SRR3111602</td><td>-</td><td>5 259 104</td><td></td><td></td></tr><tr><td>GSM2043751</td><td>220814</td><td>0444-312</td><td>ILC3</td><td>HIV+</td><td>41</td><td>ART</td><td>SRR3111603</td><td>-</td><td>5 241 712</td><td></td><td></td></tr><tr><td>GSM2043752</td><td>10814</td><td>0500-355neg</td><td>ILC2</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111604</td><td>-</td><td>802 632</td><td>11 691 304</td><td>-</td></tr><tr><td>GSM2043753</td><td>10814</td><td>0500-355neg</td><td>ILC3</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111605</td><td>-</td><td>10 888 672</td><td></td><td></td></tr><tr><td>GSM2043754</td><td>80814</td><td>0292-xxxneg</td><td>ILC2</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111606</td><td>-</td><td>5 418 958</td><td>19 222 152</td><td>-</td></tr><tr><td>GSM2043755</td><td>80814</td><td>0292-xxxneg</td><td>ILC3</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111607</td><td>-</td><td>13 803 194</td><td></td><td></td></tr><tr><td>GSM2043756</td><td>90714</td><td>0394-274</td><td>ILC2</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111608</td><td>-</td><td>13 779 570</td><td>15 991 428</td><td>-</td></tr><tr><td>GSM2043757</td><td>90714</td><td>0394-274</td><td>ILC3</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111609</td><td>-</td><td>2 211 858</td><td></td><td></td></tr><tr><td>GSM2043758</td><td>170714</td><td>0218-162neg</td><td>ILC2</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111610</td><td>-</td><td>9 838 776</td><td>18 939 560</td><td>-</td></tr><tr><td>GSM2043759</td><td>170714</td><td>0218-162neg</td><td>ILC3</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111611</td><td>-</td><td>9 100 784</td><td></td><td></td></tr></tbody></table>

(Continued)
<table><caption>Table 2. (Continued)</caption><thead><tr><th># GSM ID</th><th>ID-1</th><th>Patient</th><th>Treatment (SRR annotations)</th><th>HIV status</th><th>Days post HIV infection</th><th>Treat. status</th><th>SRR ID</th><th>Metavisor HIV detection by sample</th><th>number of raw reads</th><th>number of raw reads by patient</th><th>History for Use Case 3-1</th></tr></thead><tbody><tr><td>GSM2043760</td><td>180314</td><td>0311-217HIVneg</td><td>ILC2</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111612</td><td>-</td><td>2 281 560</td><td>7 490 832</td><td>-</td></tr><tr><td>GSM2043761</td><td>180314</td><td>0311-217HIVneg</td><td>ILC3</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111613</td><td>-</td><td>5 209 272</td><td>-</td><td>-</td></tr><tr><td>GSM2043762</td><td>230514</td><td>0440-307neg</td><td>ILC2</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111614</td><td>-</td><td>11 816 186</td><td>21 714 164</td><td>-</td></tr><tr><td>GSM2043763</td><td>230514</td><td>0440-307neg</td><td>ILC3</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111616</td><td>-</td><td>9 897 978</td><td>-</td><td>-</td></tr><tr><td>GSM2043764</td><td>240614</td><td>0518-370neg</td><td>ILC2</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111617</td><td>-</td><td>16 135 602</td><td>16 671 200</td><td>-</td></tr><tr><td>GSM2043765</td><td>240614</td><td>0518-370neg</td><td>ILC3</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111618</td><td>-</td><td>535 598</td><td>-</td><td>-</td></tr><tr><td>GSM2043766</td><td>290714</td><td>0560-420neg</td><td>ILC2</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111619</td><td>-</td><td>1 235 766</td><td>12 912 002</td><td>-</td></tr><tr><td>GSM2043767</td><td>290714</td><td>0560-420neg</td><td>ILC3</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111620</td><td>-</td><td>11 676 236</td><td>-</td><td>-</td></tr><tr><td>GSM2043768</td><td>290714</td><td>0575-419neg</td><td>ILC2</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111621</td><td>-</td><td>8 713 816</td><td>11 833 416</td><td>-</td></tr><tr><td>GSM2043769</td><td>290714</td><td>0575-419neg</td><td>ILC3</td><td>HIV-</td><td>uninfected</td><td>none</td><td>SRR3111622</td><td>-</td><td>3 119 600</td><td>-</td><td>-</td></tr></tbody></table>

doi:10.1371/journal.pone.0168397.t002
<table><caption>**Table 3.** Summary of virus detection in 36 traceable patients of the Use Case 3–2. The Data of this table were extracted from the Metavisor report file available as S15 File. Values of the column "Coverage of complete viral genome (%)" correspond to the fractions (in %) of the complete viral genomes that are covered by blast hits of viral contigs to these genomes and values of the column "Mean blast bit score" correspond to the mean values of the bit scores observed for these blast hits. Note that blast alignments to incomplete viral genomes were not taken into account. For detection of false positives, reads were aligned to the bowtie2 vir1 index before de novo assembly and counts of these reads were reported in the column "Read mapping to vir1 using bowtie2".</caption><thead><tr><th rowspan="2">Patient ID</th><th colspan="3">Extracted from Yozwiak et al. Table 1</th><th colspan="4">Metavisor</th></tr><tr><th># virus reads</th><th># initial reads</th><th>Fraction virus reads</th><th>Yozwiak et al. Virus detection</th><th># reads in NGS datasets</th><th>Metavisor Virus detection</th><th>Coverage of complete viral genome (%)</th><th>Mean blast bit score</th><th>Read mapping to vir1 using bowtie2</th><th>ENA-RUN</th></tr></thead><tbody><tr><td>566</td><td>206</td><td>1.90E+06</td><td>1.08E-04</td><td>Torque teno mini virus 4</td><td>1.07E-06</td><td>none</td><td>-</td><td>-</td><td>No Significant alignments</td><td>SRR453487</td></tr><tr><td>438</td><td>72</td><td>4.40E+06</td><td>1.64E-05</td><td>Human herpesvirus 6</td><td>3.20E-06</td><td>Human herpesvirus 6</td><td>0.36</td><td>239.3</td><td>49 Human_herpesvirus_6</td><td>SRR453437</td></tr><tr><td>401</td><td>2164</td><td>1.80E+05</td><td>1.20E-02</td><td>Hepatitis A virus</td><td>9.94E-05</td><td>Hepatitis A virus</td><td>69.74</td><td>735.2</td><td>6154 reads Hepatitis_A_virus</td><td>SRR453443, SRR453458</td></tr><tr><td>382</td><td>44</td><td>9.60E+05</td><td>4.58E-05</td><td>Human herpesvirus 4</td><td>2.33E-06</td><td>none</td><td>-</td><td>-</td><td>38 Human_herpesvirus_4</td><td>SRR453430</td></tr><tr><td>377</td><td>81</td><td>6.60E+05</td><td>1.23E-05</td><td>Cyclovirus PK5034</td><td>4.91E-06</td><td>Circovirus-like_NI/2007-3</td><td>20.45</td><td>212.0</td><td>58 reads Circovirus-like_NI/2007-376</td><td>SRR453491</td></tr><tr><td>375</td><td>53</td><td>3.80E+06</td><td>1.39E-05</td><td>Porcine circovirus 1</td><td>2.65E-06</td><td>Circovirus-like_NI/2007-3</td><td>20.63</td><td>214.0</td><td>66 reads Circovirus-like_NI/2007-3</td><td>SRR453499</td></tr><tr><td>350</td><td>48</td><td>3.00E+06</td><td>1.60E-05</td><td>Human herpesvirus 6, Torque teno mini virus 4</td><td>1.93E-06</td><td>none</td><td>-</td><td>-</td><td>29 Human_herpesvirus_6128 Torque teno mini virus 4</td><td>SRR453484</td></tr><tr><td>349</td><td>47</td><td>1.60E+06</td><td>2.94E-05</td><td>Torque teno midi virus</td><td>1.03E+06</td><td>none</td><td>-</td><td>-</td><td>No Significant alignments</td><td>SRR453464</td></tr><tr><td>345</td><td>62</td><td>2.20E+06</td><td>2.82E-05</td><td>Beak and feather disease virus</td><td>1.29E+06</td><td>Dengue_virus_4</td><td>0.84</td><td>156.0</td><td>159 Dengue_virus_465 Circovirus-like_NI/2007</td><td>SRR453506</td></tr><tr><td>344</td><td>303</td><td>1.30E+06</td><td>2.33E-04</td><td>Human herpesvirus 6</td><td>7.48E-05</td><td>Human herpesvirus 6</td><td>0.17</td><td>304.0</td><td>184 Human_herpesvirus_6</td><td>SRR453417</td></tr><tr><td>335</td><td>47</td><td>1.70E+06</td><td>2.76E-05</td><td>Torque teno virus</td><td>9.77E-05</td><td>none</td><td>-</td><td>-</td><td>24 Torque teno virus</td><td>SRR453490</td></tr><tr><td>331</td><td>113</td><td>1.40E+06</td><td>8.07E-05</td><td>Torque teno midi virus 2</td><td>8.99E-05</td><td>UNVERIFIED: _Torque_teno_virus_isolate_S55, _complete_genome</td><td>20.93</td><td>353.4</td><td>772 Torque teno virus</td><td>SRR453478</td></tr><tr><td>330</td><td>77</td><td>1.80E+06</td><td>4.28E-05</td><td>TTV-like mini virus</td><td>1.85E-06</td><td>AF191073_Stealth_virus_1_clone_3B43</td><td>2.68</td><td>143.0</td><td>14 reads Stealth_virus_1_clone_C16130_T324 reads Dengue virus 2</td><td>SRR453465, SRR453480</td></tr><tr><td>329</td><td>14</td><td>2.20E+06</td><td>6.36E-06</td><td>Gull circovirus</td><td>2.93E+06</td><td>Circovirus-like_NI/2007-3</td><td>53.64</td><td>339.5</td><td>5 reads Circovirus-like_NI/2007-312 reads Dengue virus 2</td><td>SRR453489, SRR453505</td></tr><tr><td>322</td><td>206</td><td>3.80E+06</td><td>5.42E-05</td><td>Cyclovirus PK5222</td><td>2.63E+06</td><td>Circovirus-like_NI/2007-3</td><td>25.35</td><td>257.0</td><td>208 reads Circovirus-like_NI/2007-3</td><td>SRR453498</td></tr><tr><td>321</td><td>30</td><td>4.60E+06</td><td>6.52E-06</td><td>Porcine circovirus 1</td><td>3.54E+06</td><td>none</td><td>-</td><td>-</td><td>10 reads Circovirus-like_NI/2007-317 reads Dengue virus 2</td><td>SRR453446</td></tr><tr><td>315</td><td>42</td><td>1.90E+06</td><td>2.21E-05</td><td>African swine fever virus</td><td>1.59E+06</td><td>none</td><td>-</td><td>-</td><td>22 reads Dengue virus 2</td><td>SRR453427, SRR453440</td></tr><tr><td>282</td><td>699</td><td>1.60E+06</td><td>4.20E-04</td><td>Dengue virus 2</td><td>9.66E+05</td><td>Dengue virus 2</td><td>3.63</td><td>495.0</td><td>651 reads Dengue virus 2</td><td>SRR453438</td></tr><tr><td>275</td><td>1511</td><td>1.60E+06</td><td>9.70E-04</td><td>Dengue virus 2</td><td>1.11E+06</td><td>Dengue virus 2</td><td>10.23</td><td>539.2</td><td>1436 reads Dengue virus 2</td><td>SRR453450</td></tr></tbody></table>

(Continued)
<table><caption>Table 3. (Continued)</caption><caption>Extracted from Yozwiak et al. Table 1</caption><thead><tr><th rowspan="2">Patient ID</th><th rowspan="2"># virus reads</th><th rowspan="2"># initial reads</th><th rowspan="2">Fraction virus reads</th><th rowspan="2">Yozwiak et al. Virus detection</th><th rowspan="2"># reads in NGS datasets</th><th rowspan="2">Metavisor Virus detection</th><th colspan="3">Metavisor</th><th rowspan="2">Read mapping to vir1 using bowtie2</th><th rowspan="2">ENA-RUN</th></tr><tr><th>Coverage of complete viral genome (%)</th><th>Mean blast bit score</th></tr></thead><tbody><tr><td>274</td><td>27</td><td>1.20E+06</td><td>2.30E-05</td><td>Dengue virus 1</td><td>6.24E+05</td><td>none</td><td>-</td><td>-</td><td>28 reads Dengue virus 1/2</td><td>SRR453460</td></tr><tr><td>270</td><td>28</td><td>1.20E+06</td><td>2.33E-05</td><td>Human herpesvirus 6</td><td>6.79E+05</td><td>none</td><td>-</td><td>-</td><td>20 Human_herpesvirus_6</td><td>SRR453485</td></tr><tr><td>266</td><td>135749</td><td>4.80E+06</td><td>2.80E-02</td><td>Dengue virus 2</td><td>3.36E+06</td><td>Dengue virus 2</td><td>98.65</td><td>1852.3</td><td>121347 reads Dengue virus 2</td><td>SRR453448</td></tr><tr><td>263</td><td>ND</td><td>ND</td><td>ND</td><td>TTV (virochip)</td><td>3142332</td><td>Dengue virus 2</td><td>39.19</td><td>302.0</td><td>75 densovirus92 Dengue virus</td><td>SRR453424, SRR453457</td></tr><tr><td>193</td><td>56</td><td>1.60E+06</td><td>3.50E-05</td><td>Torque teno mini virus 2</td><td>9.13E+05</td><td>none</td><td>-</td><td>-</td><td>6 reads Torque_teno_virus_isolate_TTV-S34</td><td>SRR453510</td></tr><tr><td>187</td><td>4280</td><td>1.10E+06</td><td>3.90E-03</td><td>Dengue virus 2</td><td>5.55E+05</td><td>Dengue virus 2</td><td>39.30</td><td>496.8</td><td>3970 reads Dengue virus 2</td><td>SRR453456</td></tr><tr><td>186</td><td>1701</td><td>2.00E+06</td><td>8.51E-04</td><td>Torque teno virus 15</td><td>1.32E+06</td><td>Torque teno virus (SEN virus)</td><td>55.85</td><td>389.4</td><td>541 SEN virus AY449524.127 reads Torque_teno_virus_15</td><td>SRR453425, SRR453469</td></tr><tr><td>183</td><td>66</td><td>3.00E+06</td><td>2.20E-05</td><td>Human herpesvirus 6</td><td>1.97E+06</td><td>none</td><td>-</td><td>-</td><td>57 Human_herpesvirus_6B</td><td>SRR453481</td></tr><tr><td>180</td><td>42</td><td>8.00E+05</td><td>5.25E-05</td><td>GB virus C</td><td>4.91E+06</td><td>GB virus C</td><td>1.26</td><td>185.0</td><td>41 reads GB virus</td><td>SRR453531</td></tr><tr><td>179</td><td>17</td><td>1.80E+06</td><td>9.44E-06</td><td>Torque teno mini virus 1</td><td>1.31E+06</td><td>none</td><td>-</td><td>-</td><td>No Significant alignments</td><td>SRR453474</td></tr><tr><td>171</td><td>18</td><td>1.20E+06</td><td>1.50E-05</td><td>Torque teno mini virus 2</td><td>7.81E+05</td><td>none</td><td>-</td><td>-</td><td>No Significant alignments</td><td>SRR453509</td></tr><tr><td>168</td><td>ND</td><td>ND</td><td>ND</td><td>TTV (virochip)</td><td>135412</td><td>none</td><td>-</td><td>-</td><td>No Significant alignments</td><td>SRR453451</td></tr><tr><td>161</td><td>14</td><td>3.00E+06</td><td>4.67E-06</td><td>Human parvovirus B19</td><td>2.75E+06</td><td>Human parvovirus B19</td><td>2.89</td><td>242.0</td><td>79 densovirus12 reads Human parvovirus B19</td><td>SRR453495, SRR453504</td></tr><tr><td>159</td><td>143</td><td>2.60E+06</td><td>5.50E-05</td><td>Torque teno mini virus 5</td><td>1.73E+06</td><td>Uncultured virus DNA</td><td>88.14</td><td>339.0</td><td>No Significant alignments</td><td>SRR453500</td></tr><tr><td>156</td><td>213</td><td>2.30E+06</td><td>9.26E-05</td><td>Torque teno midi virus 1</td><td>1.54E+06</td><td>Torque_teno_virus</td><td>31.01</td><td>314.7</td><td>550 reads Torque_teno_virus_isolate_S54</td><td>SRR453493</td></tr><tr><td>131</td><td>24</td><td>1.20E+06</td><td>2.00E-05</td><td>Human herpesvirus 6</td><td>5.46E+05</td><td>Human herpesvirus 6</td><td>0.18</td><td>263.0</td><td>16 reads Human_herpesvirus_6</td><td>SRR453444</td></tr><tr><td>78</td><td>113</td><td>1.20E+06</td><td>9.42E-05</td><td>Human herpesvirus 6</td><td>5.91E+05</td><td>none</td><td>-</td><td>-</td><td>68 reads Human_herpesvirus_6B</td><td>SRR453426</td></tr></tbody></table>

doi:10.1371/journal.pone.0168397.t003
**Use case 3–3.** Matranga *et al*. recently improved library preparation methods for deep sequencing of Lassa and Ebola viral RNAs in clinical and biological samples [31]. Accordingly, they were able to generate sequence datasets of 150 nt reads providing high coverage of the viral genomes. We used these datasets, relevant in the context of Lassa and Ebola outbreak and epidemic response, to demonstrate the versatility of Metavisor as well as its ability to generate high throughput reconstruction of viral genomes.

In order to take into account longer reads and higher viral sequencing depths in the available datasets [31], we adapted a Metavisor workflow for Use Case 3–3 (S17 Fig) as follows: (i) sequencing reads were directly aligned to vir1 sequences using bowtie2, without prior depletion by alignment to the human or rodent hosts; (ii) the Trinity *de novo* assembler [23] that performs well with longer reads was used instead of Oases (S1 Table); (iii) reconstruction of Lassa and Ebola genomes from the sequences of the blast hits with the nucleotide viral blast database was directly performed with the “blast to scaffold” tool without CAP3 assembly since the Trinity contigs were already covering a significant part of the viral genomes; (iv) the reports generated by our “Parse blast output and compile hits” tool as well as the reconstructed genome generated for each sample were merged in single datasets for easier browsing and subsequent phylogenetic or variant analyses; (v) for adaptability of this workflow to any type of virus, we allowed users to specify two input variables at runtime: the name of the virus to be searched for in the analysis and the identifier of the sequence to be used as guide in genome reconstruction steps.

We imported 63 sequence datasets available in the EBI SRA PRJNA254017 and PRJNA 257197 archives [31] and grouped these datasets in Lassa virus (55 fastq files) and Ebola virus (8 fastq files) dataset collections (see Table 4 for description of the sequence datasets). On the one hand, we executed the workflow (S17 Fig) taking the Lassa virus dataset collection as input sequences, “Lassa” as a filter term for the “Parse blast output and compile hits” tool and the NCBI sequence NC_004297.1 as a guide for reconstruction of the Lassa virus segment L. On the other hand, we executed the workflow taking the Ebola virus dataset collection as input sequences, “Ebola” as a filter term for the “Parse blast output and compile hits” tool and the NCBI sequence NC_002549.1 as a guide for reconstruction of the Ebola virus genome.

The results of both analyses are summarized in Table 4. Metavisor was able to detect Ebola virus in all corresponding sequence datasets (S16 File) as well as Lassa virus in 53 out the 55 sequence datasets generated from Lassa virus samples (S17 File). Consistently, Matranga *et al* did not report reconstructed Lassa genomic segments from the two remaining datasets, which likely reflects high read duplication levels in the corresponding libraries [31]. The reconstructed Lassa virus L segments and Ebola virus genomes are compiled in S18 and S19 Files, respectively. In these sequences, *de novo* assembled segments in uppercases are integrated in the reference guide sequence (lowercase) used for the reconstruction. To note, for viruses with segmented genomes such as Lassa virus, the workflow has to be used separately with appropriate guide sequences for the segment to be reconstructed. As an example, we used this workflow with the filter term “Lassa” for the “Parse blast output and compile hits” tool and the Lassa S segment NC_004296.1 for guiding the reconstruction (S20 File).

At this stage, users can use the genomic fasta sequences for further analyses. For instance, multiple sequence alignments can be performed for phylogenetics or variant analyses, or reads in the original datasets can be realigned to the viral genomes to visualize their coverages, as has been done in Use Cases 1 and 2.

## Discussion

Metavisor performs *de novo* assembly of sequencing reads and detects contigs of viral origin through blast alignments, which then can be clustered to reconstruct a viral genome.
<table><caption>Table 4. Summary of detection of Ebola and Lassa viruses in Use Case 3–3. The table summarizes the Metavisor report files available as S16 and S17 Files.</caption><thead><tr><th>Virus</th><th>BioProject</th><th>BioSample id</th><th>SRX number</th><th>SRR number</th><th>Sample ID</th><th>BAM file name</th><th>Source</th><th>Data Type, Selection</th><th>Figure, Table from Matranga et al.</th><th>Metavisor detection (Trinity)</th></tr></thead><tbody><tr><td>EBOV</td><td>PRJNA257197</td><td>SAMN03099684</td><td>SRX733660</td><td>SRR1613381</td><td>G3676-2</td><td>G3676-2_S6_L001_001.bam</td><td>Human</td><td>RNase H</td><td>Figure 5</td><td>+</td></tr><tr><td>EBOV</td><td>PRJNA257197</td><td>SAMN03099684</td><td>SRX733656</td><td>SRR1613377</td><td>G3676-2</td><td>G3676-2-std_S13_L001_001.bam</td><td>Human</td><td>RNase H</td><td>Figure 5</td><td>+</td></tr><tr><td>EBOV</td><td>PRJNA257197</td><td>SAMN03099685</td><td>SRX733661</td><td>SRR1613382</td><td>G3677-1</td><td>G3677-1_S3_L001_001.bam</td><td>Human</td><td>RNase H</td><td>Figure 5</td><td>+</td></tr><tr><td>EBOV</td><td>PRJNA257197</td><td>SAMN03099685</td><td>SRX733657</td><td>SRR1613378</td><td>G3677-1</td><td>G3677-1-std_S10_L001_001.bam</td><td>Human</td><td>RNase H</td><td>Figure 5</td><td>+</td></tr><tr><td>EBOV</td><td>PRJNA257197</td><td>SAMN03099686</td><td>SRX733662</td><td>SRR1613383</td><td>G3677-2</td><td>G3677-2_S2_L001_001.bam</td><td>Human</td><td>RNase H</td><td>Figure 5</td><td>+</td></tr><tr><td>EBOV</td><td>PRJNA257197</td><td>SAMN03099686</td><td>SRX733658</td><td>SRR1613379</td><td>G3677-2</td><td>G3677-2-std_S9_L001_001.bam</td><td>Human</td><td>RNase H</td><td>Figure 5</td><td>+</td></tr><tr><td>EBOV</td><td>PRJNA257197</td><td>SAMN03099687</td><td>SRX733663</td><td>SRR1613384</td><td>G3682-1</td><td>G3682-1_S4_L001_001.bam</td><td>Human</td><td>RNase H</td><td>Figure 5</td><td>+</td></tr><tr><td>EBOV</td><td>PRJNA257197</td><td>SAMN03099687</td><td>SRX733659</td><td>SRR1613380</td><td>G3682-1</td><td>G3682-1-std_S11_L001_001.bam</td><td>Human</td><td>RNase H</td><td>Figure 5</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927412</td><td>SRX719120</td><td>SRR1595772</td><td>G2431</td><td>LASV678_ERCC117</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927412</td><td>SRX719079</td><td>SRR1595696</td><td>G2431</td><td>LASV678_ERCC12</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927488</td><td>SRX719056</td><td>SRR1595665</td><td>ISTH1003</td><td>LASV347_ERCC126</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927488</td><td>SRX718926</td><td>SRR1595500</td><td>ISTH1003</td><td>LASV347_ERCC17</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927485</td><td>SRX718761</td><td>SRR1594619</td><td>ISTH0531</td><td>LASV334_ERCC136</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927485</td><td>SRX719205</td><td>SRR1595943</td><td>ISTH0531</td><td>LASV334_ERCC31</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927498</td><td>SRX719063</td><td>SRR1595673</td><td>ISTH1121</td><td>LASV363_ERCC69</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927498</td><td>SRX719134</td><td>SRR1595797</td><td>ISTH1121</td><td>LASV363_ERCC43</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927498</td><td>SRX719117</td><td>SRR1595763</td><td>ISTH1038</td><td>LASV349_ERCC62</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927489</td><td>SRX718979</td><td>SRR1595558</td><td>ISTH1038</td><td>LASV349_ERCC42</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927510</td><td>SRX718802</td><td>SRR1594664</td><td>ISTH2050</td><td>LASV386_ERCC84</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927503</td><td>SRX719192</td><td>SRR1595909</td><td>ISTH2020</td><td>LASV368_ERCC112</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927484</td><td>SRX718789</td><td>SRR1594651</td><td>ISTH0230</td><td>LASV435_ERCC96</td><td>Human</td><td>RNase H</td><td>Figure 2</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927592</td><td>SRX719159</td><td>SRR1595835</td><td>LM032 Dep</td><td>LM032_Depleted</td><td>Mastomys</td><td>RNase H</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927592</td><td>SRX718836</td><td>SRR1594698</td><td>LM032 Std</td><td>LM032_Standard</td><td>Mastomys</td><td>RNase H</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099734</td><td>SRX733666</td><td>SRR1613388</td><td>NHP_DK9W-AG. dep</td><td>728_Depleted</td><td>Macaque</td><td>RNase H</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099735</td><td>SRX733667</td><td>SRR1613389</td><td>NHP_DK9W-AG. std</td><td>728_Standard</td><td>Macaque</td><td>RNase H</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099736</td><td>SRX733668</td><td>SRR1613390</td><td>NHP_DK9W-AL. dep</td><td>729_Depleted</td><td>Macaque</td><td>RNase H</td><td>Figure 3</td><td>+</td></tr></tbody></table>

(Continued)
<table><caption>Table 4. (Continued)</caption><thead><tr><th>Virus</th><th>BioProject</th><th>BioSample id</th><th>SRX number</th><th>SRR number</th><th>Sample ID</th><th>BAM file name</th><th>Source</th><th>Data Type, Selection</th><th>Figure, Table from Matranga et al.</th><th>Metavisor detection (Trinity)</th></tr></thead><tbody><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099737</td><td>SRX733669</td><td>SRR1613391</td><td>NHP_DK9W-AL.<br>std</td><td>729_Standard</td><td>Macaque</td><td>RNA seq</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099738</td><td>SRX733670</td><td>SRR1613392</td><td>NHP_DK9W-B.<br>dep</td><td>734_Depleted</td><td>Macaque</td><td>RNase H</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099739</td><td>SRX733671</td><td>SRR1613393</td><td>NHP_DK9W-B.<br>std</td><td>734_Standard</td><td>Macaque</td><td>RNA seq</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099740</td><td>SRX733672</td><td>SRR1613394</td><td>NHP_DK9W-K.<br>dep</td><td>733_Depleted</td><td>Macaque</td><td>RNase H</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099741</td><td>SRX733673</td><td>SRR1613395</td><td>NHP_DK9W-K.<br>std</td><td>733_Standard</td><td>Macaque</td><td>RNA seq</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099742</td><td>SRX733674</td><td>SRR1613396</td><td>NHP_DK9W-L.<br>dep</td><td>731_Depleted</td><td>Macaque</td><td>RNase H</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099743</td><td>SRX733675</td><td>SRR1613397</td><td>NHP_DK9W-L.<br>std</td><td>731_Standard</td><td>Macaque</td><td>RNA seq</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099744</td><td>SRX733676</td><td>SRR1613398</td><td>NHP_DK9W-S.<br>dep</td><td>732_Depleted</td><td>Macaque</td><td>RNase H</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099745</td><td>SRX733677</td><td>SRR1613399</td><td>NHP_DK9W-S.<br>std</td><td>732_Standard</td><td>Macaque</td><td>RNA seq</td><td>Figure 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927592</td><td>SRX719168</td><td>SRR1595853</td><td>LM032</td><td>LASV68_BLC</td><td>Mastomys</td><td>RNA seq</td><td>Figure 4,<br>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927476</td><td>SRX727329</td><td>SRR1606288</td><td>G733</td><td>LASV_90</td><td>Human</td><td>RNA seq</td><td>Figure 4,<br>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927592</td><td>SRX733690</td><td>SRR1613412</td><td>LM032</td><td>LM032_HS</td><td>Mastomys</td><td>Hybrid<br>Selection</td><td>Figure 4,<br>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927476</td><td>SRX733681</td><td>SRR1613403</td><td>G733</td><td>G733_HS</td><td>Human</td><td>Hybrid<br>Selection</td><td>Figure 4,<br>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927593</td><td>SRX727318</td><td>SRR1606277</td><td>LM222</td><td>LASV_74</td><td>Mastomys</td><td>RNA seq</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099732</td><td>SRX733664</td><td>SRR1613386</td><td>Z002</td><td>LASV_77</td><td>Mastomys</td><td>RNA seq</td><td>Table 1</td><td>-</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099733</td><td>SRX733665</td><td>SRR1613387</td><td>G090</td><td>LASV_79</td><td>Human</td><td>RNA seq</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927477</td><td>SRX727310</td><td>SRR1606267</td><td>G771</td><td>LASV94</td><td>Human</td><td>RNA seq</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927399</td><td>SRX734464</td><td>SRR1614275</td><td>G2230</td><td>Solexa-100929.<br>tagged_332</td><td>Human</td><td>RNA seq</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927483</td><td>SRX731079</td><td>SRR1610580</td><td>ISTH0073</td><td>Solexa-106870.<br>tagged_851</td><td>Human</td><td>RNA seq</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927500</td><td>SRX719163</td><td>SRR1595846</td><td>ISTH1137</td><td>LASV353_BLC</td><td>Human</td><td>RNA seq</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927503</td><td>SRX718749</td><td>SRR1594606</td><td>ISTH2020</td><td>LASV368_ERCC03</td><td>Human</td><td>RNA seq</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927504</td><td>SRX727274</td><td>SRR1606236</td><td>ISTH2025</td><td>LASV374_ERCC58</td><td>Human</td><td>RNA seq</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927510</td><td>SRX718860</td><td>SRR1594723</td><td>ISTH2050</td><td>LASV386_ERCC48</td><td>Human</td><td>RNA seq</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927484</td><td>SRX718809</td><td>SRR1594671</td><td>ISTH0230</td><td>LASV435_ERCC53</td><td>Human</td><td>RNA seq</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927593</td><td>SRX733692</td><td>SRR1613414</td><td>LM222</td><td>LM222_HS</td><td>Mastomys</td><td>Hybrid<br>Selection</td><td>Table 1</td><td>+</td></tr></tbody></table>

(Continued)
<table><caption>Table 4. (Continued)</caption><thead><tr><th>Virus</th><th>BioProject</th><th>BioSample id</th><th>SRX number</th><th>SRR number</th><th>Sample ID</th><th>BAM file name</th><th>Source</th><th>Data Type, Selection</th><th>Figure, Table from Matranga et al.</th><th>Metavisor detection (Trinity)</th></tr></thead><tbody><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099732</td><td>SRX733678</td><td>SRR1613400</td><td>Z002</td><td>Z002_HS</td><td>Mastomys</td><td>Hybrid Selection</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN03099733</td><td>SRX733679</td><td>SRR1613401</td><td>G090</td><td>G090_HS</td><td>Human</td><td>Hybrid Selection</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927477</td><td>SRX733682</td><td>SRR1613404</td><td>G771</td><td>G771_HS</td><td>Human</td><td>Hybrid Selection</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927399</td><td>SRX733680</td><td>SRR1613402</td><td>G2230</td><td>G2230_HS</td><td>Human</td><td>Hybrid Selection</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927483</td><td>SRX733683</td><td>SRR1613405</td><td>ISTH0073</td><td>ISTH0073_HS</td><td>Human</td><td>Hybrid Selection</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927500</td><td>SRX733685</td><td>SRR1613407</td><td>ISTH1137</td><td>ISTH1137_HS</td><td>Human</td><td>Hybrid Selection</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927503</td><td>SRX733686</td><td>SRR1613408</td><td>ISTH2020</td><td>ISTH2020_HS</td><td>Human</td><td>Hybrid Selection</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927504</td><td>SRX733687</td><td>SRR1613409</td><td>ISTH2025</td><td>ISTH2025_HS</td><td>Human</td><td>Hybrid Selection</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927510</td><td>SRX733688</td><td>SRR1613410</td><td>ISTH2050</td><td>ISTH2050_HS</td><td>Human</td><td>Hybrid Selection</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927484</td><td>SRX733684</td><td>SRR1613406</td><td>ISTH0230</td><td>ISTH0230_HS</td><td>Human</td><td>Hybrid Selection</td><td>Table 1</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927592</td><td>SRX733689</td><td>SRR1613411</td><td>LM032</td><td>LM032_Depleted</td><td>Mastomys</td><td>cDNA</td><td>ND, manually added to the original sup file 3</td><td>+</td></tr><tr><td>LASV</td><td>PRJNA254017</td><td>SAMN02927592</td><td>SRX733691</td><td>SRR1613413</td><td>LM032</td><td>LM032_Standard</td><td>Mastomys</td><td>cDNA</td><td>ND, manually added to the original sup file 3</td><td>-</td></tr></tbody></table>

doi:10.1371/journal.pone.0168397.t004
On the one hand, this strategy reduces the rate of false positives since the ability to form contigs that align to known viral sequences is a strong evidence of the presence of a full viral genome in the analyzed datasets. In addition, we advise Metavisor users to remove sequence reads that align to genomes of hosts, symbionts or parasites, if these are known and available (see Fig 1). Although this treatment can be skipped (as in Use Case 3–3), it avoids chimeric assemblies of viral and nonviral sequences, while speeding up the assembly of contigs of potential viral origin. It also ensures that sequences of the host genome that have been annotated as Endogenous Viral Elements (EVEs) are not retained for viral contigs assembly. Users should keep in mind that EVEs that have not yet be identified as such may be retained by Metavisor as potential viral contigs. Should this happen, Webster et al. [13] have demonstrated that, when available, mapping of host genome sequencing reads to these contigs allows to discriminate between EVE and virus sequences. As illustrated in Use Case 2–1, when a host is known for having antiviral RNAi pathways, re-mapping small RNA reads and plotting their length distribution can also add support to the infectious origin of candidate viral contigs (21nt read peak, sense and antisense reads aligning along the contig).

On the other hand, Metavisor workflows may fail to assemble viral contigs when the abundance of viral reads is too low in sequenced samples, or when these reads align to short, scattered regions of viral genomes. However, as illustrated in the Use Case 3–2, it is possible to keep track of these false negatives by implementing a workflow that annotates and counts viral reads before the *de novo* contig assembly steps.

We developed Metavisor in Galaxy in order to benefit from a well supported framework allowing execution of computational tools and workflows through a user-friendly web interface. In addition, the advanced Galaxy functionalities ensure the highest levels of computational analyses, through rigorous recording of the produced data and metadata and of the used parameters as well as the ability to share, publish and reproduce these analyses (see Metavisor availability section in Experimental Procedures). Another major benefit is that, as any Galaxy workflow, the Metavisor workflows may be modified or extended by users. If they are already available in a Galaxy tool shed, integration of new tools in a workflow is straightforward, thanks to the Galaxy workflow editor. Although it requires coding skills, any other freely available software can be adapted to the Galaxy framework and used in a Metavisor workflow.

Through use cases, we have shown that Metavisor is adaptable: short or longer reads from small RNAseq, RNAseq or DNAseq can be used as input data with or without adapter clipping; read datasets can be used as is, or compressed using reads-to-sequences or normalization by median procedures [27]; a variety of alignment and *de novo* assembly tools can be used, provided that they have been adapted for their execution in the Galaxy framework; finally, although we provide the vir1 nucleotide and protein references to identify sequences of viral origin users are free to upload and work with their own viral references. Thus, Metavisor provides biologists and clinicians with an accessible framework for detection, reconstruction or discovery of viruses.

Viral sequences reconstructed by Metavisor can be used in a large range of subsequent analyses, including phylogenetic or genetic drift analyses in contexts of epidemics or virus surveillance in field insect vectors, animal or human populations, or systematic identification of viruses for evaluation of their morbidity. In Use Cases 3–1 to 3–3, we have shown that Metavisor allows analysis of numerous datasets in batch with consistent tracking of individual samples. Thus, we are confident that Metavisor is scalable to large epidemiological studies or to clinical diagnosis in hospital environments. For instance, it could be used to analyse RNAseq data from Zika infected patients [32,33]. Finally, we wish to stress that Metavisor has the potential for integrating detection or diagnosis of non-viral, microbial components in biological samples.
Eukaryotic parasites or symbionts and bacteria are mostly detectable in sequencing datasets from their abundant ribosomal RNAs whose sequences are strongly conserved in the main kingdoms. This raises specific issues for their accurate identification and their taxonomic resolution that are not currently addressed by Metavisor. However, many tools and databases [34] addressing these metagenomics challenges can be adapted, when not already, to the Galaxy framework. For instance, Qiime [35] and the SILVA database of ribosomal RNAs [36] can be used within Galaxy and could thus be integrated in future Metavisor workflows aiming at detection and discovery of non virus organisms in deep sequence datasets.

## Supporting Information

**S1 Fig.** Screenshot of the "Retrieve FASTA from NCBI" tool from NCBI to retrieve viral nucleotide (A) or protein (B) vir1 sequences. The query string "txid10239[orgn] NOT txid131567[orgn] NOT phage" retrieves viruses sequences (txid10239) while filtering out cellular organisms sequences (txid131567) and phage sequences.
(PDF)

**S2 Fig.** Screenshot of an output produced by the "Parse blast output and compile hits" Metavisor tool.
(PDF)

**S3 Fig.** Screenshot of Metavisor workflow for Use Case 1-1.
(PDF)

**S4 Fig.** Screenshot of Metavisor workflow for Use Case 1-2.
(PDF)

**S5 Fig.** Screenshot of Metavisor workflow for Use Case 1-3.
(PDF)

**S6 Fig.** Screenshot of Metavisor workflow for remapping for use cases 1-1, 1-2, 1-3.
(PDF)

**S7 Fig.** Screenshot of Metavisor workflow for Use Case 1-4.
(PDF)

**S8 Fig.** Screenshot of Metavisor workflow for Use Case 2-1.
(PDF)

**S9 Fig.** Alignments of small RNA sequence reads to the partially reconstructed Anopheles C virus genome (Use Case 2-1). Plot shows the abundance of 18-30-nucleotide (nt) small RNA sequence reads matching the genome sequences and histogram shows length distributions of these reads. Positive and negative values correspond to sense and antisense reads, respectively.
(PDF)

**S10 Fig.** Screenshot of Metavisor workflow for small RNA profiling of contigs.
(PDF)

**S11 Fig.** Screenshot of Metavisor workflow for Use Case 2-2.
(PDF)

**S12 Fig.** Screenshot of Metavisor workflow for remapping in use cases 2-1 and 2-2.
(PDF)
S13 Fig. Screenshot of Metavisitor workflow for Trinity test in Use Case 2-2.
(PDF)

S14 Fig. Screenshot of Metavisitor workflow for SPAdes test in Use Case 2-2.
(PDF)

S15 Fig. Screenshot of Metavisitor workflow for Use Case 3-1.
(PDF)

S16 Fig. Screenshot of Metavisitor workflow for Use Case 3-2.
(PDF)

S17 Fig. Screenshot of Metavisitor workflow for Use Case 3-3.
(PDF)

S1 File. Nora_MV sequence of Nora virus reconstructed by Metavisitor using reads col- lapsed to unique sequences (Use Case 1-1).
(TXT)

S2 File. Nora_raw_reads sequence of Nora virus reconstructed by Metavisitor using raw reads (Use Case 1-2).
(TXT)

S3 File. Nora_Median-Norm-reads sequence of Nora virus reconstructed by Metavisitor using normalisation of read abundance by median procedure (Use Case 1-3).
(TXT)

S4 File. MAFFT (http://www.ebi.ac.uk/Tools/msa/mafft/) Multiple Alignment of the Nora virus genome sequences published (JX220408.1 and NC_007919.3) or generated in Use Cases 1-1 to 1-3 (Nora_MV, Nora_raw_reads and Nora_Median-Norm-reads). A view of the alignments was produced by MView (http://www.ebi.ac.uk/Tools/msa/mview/). The html file can be visualized by opening it locally with a web browser.
(HTML)

S5 File. Output of the "parse blast output and compile hits" tool in Use Case 1-4.
(TXT)

S6 File. Output of the "parse blast output and compile hits" tool in Use Case 2-1.
(TXT)

S7 File. Output of the "Pick Fasta Sequences" tool in Use Case 2-1.
(TXT)

S8 File. Sequences of the 4 contigs generated by the "CAP3 sequence assembly" tool in Use Case 2-1.
(TXT)

S9 File. Integration of the 4 assembled contigs (S8 File) in the DCV genome scaffold NC_001834.1 by the "blast_to_scaffold" tool. Lowercase correspond to NC_001834.1 sequences while uppercase correspond to contig sequences.
(TXT)

S10 File. siRNA profiling of *de novo* assembled contigs in Use Case 2-1. Small RNA sequences reads were aligned to the contigs and size distribution and read maps were generated using the "Generate readmap and histograms from alignment files" tool. Plots show the map and abundance of 18-30 nt small RNA reads for indicated contigs and histograms show
length distributions of these reads. Positive and negative values correspond to sense and anti-sense reads, respectively.

(PDF)

S11 File. Parsing of blastx alignments with the "blast analysis, by subjects" tool in Use Case 2-2.

(TXT)

S12 File. Sequence of the 8919 nt contig in Use Case 2-2.

(TXT)

S13 File. MAFFT (http://www.ebi.ac.uk/Tools/msa/mafft/) Multiple Alignment in Clustal format of 3 AnCV genomes reconstructed with Oases, Trinity and SPAdes assembly programs.

(TXT)

S14 File. Merge of all reports generated by the "Parse blast output and compile hits" tool in Use Case 3-1.

(TXT)

S15 File. Merge of all reports generated by the "Parse blast output and compile hits" tool in Use Case 3-2. These reports are summarized in Table 3.

(TXT)

S16 File. Merge of all reports for Ebola virus generated by the "Parse blast output and compile hits" tool in Use Case 3-3. These reports are summarized in Table 4.

(TXT)

S17 File. Merge of all reports for Lassa virus generated by the "Parse blast output and compile hits" tool in Use Case 3-3. These reports are summarized in Table 4.

(TXT)

S18 File. Lassa virus segment L reconstructed sequences in NC_004297.1 scaffold in Use Case 3-3.

(TXT)

S19 File. Ebola virus reconstructed sequences in NC_002549.1 scaffold in Use Case 3-3.

(TXT)

S20 File. Lassa virus segment S reconstructed sequences in NC_004296.1 scaffold in Use Case 3-3.

(TXT)

S1 Table. Galaxy tools used in Metavisor.

(PDF)

S2 Table. Duration of execution of the Metavisor workflows. The times given correspond to execution of the workflows on a 16-core (2GHz) machine with 96 Mo RAM, Galaxy release 16.04.

(PDF)

S3 Table. Sequence-independent strategy to identify candidate viral contigs (Use Case 2-1). Set of contigs (Loci) with clear (+), unclear (?) or no siRNA signature were manually selected from S10 Fig and tested for significant blastx alignment against the vir1 index and the Non-redundant NCBI protein database (october 201).

(PDF)
**Acknowledgments**

We thank the Galaxy community for their support, Juliana Pegoraro, Eugeni Belda and Emmanuel Bischoff for helpful discussions and Julie Reveillaud for critical reading of the manuscript. GC, MvdB and CA conceived the project. CA and MvdB developed and implemented tools in the Galaxy framework. GC, MvdB and CA performed bioinformatics analysis. GC, MvdB, KV and CA wrote the manuscript. CA and KV provided funding. All authors read and approved the final manuscript.

**Author Contributions**

**Conceptualization:** GC MvdB CA.

**Data curation:** MvdB CA.

**Formal analysis:** GC MvdB CA.

**Funding acquisition:** KV CA.

**Investigation:** GC CA.

**Methodology:** GC MvdB CA.

**Project administration:** CA.

**Resources:** GC MvdB CA.

**Software:** MvdB CA.

**Supervision:** CA.

**Validation:** GC CA.

**Visualization:** MvdB CA.

**Writing – original draft:** GC MvdB KV CA.

**Writing – review & editing:** GC KV CA.

**References**

1. Kingsolver MB, Huang Z, Hardy RW. Insect antiviral innate immunity: pathways, effectors, and connections. J Mol Biol. Elsevier Ltd; 2013; 425: 4921–4936. doi: 10.1016/j.jmb.2013.10.006 PMID: 24120681

2. Ding S-W, Voinnet O. Antiviral immunity directed by small RNAs. Cell. Elsevier Inc.; 2007; 130: 413–426. doi: 10.1016/j.cell.2007.07.039 PMID: 17693253

3. van Mierlo JT, Bronkhorst AW, Overheul GJ, Sadanandan SA, Ekström J-O, Heestermans M, et al. Convergent evolution of argonaute-2 slicer antagonism in two distinct insect RNA viruses. Schneider DS, editor. PLoS Pathog. Public Library of Science; 2012; 8: e1002872. doi: 10.1371/journal.ppat.1002872 PMID: 22916019

4. Vodovar N, Goic B, Blanc H, Saleh M-C. In silico reconstruction of viral genomes from small RNAs improves virus-derived small interfering RNA profiling. J Virol. 2011; 85: 11016–11021. doi: 10.1128/JVI.05647-11 PMID: 21880776

5. Warren RL, Sutton GG, Jones SJM, Holt RA. Assembling millions of short DNA sequences using SSAKE. Bioinformatics. 2007; 23: 500–501. doi: 10.1093/bioinformatics/btl629 PMID: 17158514

6. de Andrade RRS, Vaslin MFS. SearchSmallRNA: a graphical interface tool for the assemblage of viral genomes using small RNA libraries data. Virol J. 2014; 11: 45. doi: 10.1186/1743-422X-11-45 PMID: 24607237

7. Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008; 18: 821–829. doi: 10.1101/gr.074492.107 PMID: 18349386
8. Kreuze JF, Perez A, Untiveros M, Quispe D, Fuentes S, Barker I, et al. Complete viral genome sequence and discovery of novel viruses by deep sequencing of small RNAs: a generic method for diagnosis, discovery and sequencing of viruses. Virology. 2009; 388: 1-7. doi: 10.1016/j.virol.2009.03.024 PMID: 19394993

9. Wu Q, Luo Y, Lu R, Lau N, Lai EC, Li W-X, et al. Virus discovery by deep sequencing and assembly of virus-derived small silencing RNAs. Proc Natl Acad Sci U S A. 2010; 107: 1606-1611. doi: 10.1073/pnas.0911353107 PMID: 20080648

10. Seguin J, Rajeswaran R, Malpica-López N, Martin RR, Kasschau K, Dolja VV, et al. De novo reconstruction of consensus master genomes of plant RNA and DNA viruses from siRNAs. Pappu H, editor. PLoS One. Public Library of Science; 2014; 9: e88513. doi: 10.1371/journal.pone.0088513 PMID: 24523907

11. Ho T, Tzanetakis IE. Development of a virus detection and discovery pipeline using next generation sequencing. Virology. 2014; 471-473: 54-60. doi: 10.1016/j.virol.2014.09.019 PMID: 25461531

12. Aguiar ERGR, Olmo RP, Paro S, Ferreira FV, de Faria IJ da S, Todjro YMH, et al. Sequence-independent characterization of viruses based on the pattern of viral small RNAs produced by the host. Nucleic Acids Res. 2015

13. Webster CL, Waldron FM, Robertson S, Crowson D, Ferrari G, Quintana JF, et al. The Discovery, Distribution, and Evolution of Viruses Associated with Drosophila melanogaster. PLoS Biol. 2015; 13: e1002210. doi: 10.1371/journal.pbio.1002210 PMID: 26172158

14. Surget-Groba Y, Montoya-Burgos JI. Optimization of de novo transcriptome assembly from next-generation sequencing data. Genome Res. Cold Spring Harbor Lab; 2010; 20: 1432-1440. doi: 10.1101/gr.103846.109 PMID: 20693479

15. Goecks J, Nekrutenko A, Taylor J, Galaxy Team. Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. BioMed Central Ltd; 2010; 11: R86. doi: 10.1186/gb-2010-11-8-r86 PMID: 20738864

16. Blankenberg D, Von Kuster G, Corao N, Ananda G, Lazarus R, Mangan M, et al. Galaxy: a web-based genome analysis tool for experimentalists. Curr Protoc Mol Biol. Hoboken, NJ, USA: John Wiley & Sons, Inc.; 2010; Chapter 19: Unit 19.10.1-21.

17. Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, et al. GenBank. Nucleic Acids Res. 2013; 41: D36-42. doi: 10.1093/nar/gks1195 PMID: 23193287

18. Cock PJA, Chilton JM, Grüning B, Johnson JE, Soranzo N. NCBI BLAST+ integrated into Galaxy. GigaScience. 2015; 4: 39. doi: 10.1186/s13742-015-0080-7 PMID: 26336600

19. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009; 10: R25. doi: 10.1186/gb-2009-10-3-r25 PMID: 19261174

20. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012; 9: 357-359. doi: 10.1038/nmeth.1923 PMID: 22388286

21. Katzourakis A, Gifford RJ. Endogenous viral elements in animal genomes. PLoS Genet. 2010; 6: e1001191. doi: 10.1371/journal.pgen.1001191 PMID: 21124940

22. Schulz MH, Zerbino DR, Vingron M, Birney E. Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012; 28: 1086-1092. doi: 10.1093/bioinformatics/bts094 PMID: 22368243

23. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011; 29: 644-652. doi: 10.1038/nbt.1883 PMID: 21572440

24. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012; 19: 455-477. doi: 10.1089/cmb.2012.0021 PMID: 22506599

25. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009; 10: 421. doi: 10.1186/1471-2105-10-421 PMID: 20003500

26. Huang X, Madan A. CAP3: A DNA sequence assembly program. Genome Res. 1999; 9: 868-877. doi: 10.1038/nmeth.1923 PMID: 22388286

27. Titus Brown C, Howe A, Zhang Q, Pyrkosz AB, Brom TH. A Reference-Free Algorithm for Computational Normalization of Shotgun Sequencing Data [Internet]. arXiv [q-bio.GN]. 2012. Available: http://arxiv.org/abs/1203.4802

28. Carissimo G, Eiglmeier K, Reveillaud J, Holm I, Diallo M, Diallo D, et al. Identification and Characterization of Two Novel RNA Viruses from Anopheles gambiae Species Complex Mosquitoes. PLoS One. 2016; 11: e0153881. doi: 10.1371/journal.pone.0153881 PMID: 27138938
29. Yozwiak NL, Skewes-Cox P, Stenglein MD, Balmaseda A, Harris E, DeRisi JL. Virus identification in unknown tropical febrile illness cases using deep sequencing. PLoS Negl Trop Dis. 2012; 6: e1485. doi: 10.1371/journal.pntd.0001485 PMID: 22347512

30. Kløverpris HN, Kazer SW, Mjösberg J, Mabuka JM, Wellmann A, Ndhlovu Z, et al. Innate Lymphoid Cells Are Depleted Irreversibly during Acute HIV-1 Infection in the Absence of Viral Suppression. Immunity. 2016; 44: 391–405. doi: 10.1016/j.immuni.2016.01.006 PMID: 26850658

31. Matranga CB, Andersen KG, Winnicki S, Busby M, Gladden AD, Tewhey R, et al. Enhanced methods for unbiased deep sequencing of Lassa and Ebola RNA viruses from clinical and biological samples. Genome Biol. 2014; 15: 1–12.

32. Sardi SI, Somasekar S, Naccache SN, Bandeira AC, Tauro LB, Campos GS, et al. Coinfections of Zika and Chikungunya Viruses in Bahia, Brazil, Identified by Metagenomic Next-Generation Sequencing. J Clin Microbiol. 2016; 54: 2348–2353. doi: 10.1128/JCM.00877-16 PMID: 27413190

33. Wang Z, Ma'ayan A. An open RNA-Seq data analysis pipeline tutorial with an example of reprocessing data from a recent Zika virus study. F1000Res. 2016; 5: 1574. doi: 10.12688/f1000research.9110.1 PMID: 27583132

34. Kim M, Lee K-H, Yoon S-W, Kim B-S, Chun J, Yi H. Analytical tools and databases for metagenomics in the next-generation sequencing era. Genomics Inform. 2013; 11: 102–113. doi: 10.5808/GI.2013.11.3. 102 PMID: 24124405

35. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010; 7: 335–336. doi: 10.1038/ nmeth.f.303 PMID: 20383131

36. Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013; 41: D590–6. doi: 10.1093/nar/gks1219 PMID: 23193283