UnCoVar: a reproducible and scalable workflow for transparent and robust virus variant calling and lineage assignment using SARS-CoV-2 as an example.

BACKGROUND: At a global scale, the SARS-CoV-2 virus did not remain in its initial genotype for a long period of time, with the first global reports of variants of concern (VOCs) in late 2020. Subsequently, genome sequencing has become an indispensable tool for characterizing the ongoing pandemic, particularly for typing SARS-CoV-2 samples obtained from patients or environmental surveillance. For such SARS-CoV-2 typing, various in vitro and in silico workflows exist, yet to date, no systematic cross-platform validation has been reported. RESULTS: In this work, we present the first comprehensive cross-platform evaluation and validation of in silico SARS-CoV-2 typing workflows. The evaluation relies on a dataset of 54 patient-derived samples sequenced with several different in vitro approaches on all relevant state-of-the-art sequencing platforms. Moreover, we present UnCoVar, a robust, production-grade reproducible SARS-CoV-2 typing workflow that outperforms all other tested approaches in terms of precision and recall. CONCLUSIONS: In many ways, the SARS-CoV-2 pandemic has accelerated the development of techniques and analytical approaches. We believe that this can serve as a blueprint for dealing with future pandemics. Accordingly, UnCoVar is easily generalizable towards other viral pathogens and future pandemics. The fully automated workflow assembles virus genomes from patient samples, identifies existing lineages, and provides high-resolution insights into individual mutations. UnCoVar includes extensive quality control and automatically generates interactive visual reports. UnCoVar is implemented as a Snakemake workflow. The open-source code is available under a BSD 2-clause license at github.com/IKIM-Essen/uncovar.

Phylobook: a tool for display, clade annotation and extraction of sequences from molecular phylogenies.

As the volume of sequence data from variable pathogens increases, means of analyzing, annotating and extracting specific taxa for study becomes more difficult. To meet these challenges for datasets with hundreds to thousands of taxa, 'Phylobook' was developed. Starting with a sequence alignment file, Phylobook generates and displays phylogenetic trees adjacent to highlighter plots showing the position of mutations, and allows the user to identify lineages and recombinants, annotate and export selected subsets of sequences for downstream analysis. Accurate lineage assignment, which is difficult to automate, is aided using annotations created by different clustering methods. Phylobook provides web-based display combined with automated clustering and manual editing to allow for expert assessment and correction of lineage assignments and extraction for downstream analysis.

Phylobook generates and displays phylogenetic trees adjacent to highlighter plots showing the position of mutations, allows the user to identify lineages and recombinants, annotate and export selected subsets of sequences for downstream analysis.

Molecular phylogeny and evolutionary dynamics of influenza A nonstructural (NS) gene.

While the nonstructural gene (NS) of the influenza A virus plays a crucial role in viral virulence and replication, the complete understanding of its molecular phylogeny and evolutionary dynamics remains lacking. In this study, the phylogenetic analysis of 7581 NS sequences revealed ten distinct lineages within alleles A and B: three host-specific (human, classical swine and equine), two reassortment-originated (A(H1N1)pdm09 and triple reassortment swine), one transmission-originated (Eurasian swine), and two geographically isolated avian (Eurasian/Oceanian and North American) for allele A and two geographically isolated avian (Eurasian/Oceanian and North American) for allele B. The average nucleotide substitution rates of the lineages range from 1.24×10(-3) (equine) to 4.34×10(-3) (A(H1N1)pdm09) substitutions per site per year. The selection pressure analysis demonstrated that the dN/dS ratio of the NS gene in A(H1N1)pdm09 lineage was higher than its closely related triple reassortant swine, which could be attributed to the adaptation to the new host and/or intensive surveillance after the inter-species transmission from swine to human. The positive selection sites were found in all lineages except the equine lineage and mostly in the NS1 region. The positive selection sites 22, 26, 226, 227 and 230 of the human lineage are significant because these residues participate in either forming the dimerization of the two RNA binding domain (RBD) monomers or blocking the replication of host genes. Residues at position 171 provide hydrophobic interactions with hydrophobic residues at p85ß and thus induce viral cell growth. The lineages and evolutionary dynamics of influenza A NS gene obtained in this study, along with the studies of other gene segments, are expected to improve the early detection of new viruses and thus have the potential to enhance influenza surveillance.