Partitioning the genetic diversity of a virus family: approach and evaluation through a case study of picornaviruses

Rapid taxonomic categorization of short, abundant virus sequences for ecological analyses

Public health concerns about recent viral epidemics have motivated researchers to seek novel ways to understand pathogen infection in native, wildlife hosts. With its deep history of tools and perspectives for understanding the abundance and distribution of organisms, ecology can shed new light on viral infection dynamics. However, datasets allowing deep explorations of viral communities from an ecological perspective are lacking. We sampled 1086 bats from two, adjacent Puerto Rican caves and tested them for infection by herpesviruses, resulting in 3131 short, viral sequences. Using percent identity of nucleotides and a machine learning algorithm (affinity propagation), we categorized herpesviruses into 43 operational taxonomic units (OTUs) to be used in place of species in subsequent ecological analyses. Herpesvirus metacommunities demonstrated long-tailed rank frequency distributions at all analyzed levels of host organization (i.e., individual, population, and community). Although 13 herpesvirus OTUs were detected in more than one host species, OTUs generally exhibited host specificity by infecting a single core host species at a significantly higher prevalence than in all satellite species combined. We describe the natural history of herpesvirus metacommunities in Puerto Rican bats and suggest that viruses follow the general law that communities comprise few common and many rare species. To guide future efforts in the field of viral ecology, hypotheses are presented regarding mechanisms that contribute to these patterns.

Deep mining of the Sequence Read Archive reveals major genetic innovations in coronaviruses and other nidoviruses of aquatic vertebrates

Virus discovery by genomics and metagenomics empowered studies of viromes, facilitated characterization of pathogen epidemiology, and redefined our understanding of the natural genetic diversity of viruses with profound functional and structural implications. Here we employed a data-driven virus discovery approach that directly queries unprocessed sequencing data in a highly parallelized way and involves a targeted viral genome assembly strategy in a wide range of sequence similarity. By screening more than 269,000 datasets of numerous authors from the Sequence Read Archive and using two metrics that quantitatively assess assembly quality, we discovered 40 nidoviruses from six virus families whose members infect vertebrate hosts. They form 13 and 32 putative viral subfamilies and genera, respectively, and include 11 coronaviruses with bisegmented genomes from fishes and amphibians, a giant 36.1 kilobase coronavirus genome with a duplicated spike glycoprotein (S) gene, 11 tobaniviruses and 17 additional corona-, arteri-, cremega-, nanhypo- and nangoshaviruses. Genome segmentation emerged in a single evolutionary event in the monophyletic lineage encompassing the subfamily Pitovirinae. We recovered the bisegmented genome sequences of two coronaviruses from RNA samples of 69 infected fishes and validated the presence of poly(A) tails at both segments using 3'RACE PCR and subsequent Sanger sequencing. We report a genetic linkage between accessory and structural proteins whose phylogenetic relationships and evolutionary distances are incongruent with the phylogeny of replicase proteins. We rationalize these observations in a model of inter-family S recombination involving at least five ancestral corona- and tobaniviruses of aquatic hosts. In support of this model, we describe an individual fish co-infected with members from the families Coronaviridae and Tobaniviridae. Our results expand the scale of the known extraordinary evolutionary plasticity in nidoviral genome architecture and call for revisiting fundamentals of genome expression, virus particle biology, host range and ecology of vertebrate nidoviruses.

Translating genomic exploration of the family Polyomaviridae into confident human polyomavirus detection

The Polyomaviridae is a family of ubiquitous dsDNA viruses that establish persistent infection early in life. Screening for human polyomaviruses (HPyVs), which comprise 14 diverse species, relies upon species-specific qPCRs whose validity may be challenged by accelerating genomic exploration of the virosphere. Using this reasoning, we tested 64 published HPyV qPCR assays in silico against the 1781 PyV genome sequences that were divided in targets and nontargets, based on anticipated species specificity of each qPCR. We identified several cases of problematic qPCR performance that were confirmed in vitro and corrected through using degenerate oligos. Furthermore, our study ranked 8 out of 52 tested BKPyV qPCRs as remaining of consistently high quality in the wake of recent PyV discoveries and showed how sensitivity of most other qPCRs could be rescued by annealing temperature adjustment. This study establishes an efficient framework for ensuring confidence in available HPyV qPCRs in the genomic era.

The Neighborhood of the Spike Gene Is a Hotspot for Modular Intertypic Homologous and Nonhomologous Recombination in Coronavirus Genomes

Coronaviruses (CoVs) have very large RNA viral genomes with a distinct genomic architecture of core and accessory open reading frames (ORFs). It is of utmost importance to understand their patterns and limits of homologous and nonhomologous recombination, because such events may affect the emergence of novel CoV strains, alter their host range, infection rate, tissue tropism pathogenicity, and their ability to escape vaccination programs. Intratypic recombination among closely related CoVs of the same subgenus has often been reported; however, the patterns and limits of genomic exchange between more distantly related CoV lineages (intertypic recombination) need further investigation. Here, we report computational/evolutionary analyses that clearly demonstrate a substantial ability for CoVs of different subgenera to recombine. Furthermore, we show that CoVs can obtain-through nonhomologous recombination-accessory ORFs from core ORFs, exchange accessory ORFs with different CoV genera, with other viruses (i.e., toroviruses, influenza C/D, reoviruses, rotaviruses, astroviruses) and even with hosts. Intriguingly, most of these radical events result from double crossovers surrounding the Spike ORF, thus highlighting both the instability and mobile nature of this genomic region. Although many such events have often occurred during the evolution of various CoVs, the genomic architecture of the relatively young SARS-CoV/SARS-CoV-2 lineage so far appears to be stable.

The Remarkable Evolutionary Plasticity of Coronaviruses by Mutation and Recombination: Insights for the COVID-19 Pandemic and the Future Evolutionary Paths of SARS-CoV-2

Coronaviruses (CoVs) constitute a large and diverse subfamily of positive-sense single-stranded RNA viruses. They are found in many mammals and birds and have great importance for the health of humans and farm animals. The current SARS-CoV-2 pandemic, as well as many previous epidemics in humans that were of zoonotic origin, highlights the importance of studying the evolution of the entire CoV subfamily in order to understand how novel strains emerge and which molecular processes affect their adaptation, transmissibility, host/tissue tropism, and patho non-homologous genicity. In this review, we focus on studies over the last two years that reveal the impact of point mutations, insertions/deletions, and intratypic/intertypic homologous and non-homologous recombination events on the evolution of CoVs. We discuss whether the next generations of CoV vaccines should be directed against other CoV proteins in addition to or instead of spike. Based on the observed patterns of molecular evolution for the entire subfamily, we discuss five scenarios for the future evolutionary path of SARS-CoV-2 and the COVID-19 pandemic. Finally, within this evolutionary context, we discuss the recently emerged Omicron (B.1.1.529) VoC.

Bioinformatics of virus taxonomy: foundations and tools for developing sequence-based hierarchical classification

The genome sequence is the only characteristic readily obtainable for all known viruses, underlying the growing role of comparative genomics in organizing knowledge about viruses in a systematic evolution-aware way, known as virus taxonomy. Overseen by the International Committee on Taxonomy of Viruses (ICTV), development of virus taxonomy involves taxa demarcation at 15 ranks of a hierarchical classification, often in host-specific manner. Outside the ICTV remit, researchers assess fitting numerous unclassified viruses into the established taxa. They employ different metrics of virus clustering, basing on conserved domain(s), separation of viruses in rooted phylogenetic trees and pair-wise distance space. Computational approaches differ further in respect to methodology, number of ranks considered, sensitivity to uneven virus sampling, and visualization of results. Advancing and using computational tools will be critical for improving taxa demarcation across the virosphere and resolving rank origins in research that may also inform experimental virology.

A k-mer based approach for classifying viruses without taxonomy identifies viral associations in human autism and plant microbiomes

Viruses are an underrepresented taxa in the study and identification of microbiome constituents; however, they play an essential role in health, microbiome regulation, and transfer of genetic material. Only a few thousand viruses have been isolated, sequenced, and assigned a taxonomy, which limits the ability to identify and quantify viruses in the microbiome. Additionally, the vast diversity of viruses represents a challenge for classification, not only in constructing a viral taxonomy, but also in identifying similarities between a virus' genotype and its phenotype. However, the diversity of viral sequences can be leveraged to classify their sequences in metagenomic and metatranscriptomic samples, even if they do not have a taxonomy. To identify and quantify viruses in transcriptomic and genomic samples, we developed a dynamic programming algorithm for creating a classification tree out of 715,672 metagenome viruses. To create the classification tree, we clustered proportional similarity scores generated from the k-mer profiles of each of the metagenome viruses to create a database of metagenomic viruses. The resulting Kraken2 database of the metagenomic viruses can be found here: https://www.osti.gov/biblio/1615774 and is compatible with Kraken2. We then integrated the viral classification database with databases created with genomes from NCBI for use with ParaKraken (a parallelized version of Kraken provided in Supplemental Zip 1), a metagenomic/transcriptomic classifier. To illustrate the breadth of our utility for classifying metagenome viruses, we analyzed data from a plant metagenome study identifying genotypic and compartment specific differences between two Populus genotypes in three different compartments. We also identified a significant increase in abundance of eight viral sequences in post mortem brains in a human metatranscriptome study comparing Autism Spectrum Disorder patients and controls. We also show the potential accuracy for classifying viruses by utilizing both the JGI and NCBI viral databases to identify the uniqueness of viral sequences. Finally, we validate the accuracy of viral classification with NCBI databases containing viruses with taxonomy to identify pathogenic viruses in known COVID-19 and cassava brown streak virus infection samples. Our method represents the compulsory first step in better understanding the role of viruses in the microbiome by allowing for a more complete identification of sequences without taxonomy. Better classification of viruses will improve identifying associations between viruses and their hosts as well as viruses and other microbiome members. Despite the lack of taxonomy, this database of metagenomic viruses can be used with any tool that utilizes a taxonomy, such as Kraken, for accurate classification of viruses.

Novel Filoviruses, Hantavirus, and Rhabdovirus in Freshwater Fish, Switzerland, 2017

European perch (Perca fluviatilis) are increasingly farmed as a human food source. Viral infections of European perch remain largely unexplored, thereby putting farm populations at incalculable risk for devastating fish epizootics and presenting a potential hazard to consumers. To address these concerns, we applied metatranscriptomics to identify disease-associated viruses in European perch farmed in Switzerland. Unexpectedly, in clinically diseased fish we detected novel freshwater fish filoviruses, a novel freshwater fish hantavirus, and a previously unknown rhabdovirus. Hantavirus titers were high, and we demonstrated virus in macrophages and gill endothelial cells by using in situ hybridization. Rhabdovirus titers in organ samples were low, but virus could be isolated on cell culture. Our data add to the hypothesis that filoviruses, hantaviruses, and rhabdoviruses are globally distributed common fish commensals, pathogens, or both. Our findings shed new light on negative-sense RNA virus diversity and evolution.

Perspective on taxonomic classification of uncultivated viruses

Historically, virus taxonomy has been limited to describing viruses that were readily cultivated in the laboratory or emerging in natural biomes. Metagenomic analyses, single-particle sequencing, and database mining efforts have yielded new sequence data on an astounding number of previously unknown viruses. As metagenomes are relatively free of biases, these data provide an unprecedented insight into the vastness of the virosphere, but to properly value the extent of this diversity it is critical that the viruses are taxonomically classified. Inclusion of uncultivated viruses has already improved the process as well as the understanding of the taxa, viruses, and their evolutionary relationships. The continuous development and testing of computational tools will be required to maintain a dynamic virus taxonomy that can accommodate the new discoveries.

Taxonomy classification using genomic footprint of mitochondrial sequences

BACKGROUND

Advancement in the sequencing technology yields a huge number of genomes of a multitude of organisms in our planet. One of the fundamental tasks for processing and analyzing these sequences is to organize them in the existing taxonomic orders. <P> Method: Recently we proposed a novel approach, GenFooT, of taxonomy classification using the concept of genomic footprint (GFP). The technique is further refined and enhanced in this work leading to improved accuracies in the task of taxonomic classification on various benchmark datasets. GenFooT maps a genome sequence in a 2D coordinate space and extracts features from that representation. It uses two hyper-parameters, namely block size and number of fragments of genomic sequence while computing the feature. In this work, we propose an analysis for choosing values of those parameters adaptively from the sequences. The enhanced version of GenFooT is named GenFooT2. <P> Results and Conclusion: We have experimented GenFooT2 on ten different biological datasets of genomic sequences of various organisms belonging to different taxonomy ranks. Our experimental results indicate more than 3% improved classification performance of the proposed features with Logistic regression classifier than the GenFooT. We also performed the statistical test to compare the performance of GenFooT2 with the state-of-the-art methods including our previous method GenFooT. The experimental results as well as the statistical test exhibit that the performance of the proposed GenFooT2 is significantly better.

SARS-CoV-2 gene content and COVID-19 mutation impact by comparing 44 Sarbecovirus genomes

Despite its clinical importance, the SARS-CoV-2 gene set remains unresolved, hindering dissection of COVID-19 biology. We use comparative genomics to provide a high-confidence protein-coding gene set, characterize evolutionary constraint, and prioritize functional mutations. We select 44 Sarbecovirus genomes at ideally-suited evolutionary distances, and quantify protein-coding evolutionary signatures and overlapping constraint. We find strong protein-coding signatures for ORFs 3a, 6, 7a, 7b, 8, 9b, and a novel alternate-frame gene, ORF3c, whereas ORFs 2b, 3d/3d-2, 3b, 9c, and 10 lack protein-coding signatures or convincing experimental evidence of protein-coding function. Furthermore, we show no other conserved protein-coding genes remain to be discovered. Mutation analysis suggests ORF8 contributes to within-individual fitness but not person-to-person transmission. Cross-strain and within-strain evolutionary pressures agree, except for fewer-than-expected within-strain mutations in nsp3 and S1, and more-than-expected in nucleocapsid, which shows a cluster of mutations in a predicted B-cell epitope, suggesting immune-avoidance selection. Evolutionary histories of residues disrupted by spike-protein substitutions D614G, N501Y, E484K, and K417N/T provide clues about their biology, and we catalog likely-functional co-inherited mutations. Previously reported RNA-modification sites show no enrichment for conservation. Here we report a high-confidence gene set and evolutionary-history annotations providing valuable resources and insights on SARS-CoV-2 biology, mutations, and evolution.

Recognizing species as a new focus of virus research

Species taxa are the units of taxonomy most suited to measure virus diversity, and they account for more than 70% of all virus taxa. Yet, as evidenced by the content of GenBank entries and illustrated by the recent literature on SARS-CoV-2, they are the most neglected taxa of virus research. To correct this disparity, we propose to make species taxa a first choice for communicating virus taxonomy in publications concerning viruses. We see it as a key step toward promoting research on diverse viruses, including pathogens, at this fundamental level of biology.

Coronavirus replication-transcription complex: Vital and selective NMPylation of a conserved site in nsp9 by the NiRAN-RdRp subunit

RNA-dependent RNA polymerases (RdRps) of the Nidovirales (Coronaviridae, Arteriviridae, and 12 other families) are linked to an amino-terminal (N-terminal) domain, called NiRAN, in a nonstructural protein (nsp) that is released from polyprotein 1ab by the viral main protease (M^pro). Previously, self-GMPylation/UMPylation activities were reported for an arterivirus NiRAN-RdRp nsp and suggested to generate a transient state primed for transferring nucleoside monophosphate (NMP) to (currently unknown) viral and/or cellular biopolymers. Here, we show that the coronavirus (human coronavirus [HCoV]-229E and severe acute respiratory syndrome coronavirus 2) nsp12 (NiRAN-RdRp) has Mn²⁺-dependent NMPylation activity that catalyzes the transfer of a single NMP to the cognate nsp9 by forming a phosphoramidate bond with the primary amine at the nsp9 N terminus (N3825) following M^pro-mediated proteolytic release of nsp9 from N-terminally flanking nsps. Uridine triphosphate was the preferred nucleotide in this reaction, but also adenosine triphosphate, guanosine triphosphate, and cytidine triphosphate were suitable cosubstrates. Mutational studies using recombinant coronavirus nsp9 and nsp12 proteins and genetically engineered HCoV-229E mutants identified residues essential for NiRAN-mediated nsp9 NMPylation and virus replication in cell culture. The data corroborate predictions on NiRAN active-site residues and establish an essential role for the nsp9 N3826 residue in both nsp9 NMPylation in vitro and virus replication. This residue is part of a conserved N-terminal NNE tripeptide sequence and shown to be the only invariant residue in nsp9 and its homologs in viruses of the family Coronaviridae The study provides a solid basis for functional studies of other nidovirus NMPylation activities and suggests a possible target for antiviral drug development.

Reassessing species demarcation criteria in viroid taxonomy by pairwise identity matrices

With a small, circular and non-protein coding RNA genome, viroids are the smallest infectious agents. They invade plants, which in turn may develop symptoms. Since their discovery about 50 years ago, more than thirty viroids have been reported and classified using as species demarcation less than 90 per cent sequence identity on the overall genome and evidence of biological divergence with respect to the closest related viroids. In the last few years, new viroids have been identified that infect latently their (frequently) woody hosts and have a narrow experimental hosts range, complicating and slowing down studies on their biology. As a consequence, several viroids are still waiting for classification. Moreover, the number of new viroids is expected to increase in the next years due to the use of high-throughput sequencing technologies with diagnostics purposes. Therefore, establishment of reliable species demarcation criteria mainly based on molecular features of viroids is needed. Here, viroid classification is reassessed and a scheme based on pairwise sequence identity matrices is developed. After identifying a threshold pairwise identity score (PWIS) for each viroid genus, to be used as a species demarcation criterion, we show that most of those yet unclassified viroids can be assigned to a known or to a new species, thus limiting the need for additional biological evidence to only a few more complex situations. The advantages of this PWIS-based method are that the proposed identity thresholds for species demarcations are not arbitrarily established and evidence for biological divergence is not mandatory. Importantly, the current classification is not essentially modified. A protocol for a tentative fast classification of new viroids according to the proposed approach is also provided.

Using affinity propagation clustering for identifying bacterial clades and subclades with whole-genome sequences of Francisella tularensis

By combining a reference-independent SNP analysis and average nucleotide identity (ANI) with affinity propagation clustering (APC), we developed a significantly improved methodology allowing resolving phylogenetic relationships, based on objective criteria. These bioinformatics tools can be used as a general ruler to determine phylogenetic relationships and clustering of bacteria, exemplary done with Francisella (F.) tularensis. Molecular epidemiology of F. tularensis is currently assessed mostly based on laboratory methods and molecular analysis. The high evolutionary stability and the clonal nature makes Francisella ideal for subtyping with single nucleotide polymorphisms (SNPs). Sequencing and real-time PCR can be used to validate the SNP analysis. We investigate whole-genome sequences of 155 F. tularensis subsp. holarctica isolates. Phylogenetic testing was based on SNPs and average nucleotide identity (ANI) as reference independent, alignment-free methods taking small-scale and large-scale differences within the genomes into account. Especially the whole genome SNP analysis with kSNP3.0 allowed deciphering quite subtle signals of systematic differences in molecular variation. Affinity propagation clustering (APC) resulted in three clusters showing the known clades B.4, B.6, and B.12. These data correlated with the results of real-time PCR assays targeting canSNPs loci. Additionally, we detected two subtle sub-clusters. SplitsTree was used with standard-setting using the aligned SNPs from Parsnps. Together APC, HierBAPS, and SplitsTree enabled us to generate hypotheses about epidemiologic relationships between bacterial clusters and describing the distribution of isolates. Our data indicate that the choice of the typing technique can increase our understanding of the pathogenesis and transmission of diseases with the eventual for prevention. This is opening perspectives to be applied to other bacterial species. The data provide evidence that Germany might be the collision zone where the clade B.12, also known as the East European clade, overlaps with the clade B.6, also known as the Iberian clade. Described methods allow generating a new, more detailed perspective for F. tularensis subsp. holarctica phylogeny. These results may encourage to determine phylogenetic relationships and clustering of other bacteria the same way.

Phlebovirus sequences detected in ticks collected in Russia: Novel phleboviruses, distinguishing criteria and high tick specificity

Phlebovirus is an abundant and rather heterogeneous genus within the Phenuiviridae family (order Bunyavirales). The genus Phlebovirus is divided into two antigenic complexes, which also correspond to the main vector: sandflies/mosquitoes and ticks. Previously, only sandfly/mosquito-borne phleboviruses were associated with human disease, such as Rift Valley fever virus, Toscana virus, Sicilian and Naples Sandfly fever viruses and others. Until recently, tick-borne phleboviruses were not considered as human pathogens. After the discovery of severe fever with thrombocytopenia syndrome, interest to tick-borne phleboviruses has increased dramatically. In the last decade, many novel phleboviruses have been reported in different regions. Despite this, the diversity, ecology and pathogenicity of these viruses still remain obscure. The aim of this work was to study the diversity of phleboviruses in ticks collected in several regions of Russia. We used pan-phlebovirus RT-PCR assays based on multiple degenerate primers targeting the polymerase gene fragment. Arthropod specimens were collected from 2005 to 2018. A total of 5901 Ixodidae ticks combined into 1116 pools were screened. A total of 160 specific amplicons were produced. In three cases RT-PCR assays amplified two distinct viruses from same tick pools. Direct sequencing of amplicons and subsequent phylogenetic analysis revealed twelve representatives of divergent phlebovirus groups. Based on the distribution of pairwise nucleotide sequence identity values, a cut-off (88%) was suggested to distinguish tick-borne phleboviruses. According to this provisional criterion, two viruses found here could be termed novel, while ten viruses have been described in previous studies. Detected phleboviruses demonstrated almost perfect specificity to a tick species or, at least, a genus. The same pattern was observed for tick-borne phleboviruses found in different studies around the world. Viruses that grouped together on a phylogenetic tree and differed less than this sequence identity threshold suggested above were hosted by ticks from the same genus.

Four point-of-care lateral flow immunoassays for diagnosis of COVID-19 and for assessing dynamics of antibody responses to SARS-CoV-2

OBJECTIVES

We aimed to evaluate the role of rapid serological tests in the management of coronavirus disease 2019 (COVID-19) patients.

METHODS

This retrospective study enrolled 16 real-time reverse transcription polymerase chain reaction-confirmed symptomatic patients with COVID-19 and 58 COVID-19 negative patients at a medical center in Taiwan over a 3-month period. Serial serum samples were collected and tested for antibody response using four point-of-care (POC) lateral flow immunoassays (LFIA) (ALLTEST 2019-nCoV IgG/IgM Rapid Test, Dynamiker 2019-nCoV IgG/IgM Rapid Test, ASK COVID-19 IgG/IgM Rapid Test, and Wondfo SARS-CoV-2 Antibody Test). Time-dependent detection sensitivity and timeliness of seroconversion were determined and compared between the four POC rapid tests.

RESULTS

The overall sensitivity and specificity of the four tests for detecting anti-SARS-CoV-2 antibodies after 3 weeks of symptom onset were 100% and 100%, respectively. There was no significant difference between the rapid tests used for detection of IgM and IgG separately and those used for detection of combined total antibody (mainly IgM/IgG). There was no significant difference between the four POC rapid tests in terms of time required for determining seroconversion of COVID-19. Patients with COVID-19 with pneumonia demonstrated shorter seroconversion time than those without pneumonia.

CONCLUSION

Though the POC antibody rapid tests based on LFIA showed reliable performance in the detection of SARS-CoV-2-specific antibodies, the results of these tests should be interpreted and applied appropriately in the context of antibody dynamic of COVID-19 infection. COVID-19 patients complicated with pneumonia exhibited earlier anti-SARS-CoV-2 antibody response than COVID-19 patients without pneumonia.

The new scope of virus taxonomy: partitioning the virosphere into 15 hierarchical ranks

Virus taxonomy emerged as a discipline in the middle of the twentieth century. Traditionally, classification by virus taxonomists has been focussed on the grouping of relatively closely related viruses. However, during the past few years, the International Committee on Taxonomy of Viruses (ICTV) has recognized that the taxonomy it develops can be usefully extended to include the basal evolutionary relationships among distantly related viruses. Consequently, the ICTV has changed its Code to allow a 15-rank classification hierarchy that closely aligns with the Linnaean taxonomic system and may accommodate the entire spectrum of genetic divergence in the virosphere. The current taxonomies of three human pathogens, Ebola virus, severe acute respiratory syndrome coronavirus and herpes simplex virus 1 are used to illustrate the impact of the expanded rank structure. This new rank hierarchy of virus taxonomy will stimulate further research on virus origins and evolution, and vice versa, and could promote crosstalk with the taxonomies of cellular organisms.

The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

The present outbreak of a coronavirus-associated acute respiratory disease called coronavirus disease 19 (COVID-19) is the third documented spillover of an animal coronavirus to humans in only two decades that has resulted in a major epidemic. The Coronaviridae Study Group (CSG) of the International Committee on Taxonomy of Viruses, which is responsible for developing the classification of viruses and taxon nomenclature of the family Coronaviridae, has assessed the placement of the human pathogen, tentatively named 2019-nCoV, within the Coronaviridae. Based on phylogeny, taxonomy and established practice, the CSG recognizes this virus as forming a sister clade to the prototype human and bat severe acute respiratory syndrome coronaviruses (SARS-CoVs) of the species Severe acute respiratory syndrome-related coronavirus, and designates it as SARS-CoV-2. In order to facilitate communication, the CSG proposes to use the following naming convention for individual isolates: SARS-CoV-2/host/location/isolate/date. While the full spectrum of clinical manifestations associated with SARS-CoV-2 infections in humans remains to be determined, the independent zoonotic transmission of SARS-CoV and SARS-CoV-2 highlights the need for studying viruses at the species level to complement research focused on individual pathogenic viruses of immediate significance. This will improve our understanding of virus–host interactions in an ever-changing environment and enhance our preparedness for future outbreaks.

The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

The present outbreak of a coronavirus-associated acute respiratory disease called coronavirus disease 19 (COVID-19) is the third documented spillover of an animal coronavirus to humans in only two decades that has resulted in a major epidemic. The Coronaviridae Study Group (CSG) of the International Committee on Taxonomy of Viruses, which is responsible for developing the classification of viruses and taxon nomenclature of the family Coronaviridae, has assessed the placement of the human pathogen, tentatively named 2019-nCoV, within the Coronaviridae. Based on phylogeny, taxonomy and established practice, the CSG recognizes this virus as forming a sister clade to the prototype human and bat severe acute respiratory syndrome coronaviruses (SARS-CoVs) of the species Severe acute respiratory syndrome-related coronavirus, and designates it as SARS-CoV-2. In order to facilitate communication, the CSG proposes to use the following naming convention for individual isolates: SARS-CoV-2/host/location/isolate/date. While the full spectrum of clinical manifestations associated with SARS-CoV-2 infections in humans remains to be determined, the independent zoonotic transmission of SARS-CoV and SARS-CoV-2 highlights the need for studying viruses at the species level to complement research focused on individual pathogenic viruses of immediate significance. This will improve our understanding of virus–host interactions in an ever-changing environment and enhance our preparedness for future outbreaks.

The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

The present outbreak of a coronavirus-associated acute respiratory disease called coronavirus disease 19 (COVID-19) is the third documented spillover of an animal coronavirus to humans in only two decades that has resulted in a major epidemic. The Coronaviridae Study Group (CSG) of the International Committee on Taxonomy of Viruses, which is responsible for developing the classification of viruses and taxon nomenclature of the family Coronaviridae, has assessed the placement of the human pathogen, tentatively named 2019-nCoV, within the Coronaviridae. Based on phylogeny, taxonomy and established practice, the CSG recognizes this virus as forming a sister clade to the prototype human and bat severe acute respiratory syndrome coronaviruses (SARS-CoVs) of the species Severe acute respiratory syndrome-related coronavirus, and designates it as SARS-CoV-2. In order to facilitate communication, the CSG proposes to use the following naming convention for individual isolates: SARS-CoV-2/host/location/isolate/date. While the full spectrum of clinical manifestations associated with SARS-CoV-2 infections in humans remains to be determined, the independent zoonotic transmission of SARS-CoV and SARS-CoV-2 highlights the need for studying viruses at the species level to complement research focused on individual pathogenic viruses of immediate significance. This will improve our understanding of virus–host interactions in an ever-changing environment and enhance our preparedness for future outbreaks.

The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

The present outbreak of a coronavirus-associated acute respiratory disease called coronavirus disease 19 (COVID-19) is the third documented spillover of an animal coronavirus to humans in only two decades that has resulted in a major epidemic. The Coronaviridae Study Group (CSG) of the International Committee on Taxonomy of Viruses, which is responsible for developing the classification of viruses and taxon nomenclature of the family Coronaviridae, has assessed the placement of the human pathogen, tentatively named 2019-nCoV, within the Coronaviridae. Based on phylogeny, taxonomy and established practice, the CSG recognizes this virus as forming a sister clade to the prototype human and bat severe acute respiratory syndrome coronaviruses (SARS-CoVs) of the species Severe acute respiratory syndrome-related coronavirus, and designates it as SARS-CoV-2. In order to facilitate communication, the CSG proposes to use the following naming convention for individual isolates: SARS-CoV-2/host/location/isolate/date. While the full spectrum of clinical manifestations associated with SARS-CoV-2 infections in humans remains to be determined, the independent zoonotic transmission of SARS-CoV and SARS-CoV-2 highlights the need for studying viruses at the species level to complement research focused on individual pathogenic viruses of immediate significance. This will improve our understanding of virus–host interactions in an ever-changing environment and enhance our preparedness for future outbreaks.

First Report of a Mesonivirus and Its Derived Small RNAs in an Aphid Species Aphis citricidus (Hemiptera: Aphididae), Implying Viral Infection Activity

We report a new positive-sense single-stranded RNA (ss RNA+) virus from the brown citrus aphid Aphis citricidus. The 20,300 nucleotide (nt)-long viral genome contains five open-reading frames and encodes six conserved domains (TM2, 3CLpro, TM3, RdRp, Zm, and HEL1). Phylogenetic analysis and amino acid sequence analysis revealed this virus might belong to an unassigned genus in the family Mesoniviridae. The presence of the virus was also confirmed in the field population. Importantly, analysis of the virus-derived small RNAs showed a 22-nt peak, implying that viral infection triggers the small interfering RNA pathway as antiviral immunity in aphids. This is the first report of a mesonivirus in invertebrates other than mosquitoes.

Re-assessing the diversity of negative strand RNA viruses in insects

The spectrum of viruses in insects is important for subjects as diverse as public health, veterinary medicine, food production, and biodiversity conservation. The traditional interest in vector-borne diseases of humans and livestock has drawn the attention of virus studies to hematophagous insect species. However, these represent only a tiny fraction of the broad diversity of Hexapoda, the most speciose group of animals. Here, we systematically probed the diversity of negative strand RNA viruses in the largest and most representative collection of insect transcriptomes from samples representing all 34 extant orders of Hexapoda and 3 orders of Entognatha, as well as outgroups, altogether representing 1243 species. Based on profile hidden Markov models we detected 488 viral RNA-directed RNA polymerase (RdRp) sequences with similarity to negative strand RNA viruses. These were identified in members of 324 arthropod species. Selection for length, quality, and uniqueness left 234 sequences for analyses, showing similarity to genomes of viruses classified in Bunyavirales (n = 86), Articulavirales (n = 54), and several orders within Haploviricotina (n = 94). Coding-complete genomes or nearly-complete subgenomic assemblies were obtained in 61 cases. Based on phylogenetic topology and the availability of coding-complete genomes we estimate that at least 20 novel viral genera in seven families need to be defined, only two of them monospecific. Seven additional viral clades emerge when adding sequences from the present study to formerly monospecific lineages, potentially requiring up to seven additional genera. One long sequence may indicate a novel family. For segmented viruses, cophylogenies between genome segments were generally improved by the inclusion of viruses from the present study, suggesting that in silico misassembly of segmented genomes is rare or absent. Contrary to previous assessments, significant virus-host codivergence was identified in major phylogenetic lineages based on two different approaches of codivergence analysis in a hypotheses testing framework. In spite of these additions to the known spectrum of viruses in insects, we caution that basing taxonomic decisions on genome information alone is challenging due to technical uncertainties, such as the inability to prove integrity of complete genome assemblies of segmented viruses.

A novel alignment-free method for HIV-1 subtype classification

HIV-1 is the most common and pathogenic strain of human immunodeficiency virus consisting of many subtypes. To study the difference among HIV-1 subtypes in infection, diagnosis and drug design, it is important to identify HIV-1 subtypes from clinical HIV-1 samples. In this work, we propose an effective numeric representation called Subsequence Natural Vector (SNV) to encode HIV-1 sequences. Using the representation, we introduce an improved linear discriminant analysis method to classify HIV-1 viruses correctly. SNV is based on distribution of nucleotides in HIV-1 viral sequences. It not only computes the number of nucleotides, but also describes the position and variance of nucleotides in viruses. To validate our alignment-free method, 6902 complete genomes and 11,668 pol gene sequences of HIV-1 subtypes were collected from the up-to-date Los Alamos HIV database. SNV outperforms the three popular methods, Kameris, Comet and REGA, with almost 100% Sensitivity and Specificity, also with much less time. Our subtyping algorithm especially works better for circulating recombinant forms (CRFs) consisting of a few sequences. Our approach is also powerful to separate unique recombinant forms (URFs) from other subtypes with 100% Sensitivity and Specificity. Moreover, phylogenetic trees based on SNV representation are constructed using full-length HIV-1 genomes and pol genes respectively, where viruses from the same subtype are clustered together correctly.

Hantaviridae: Current Classification and Future Perspectives

In recent years, negative-sense RNA virus classification and taxon nomenclature have undergone considerable transformation. In 2016, the new order Bunyavirales was established, elevating the previous genus Hantavirus to family rank, thereby creating Hantaviridae. Here we summarize affirmed taxonomic modifications of this family from 2016 to 2019. Changes involve the admission of >30 new hantavirid species and the establishment of subfamilies and novel genera based on DivErsity pArtitioning by hieRarchical Clustering (DEmARC) analysis of genomic sequencing data. We outline an objective framework that can be used in future classification schemes when more hantavirids sequences will be available. Finally, we summarize current taxonomic proposals and problems in hantavirid taxonomy that will have to be addressed shortly.

Simulation data for the estimation of numerical constants for approximating pairwise evolutionary distances between amino acid sequences

Estimating the number of substitution events per site that have occurred during the evolution of a pair of amino acid sequences is a common task in phylogenetics and comparative genomics that often requires quite slow maximum-likelihood procedures when taking into account explicit evolutionary models. Data presented in this article are large sets of numbers of substitution events and associated numbers of observed differences between pairs of aligned amino acid sequences that have been generated through a simulation procedure of sequence evolution under a broad range of evolutionary models. These data are available at https://zenodo.org/record/2653704 (doi:10.5281/zenodo.2653704). They are accompanied in this paper by figures showing the strong relationship between the corresponding evolutionary and uncorrected distances, as well as estimated numerical constants that determine non-linear functions that fit the simulated data. These numerical constants can be useful to quickly estimate pairwise evolutionary distances directly from uncorrected distances between aligned amino acid sequences.

Exploiting the Legacy of the Arbovirus Hunters

In recent years, it has become evident that a generational gap has developed in the community of arbovirus research. This apparent gap is due to the dis-investment of training for the next generation of arbovirologists, which threatens to derail the rich history of virus discovery, field epidemiology, and understanding of the richness of diversity that surrounds us. On the other hand, new technologies have resulted in an explosion of virus discovery that is constantly redefining the virosphere and the evolutionary relationships between viruses. This paradox presents new challenges that may have immediate and disastrous consequences for public health when yet to be discovered arboviruses emerge. In this review we endeavor to bridge this gap by providing a historical context for the work being conducted today and provide continuity between the generations. To this end, we will provide a narrative of the thrill of scientific discovery and excitement and the challenges lying ahead.

Modeling of Mutational Events in the Evolution of Viruses

Diverse studies of viral evolution have led to the recognition that the evolutionary rates of viral taxa observed are dependent on the time scale being investigated—with short-term studies giving fast substitution rates, and orders of magnitude lower rates for deep calibrations. Although each of these factors may contribute to this time dependent rate phenomenon, a more fundamental cause should be considered. We sought to test computationally whether the basic phenomena of virus evolution (mutation, replication, and selection) can explain the relationships between the evolutionary and phylogenetic distances. We tested, by computational inference, the hypothesis that the phylogenetic distances between the pairs of sequences are functions of the evolutionary path lengths between them. A Basic simulation revealed that the relationship between simulated genetic and mutational distances is non-linear, and can be consistent with different rates of nucleotide substitution at different depths of branches in phylogenetic trees.

Structural basis for catalysis and substrate specificity of a 3C-like cysteine protease from a mosquito mesonivirus

Cavally virus (CavV) is a mosquito-borne plus-strand RNA virus in the family Mesoniviridae (order Nidovirales). We present X-ray structures for the CavV 3C-like protease (3CL^pro), as a free enzyme and in complex with a peptide aldehyde inhibitor mimicking the P4-to-P1 residues of a natural substrate. The 3CL^pro structure (refined to 1.94 Å) shows that the protein forms dimers. The monomers are comprised of N-terminal domains I and II, which adopt a chymotrypsin-like fold, and a C-terminal α-helical domain III. The catalytic Cys-His dyad is assisted by a complex network of interactions involving a water molecule that mediates polar contacts between the catalytic His and a conserved Asp located in the domain II-III junction and is suitably positioned to stabilize the developing positive charge of the catalytic His in the transition state during catalysis. The study also reveals the structural basis for the distinct P2 Asn-specific substrate-binding pocket of mesonivirus 3CL^pros.

•

First structure of a 3CL^pro of an invertebrate RNA virus.

•

Structural basis of the unique substrate specificity defined by Asn at the P2 position of mesonivirus 3CL^pro substrates.

•

Emerging role of a conserved Asp residue that assists the Cys-His catalytic dyad in vertebrate and invertebrate 3CL^pros.

Impact of rotavirus and hepatitis A virus by worldwide climatic changes during the period between 2000 and 2013

Enteric viruses are present in the environment as a result of the discharge of poorly or untreated wastewater. The spread of enteric viruses in the environment depend to human activities like stools of infected individuals ejected in the external environment can be transmitted by water sources and back to susceptible individuals for other cycles of illness. Among the enteric viruses Rotaviruses (RV) and Hepatitis A viruses (HAV) is the most detected in wastewater causing gastroenteritis and acute hepatitis. Therefore, it is of interest to climate change, mainly temperature and carbon Dioxide (CO2) variations, on Rotavirus and Hepatitis A as a model of enteric viruses present in the aquatic environment using computational modelling tools. The results of genetic ratio showed a negative correlation between the epidemiological data and the mutation rate. However, the correlation was positive between the temperature, CO2 increase, and the rate of mutation. The positive correlation is explained by the adaptation of the viruses to the climatic changes, the RNA polymerase of the RV induces errors to adapt to the environmental conditions. The simultaneous increase in number of infections and temperature in 2010 has been demonstrated in previous studies deducing that viral pathogenicity increase with temperature increase.

A planarian nidovirus expands the limits of RNA genome size

RNA viruses are the only known RNA-protein (RNP) entities capable of autonomous replication (albeit within a permissive host environment). A 33.5 kilobase (kb) nidovirus has been considered close to the upper size limit for such entities; conversely, the minimal cellular DNA genome is in the 100–300 kb range. This large difference presents a daunting gap for the transition from primordial RNP to contemporary DNA-RNP-based life. Whether or not RNA viruses represent transitional steps towards DNA-based life, studies of larger RNA viruses advance our understanding of the size constraints on RNP entities and the role of genome size in virus adaptation. For example, emergence of the largest previously known RNA genomes (20–34 kb in positive-stranded nidoviruses, including coronaviruses) is associated with the acquisition of a proofreading exoribonuclease (ExoN) encoded in the open reading frame 1b (ORF1b) in a monophyletic subset of nidoviruses. However, apparent constraints on the size of ORF1b, which encodes this and other key replicative enzymes, have been hypothesized to limit further expansion of these viral RNA genomes. Here, we characterize a novel nidovirus (planarian secretory cell nidovirus; PSCNV) whose disproportionately large ORF1b-like region including unannotated domains, and overall 41.1-kb genome, substantially extend the presumed limits on RNA genome size. This genome encodes a predicted 13,556-aa polyprotein in an unconventional single ORF, yet retains canonical nidoviral genome organization and expression, as well as key replicative domains. These domains may include functionally relevant substitutions rarely or never before observed in highly conserved sites of RdRp, NiRAN, ExoN and 3CLpro. Our evolutionary analysis suggests that PSCNV diverged early from multi-ORF nidoviruses, and acquired additional genes, including those typical of large DNA viruses or hosts, e.g. Ankyrin and Fibronectin type II, which might modulate virus-host interactions. PSCNV's greatly expanded genome, proteomic complexity, and unique features–impressive in themselves–attest to the likelihood of still-larger RNA genomes awaiting discovery.

RNA viruses are the only known RNA-protein (RNP) entities capable of autonomous replication. The upper genome size for such entities was assumed to be <35 kb; conversely, the minimal cellular DNA genome is in the 100–300 kilobase (kb) range. This large difference presents a daunting gap for the proposed evolution of contemporary DNA-RNP-based life from primordial RNP entities. Here, we describe a nidovirus from planarians, named planarian secretory cell nidovirus (PSCNV), whose 41.1 kb genome is 23% larger than any riboviral genome yet discovered. This increase is nearly equivalent in size to the entire poliovirus genome, and it equips PSCNV with an unprecedented extra coding capacity to adapt. PSCNV has broken apparent constraints on the size of the genomic subregion that encodes core replication machinery in other nidoviruses, including coronaviruses, and has acquired genes not previously observed in RNA viruses. This virus challenges and advances our understanding of the limits to RNA genome size.

A Novel Hantavirus of the European Mole, Bruges Virus, Is Involved in Frequent Nova Virus Coinfections

Hantaviruses are zoonotic viruses with a complex evolutionary history of virus–host coevolution and cross-species transmission. Although hantaviruses have a broad reservoir host range, virus–host relationships were previously thought to be strict, with a single virus species infecting a single host species. Here, we describe Bruges virus, a novel hantavirus harbored by the European mole (Talpa europaea), which is the well-known host of Nova virus. Phylogenetic analyses of all three genomic segments showed tree topology inconsistencies, suggesting that Bruges virus has emerged from cross-species transmission and ancient reassortment events. A high number of coinfections with Bruges and Nova viruses was detected, but no evidence was found for reassortment between these two hantaviruses. These findings highlight the complexity of hantavirus evolution and the importance of further investigation of hantavirus–reservoir relationships.

Increasing the number of available ranks in virus taxonomy from five to ten and adopting the Baltimore classes as taxa at the basal rank

This opinion article makes a case for increasing the number of ranks used in virus taxonomy from the current five to ten (as are used to classify cellular life forms) and placing the Baltimore classes in the proposed basal rank of domain. These suggestions aim at initiating the process of accommodation of Baltimore classes in virus taxonomy and extension of the virus taxonomy scale to encompass also the most distant relationships.

ViCTree: an automated framework for taxonomic classification from protein sequences

Motivation

The increasing rate of submission of genetic sequences into public databases is providing a growing resource for classifying the organisms that these sequences represent. To aid viral classification, we have developed ViCTree, which automatically integrates the relevant sets of sequences in NCBI GenBank and transforms them into an interactive maximum likelihood phylogenetic tree that can be updated automatically. ViCTree incorporates ViCTreeView, which is a JavaScript-based visualization tool that enables the tree to be explored interactively in the context of pairwise distance data.

Results

To demonstrate utility, ViCTree was applied to subfamily Densovirinae of family Parvoviridae. This led to the identification of six new species of insect virus.

Availability and implementation

ViCTree is open-source and can be run on any Linux- or Unix-based computer or cluster. A tutorial, the documentation and the source code are available under a GPL3 license, and can be accessed at http://bioinformatics.cvr.ac.uk/victree_web/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Defining objective clusters for rabies virus sequences using affinity propagation clustering

Rabies is caused by lyssaviruses, and is one of the oldest known zoonoses. In recent years, more than 21,000 nucleotide sequences of rabies viruses (RABV), from the prototype species rabies lyssavirus, have been deposited in public databases. Subsequent phylogenetic analyses in combination with metadata suggest geographic distributions of RABV. However, these analyses somewhat experience technical difficulties in defining verifiable criteria for cluster allocations in phylogenetic trees inviting for a more rational approach. Therefore, we applied a relatively new mathematical clustering algorythm named ‘affinity propagation clustering’ (AP) to propose a standardized sub-species classification utilizing full-genome RABV sequences. Because AP has the advantage that it is computationally fast and works for any meaningful measure of similarity between data samples, it has previously been applied successfully in bioinformatics, for analysis of microarray and gene expression data, however, cluster analysis of sequences is still in its infancy. Existing (516) and original (46) full genome RABV sequences were used to demonstrate the application of AP for RABV clustering. On a global scale, AP proposed four clusters, i.e. New World cluster, Arctic/Arctic-like, Cosmopolitan, and Asian as previously assigned by phylogenetic studies. By combining AP with established phylogenetic analyses, it is possible to resolve phylogenetic relationships between verifiably determined clusters and sequences. This workflow will be useful in confirming cluster distributions in a uniform transparent manner, not only for RABV, but also for other comparative sequence analyses.

Rabies is one of the oldest known zoonoses, caused by lyssaviruses. In recent years, more than 21,000 nucleotide sequences for rabies viruses (RABV) have been deposited in public databases. In this study, a novel mathematical approach called affinity propagation (AP) clustering, a highly powerful tool, to verifiably divide full genome RABV sequences into genetic clusters, was used. A panel of existing and novel RABV full genome sequences was used to demonstrate the application of AP for RABV clustering. Using a combination of AP with established phylogenetic analyses is useful in resolving phylogenetic relationships between more objectively determined clusters and sequences. This workflow will help to substantiate a transparent cluster distribution, not only for RABV, but also for other comparative sequence analyses.

Genotypes of hepatitis a virus in Turkey: first report and clinical profile of children infected with sub-genotypes IA and IIIA

Background

Hepatitis A virus (HAV) is a food and water-borne virus causing clinical (mainly hepatitis) and subclinical disease in humans. It is important to characterize circulating strains of HAV in order to prevent HAV infections using efficacious vaccines. The aim of this study was the detection and characterization of the circulating strains of HAV in Turkey by performing serology, RT-PCR, sequencing and phylogenetic analysis.

Methods

In this study, 355 HAV suspected cases were analysed by ELISA for the presence of antibodies to HAV. RNA was extracted from 54 HAV IgM positive human sera. None of the suspect cases were vaccinated against HAV and they never received blood transfusions. Samples found positive by RT-PCR using primers targeting the VP1/VP2A junction and VP1/VP3 capsid region of HAV, were subjected to sequencing and phylogenetic analyses.

Results

IgM type antibodies to HAV were detected in 54 patients. Twenty one of them were students. The age of IgM positive cases was between 3 and 60 years. IgM positivity differed in age groups and was higher in the age group 3 to 10 years. Phylogenetic analysis showed that the majority of HAV strains detected in this study belong to the “HAV 1B” cluster. In addition, the HAV sub-genotypes IA (KT874461.1) and IIIA (KT222963.1) were found in 2 children. These sub-genotypes were not previously reported in Turkey. The child who carried sub-genotype IIIA travelled to Afghanistan and presented with abdominal pain, icterus and vomitus. He was positive for anti-HAV IgM and IgG but negative for hepatitis B and C. Liver enzymes like aspartate aminotransferase, alanine aminotransferase, alkaline phosphatase, gamma-glutamyl transferase and lactate dehydrogenase were severely elevated. Bilirubin levels were also increased. White blood cells, neutrophils and hemoglobin were decreased while lymphocytes and monocytes were increased. Similar clinical signs and laboratory findings were reported for the child infected with sub-genotype IA but aspartate aminotransferase and alanine aminotransferase were not severely elevated.

Conclusions

The results indicate that molecular studies determining the HAV genotype variation in Turkey are timely and warranted. The majority of IgM positive cases in 3–10 year old patients indicate that childhood vaccination is important. Sub-genotype IB is the most prevalant genotype in Turkey. Surprisingly, sub-genotype IA and IIIA are also present in Turkey; future diagnostic efforts need to include diagnostic methods which can identify this emerging HAV genotypes. Our results also show that one important risk factor for contracting hepatitis A virus is international travel since genotype IIIA was detected in a child who had travelled to Afghanistan.

A machine learning approach for viral genome classification

Background

Advances in cloning and sequencing technology are yielding a massive number of viral genomes. The classification and annotation of these genomes constitute important assets in the discovery of genomic variability, taxonomic characteristics and disease mechanisms. Existing classification methods are often designed for specific well-studied family of viruses. Thus, the viral comparative genomic studies could benefit from more generic, fast and accurate tools for classifying and typing newly sequenced strains of diverse virus families.

Results

Here, we introduce a virus classification platform, CASTOR, based on machine learning methods. CASTOR is inspired by a well-known technique in molecular biology: restriction fragment length polymorphism (RFLP). It simulates, in silico, the restriction digestion of genomic material by different enzymes into fragments. It uses two metrics to construct feature vectors for machine learning algorithms in the classification step. We benchmark CASTOR for the classification of distinct datasets of human papillomaviruses (HPV), hepatitis B viruses (HBV) and human immunodeficiency viruses type 1 (HIV-1). Results reveal true positive rates of 99%, 99% and 98% for HPV Alpha species, HBV genotyping and HIV-1 M subtyping, respectively. Furthermore, CASTOR shows a competitive performance compared to well-known HIV-1 specific classifiers (REGA and COMET) on whole genomes and pol fragments.

Conclusion

The performance of CASTOR, its genericity and robustness could permit to perform novel and accurate large scale virus studies. The CASTOR web platform provides an open access, collaborative and reproducible machine learning classifiers. CASTOR can be accessed at http://castor.bioinfo.uqam.ca.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1602-3) contains supplementary material, which is available to authorized users.

Domain Organization and Evolution of the Highly Divergent 5' Coding Region of Genomes of Arteriviruses, Including the Novel Possum Nidovirus

In five experimentally characterized arterivirus species, the 5'-end genome coding region encodes the most divergent nonstructural proteins (nsp's), nsp1 and nsp2, which include papain-like proteases (PLPs) and other poorly characterized domains. These are involved in regulation of transcription, polyprotein processing, and virus-host interaction. Here we present results of a bioinformatics analysis of this region of 14 arterivirus species, including that of the most distantly related virus, wobbly possum disease virus (WPDV), determined by a modified 5' rapid amplification of cDNA ends (RACE) protocol. By combining profile-profile comparisons and phylogeny reconstruction, we identified an association of the four distinct domain layouts of nsp1-nsp2 with major phylogenetic lineages, implicating domain gain, including duplication, and loss in the early nsp1 evolution. Specifically, WPDV encodes highly divergent homologs of PLP1a, PLP1b, PLP1c, and PLP2, with PLP1a lacking the catalytic Cys residue, but does not encode nsp1 Zn finger (ZnF) and "nuclease" domains, which are conserved in other arteriviruses. Unexpectedly, our analysis revealed that the only catalytically active nsp1 PLP of equine arteritis virus (EAV), known as PLP1b, is most similar to PLP1c and thus is likely to be a PLP1b paralog. In all non-WPDV arteriviruses, PLP1b/c and PLP1a show contrasting patterns of conservation, with the N- and C-terminal subdomains, respectively, being enriched with conserved residues, which is indicative of different functional specializations. The least conserved domain of nsp2, the hypervariable region (HVR), has its size varied 5-fold and includes up to four copies of a novel PxPxPR motif that is potentially recognized by SH3 domain-containing proteins. Apparently, only EAV lacks the signal that directs -2 ribosomal frameshifting in the nsp2 coding region.IMPORTANCE Arteriviruses comprise a family of mammalian enveloped positive-strand RNA viruses that include some of the most economically important pathogens of swine. Most of our knowledge about this family has been obtained through characterization of viruses from five species: Equine arteritis virus, Simian hemorrhagic fever virus, Lactate dehydrogenase-elevating virus, Porcine respiratory and reproductive syndrome virus 1, and Porcine respiratory and reproductive syndrome virus 2 Here we present the results of comparative genomics analyses of viruses from all known 14 arterivirus species, including the most distantly related virus, WPDV, whose genome sequence was completed in this study. Our analysis focused on the multifunctional 5'-end genome coding region that encodes multidomain nonstructural proteins 1 and 2. Using diverse bioinformatics techniques, we identified many patterns of evolutionary conservation that are specific to members of distinct arterivirus species, both characterized and novel, or their groups. They are likely associated with structural and functional determinants important for virus replication and virus-host interaction.

Consensus statement: Virus taxonomy in the age of metagenomics

The number and diversity of viral sequences that are identified in metagenomic data far exceeds that of experimentally characterized virus isolates. In a recent workshop, a panel of experts discussed the proposal that, with appropriate quality control, viruses that are known only from metagenomic data can, and should be, incorporated into the official classification scheme of the International Committee on Taxonomy of Viruses (ICTV). Although a taxonomy that is based on metagenomic sequence data alone represents a substantial departure from the traditional reliance on phenotypic properties, the development of a robust framework for sequence-based virus taxonomy is indispensable for the comprehensive characterization of the global virome. In this Consensus Statement article, we consider the rationale for why metagenomic sequence data should, and how it can, be incorporated into the ICTV taxonomy, and present proposals that have been endorsed by the Executive Committee of the ICTV.

A taxonomy update for the family Polyomaviridae

Many distinct polyomaviruses infecting a variety of vertebrate hosts have recently been discovered, and their complete genome sequence could often be determined. To accommodate this fast-growing diversity, the International Committee on Taxonomy of Viruses (ICTV) Polyomaviridae Study Group designed a host- and sequence-based rationale for an updated taxonomy of the family Polyomaviridae. Applying this resulted in numerous recommendations of taxonomical revisions, which were accepted by the Executive Committee of the ICTV in December 2015. New criteria for definition and creation of polyomavirus species were established that were based on the observed distance between large T antigen coding sequences. Four genera (Alpha-, Beta, Gamma- and Deltapolyomavirus) were delineated that together include 73 species. Species naming was made as systematic as possible - most species names now consist of the binomial name of the host species followed by polyomavirus and a number reflecting the order of discovery. It is hoped that this important update of the family taxonomy will serve as a stable basis for future taxonomical developments.

Past, present, and future of arenavirus taxonomy

Until recently, members of the monogeneric family Arenaviridae (arenaviruses) have been known to infect only muroid rodents and, in one case, possibly phyllostomid bats. The paradigm of arenaviruses exclusively infecting small mammals shifted dramatically when several groups independently published the detection and isolation of a divergent group of arenaviruses in captive alethinophidian snakes. Preliminary phylogenetic analyses suggest that these reptilian arenaviruses constitute a sister clade to mammalian arenaviruses. Here, the members of the International Committee on Taxonomy of Viruses (ICTV) Arenaviridae Study Group, together with other experts, outline the taxonomic reorganization of the family Arenaviridae to accommodate reptilian arenaviruses and other recently discovered mammalian arenaviruses and to improve compliance with the Rules of the International Code of Virus Classification and Nomenclature (ICVCN). PAirwise Sequence Comparison (PASC) of arenavirus genomes and NP amino acid pairwise distances support the modification of the present classification. As a result, the current genus Arenavirus is replaced by two genera, Mammarenavirus and Reptarenavirus, which are established to accommodate mammalian and reptilian arenaviruses, respectively, in the same family. The current species landscape among mammalian arenaviruses is upheld, with two new species added for Lunk and Merino Walk viruses and minor corrections to the spelling of some names. The published snake arenaviruses are distributed among three new separate reptarenavirus species. Finally, a non-Latinized binomial species name scheme is adopted for all arenavirus species. In addition, the current virus abbreviations have been evaluated, and some changes are introduced to unequivocally identify each virus in electronic databases, manuscripts, and oral proceedings.

Discovery of an essential nucleotidylating activity associated with a newly delineated conserved domain in the RNA polymerase-containing protein of all nidoviruses

RNA viruses encode an RNA-dependent RNA polymerase (RdRp) that catalyzes the synthesis of their RNA(s). In the case of positive-stranded RNA viruses belonging to the order Nidovirales, the RdRp resides in a replicase subunit that is unusually large. Bioinformatics analysis of this non-structural protein has now revealed a nidoviral signature domain (genetic marker) that is N-terminally adjacent to the RdRp and has no apparent homologs elsewhere. Based on its conservation profile, this domain is proposed to have nucleotidylation activity. We used recombinant non-structural protein 9 of the arterivirus equine arteritis virus (EAV) and different biochemical assays, including irreversible labeling with a GTP analog followed by a proteomics analysis, to demonstrate the manganese-dependent covalent binding of guanosine and uridine phosphates to a lysine/histidine residue. Most likely this was the invariant lysine of the newly identified domain, named nidovirus RdRp-associated nucleotidyltransferase (NiRAN), whose substitution with alanine severely diminished the described binding. Furthermore, this mutation crippled EAV and prevented the replication of severe acute respiratory syndrome coronavirus (SARS-CoV) in cell culture, indicating that NiRAN is essential for nidoviruses. Potential functions supported by NiRAN may include nucleic acid ligation, mRNA capping and protein-primed RNA synthesis, possibilities that remain to be explored in future studies.

Methods for virus classification and the challenge of incorporating metagenomic sequence data

The division of viruses into orders, families, genera and species provides a classification framework that seeks to organize and make sense of the diversity of viruses infecting animals, plants and bacteria. Classifications are based on similarities in genome structure and organization, the presence of homologous genes and sequence motifs and at lower levels such as species, host range, nucleotide and antigenic relatedness and epidemiology. Classification below the level of family must also be consistent with phylogeny and virus evolutionary histories. Recently developed methods such as PASC, DEMaRC and NVR offer alternative strategies for genus and species assignments that are based purely on degrees of divergence between genome sequences. They offer the possibility of automating classification of the vast number of novel virus sequences being generated by next-generation metagenomic sequencing. However, distance-based methods struggle to deal with the complex evolutionary history of virus genomes that are shuffled by recombination and reassortment, and where taxonomic lineages evolve at different rates. In biological terms, classifications based on sequence distances alone are also arbitrary whereas the current system of virus taxonomy is of utility precisely because it is primarily based upon phenotypic characteristics. However, a separate system is clearly needed by which virus variants that lack biological information might be incorporated into the ICTV classification even if based solely on sequence relationships to existing taxa. For these, simplified taxonomic proposals and naming conventions represent a practical way to expand the existing virus classification and catalogue our rapidly increasing knowledge of virus diversity.

Arterivirus nsp12 versus the coronavirus nsp16 2'-O-methyltransferase: comparison of the C-terminal cleavage products of two nidovirus pp1ab polyproteins

The 3'-terminal domain of the most conserved ORF1b in three of the four families of the order Nidovirales (except for the family Arteriviridae) encodes a (putative) 2'-O-methyltransferase (2'-O-MTase), known as non structural protein (nsp) 16 in the family Coronaviridae and implicated in methylation of the 5' cap structure of nidoviral mRNAs. As with coronavirus transcripts, arterivirus mRNAs are assumed to possess a 5' cap although no candidate MTases have been identified thus far. To address this knowledge gap, we analysed the uncharacterized nsp12 of arteriviruses, which occupies the ORF1b position equivalent to that of the nidovirus 2'-O-MTase (coronavirus nsp16). In our in-depth bioinformatics analysis of nsp12, the protein was confirmed to be family specific whilst having diverged much further than other nidovirus ORF1b-encoded proteins, including those of the family Coronaviridae. Only one invariant and several partially conserved, predominantly aromatic residues were identified in nsp12, which may adopt a structure with alternating α-helices and β-strands, an organization also found in known MTases. However, no statistically significant similarity was found between nsp12 and the twofold larger coronavirus nsp16, nor could we detect MTase activity in biochemical assays using recombinant equine arteritis virus (EAV) nsp12. Our further analysis established that this subunit is essential for replication of this prototypic arterivirus. Using reverse genetics, we assessed the impact of 25 substitutions at 14 positions, yielding virus phenotypes ranging from WT-like to non-viable. Notably, replacement of the invariant phenylalanine 109 with tyrosine was lethal. We concluded that nsp12 plays an essential role during EAV replication, possibly by acting as a co-factor for another enzyme.

Similarity-based codes sequentially assigned to ebolavirus genomes are informative of species membership, associated outbreaks, and transmission chains

Genome-similarity based codes were assigned to individual ebolavirus isolates. Codes were found to be informative of phylogenetic and epidemiological relationships. It is proposed that such codes should be assigned to every genome-sequenced virus to complement current viral taxonomy.

Background. Developing a universal standardized microbial typing and nomenclature system that provides phylogenetic and epidemiological information in real time has never been as urgent in public health as it is today. We previously proposed to use genome similarity as the basis for immediate and precise typing and naming of individual organisms or viruses. In this study, we tested the validity of the proposed system and applied it to the epidemiology of infectious diseases using Ebola virus disease (EVD) outbreaks as the example.

Methods. One hundred twenty-eight publicly available ebolavirus genomes were compared with each other, and average nucleotide identity (ANI) was calculated. The ANI was then used to assign unique codes, hereafter referred to as Life Identification Numbers (LINs), to every viral isolate, whereby each LIN consisted of a series of positions reflecting increasing genome similarity. Congruence of LINs with phylogenetic and epidemiological relationships was then determined.

Results. Assigned LINs correlate with phylogeny at the species and infraspecies level and can even identify some individual transmission chains during the 2014–2015 EVD epidemic in West Africa.

Conclusions. Life Identification Numbers can provide a fast, automated, standardized, and scalable approach to precisely identify and name viral isolates upon genome sequence submission, facilitating unambiguous communication during disease epidemics among clinicians, epidemiologists, and governments.