Loading...
Evaluating assembly and variant calling software for strain-resolved analysis of large DNA viruses
Deng, Zhi-Luo Dhingra, Akshay Fritz, Adrian Götting, Jasper Münch, Philipp C Steinbrück, Lars Schulz, Thomas F Ganzenmüller, Tina McHardy, Alice C
Deng, Zhi-Luo
Dhingra, Akshay
Fritz, Adrian
Götting, Jasper
Münch, Philipp C
Steinbrück, Lars
Schulz, Thomas F
Ganzenmüller, Tina
McHardy, Alice C
Citations
Altmetric:
Advisors
Editors
Other Contributors
Issue Date
2020-07-07
Submitted date
Files
Loading...
Open Access publication
Adobe PDF, 4.34 MB
Loading...
supplemental material
Adobe PDF, 3.31 MB
Other Titles
Abstract
Infection with human cytomegalovirus (HCMV) can cause severe complications in immunocompromised individuals and congenitally infected children. Characterizing heterogeneous viral populations and their evolution by high-throughput sequencing of clinical specimens requires the accurate assembly of individual strains or sequence variants and suitable variant calling methods. However, the performance of most methods has not been assessed for populations composed of low divergent viral strains with large genomes, such as HCMV. In an extensive benchmarking study, we evaluated 15 assemblers and 6 variant callers on 10 lab-generated benchmark data sets created with two different library preparation protocols, to identify best practices and challenges for analyzing such data. Most assemblers, especially metaSPAdes and IVA, performed well across a range of metrics in recovering abundant strains. However, only one, Savage, recovered low abundant strains and in a highly fragmented manner. Two variant callers, LoFreq and VarScan2, excelled across all strain abundances. Both shared a large fraction of false positive variant calls, which were strongly enriched in T to G changes in a ‘G.G’ context. The magnitude of this context-dependent systematic error is linked to the experimental protocol. We provide all benchmarking data, results and the entire benchmarking workflow named QuasiModo, Quasispecies Metric determination on omics, under the GNU General Public License v3.0 (https://github.com/hzi-bifo/Quasimodo), to enable full reproducibility and further benchmarking on these and other data.
Citation
Briefings in Bioinformatics, 2020;, bbaa123, https://doi.org/10.1093/bib/bbaa123.
Publisher
Journal
PubMed ID
PubMed Central ID
Additional Links
Embedded video
Type
Problem solving protocol
Language
en
Description
Series/Report no.
ISSN
EISSN
1477-4054
ISBN
ISMN
Gov't Doc #
Sponsors
Deutsches Zentrum für Infektionsforschung
License
Attribution 4.0 International