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.

What's what in a pandemic? Virus, disease, and societal disaster must be differentiated.

Viruses, the diseases they can trigger, and the possible associated societal disaster represent different entities. To engage with the complexities of viral pandemics, we need to recognize each entity by using a distinctive name.

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 (Mpro). 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 Mn2+-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 Mpro-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.

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.

ICTV Virus Taxonomy Profile: Arteriviridae 2021.

The family Arteriviridae comprises enveloped RNA viruses with a linear, positive-sense genome of approximately 12.7 to 15.7 kb. The spherical, pleomorphic virions have a median diameter of 50-74 nm and include eight to eleven viral proteins. Arteriviruses infect non-human mammals in a vector-independent manner. Infections are often persistent and can either be asymptomatic or produce overt disease. Some arteriviruses are important veterinary pathogens while others infect particular species of wild rodents or African non-human primates. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Arteriviridae, which is available at ictv.global/report/arteriviridae.

A nidovirus perspective on SARS-CoV-2.

Two pandemics of respiratory distress diseases associated with zoonotic introductions of the species Severe acute respiratory syndrome-related coronavirus in the human population during 21st century raised unprecedented interest in coronavirus research and assigned it unseen urgency. The two viruses responsible for the outbreaks, SARS-CoV and SARS-CoV-2, respectively, are in the spotlight, and SARS-CoV-2 is the focus of the current fast-paced research. Its foundation was laid down by studies of many corona- and related viruses that collectively form the vast order Nidovirales. Comparative genomics of nidoviruses played a key role in this advancement over more than 30 years. It facilitated the transfer of knowledge from characterized to newly identified viruses, including SARS-CoV and SARS-CoV-2, as well as contributed to the dissection of the nidovirus proteome and identification of patterns of variations between different taxonomic groups, from species to families. This review revisits selected cases of protein conservation and variation that define nidoviruses, illustrates the remarkable plasticity of the proteome during nidovirus adaptation, and asks questions at the interface of the proteome and processes that are vital for nidovirus reproduction and could inform the ongoing research of SARS-CoV-2.

LAMPA, LArge Multidomain Protein Annotator, and its application to RNA virus polyproteins.

MOTIVATION: To facilitate accurate estimation of statistical significance of sequence similarity in profile-profile searches, queries should ideally correspond to protein domains. For multidomain proteins, using domains as queries depends on delineation of domain borders, which may be unknown. Thus, proteins are commonly used as queries that complicate establishing homology for similarities close to cutoff levels of statistical significance. RESULTS: In this article, we describe an iterative approach, called LAMPA, LArge Multidomain Protein Annotator, that resolves the above conundrum by gradual expansion of hit coverage of multidomain proteins through re-evaluating statistical significance of hit similarity using ever smaller queries defined at each iteration. LAMPA employs TMHMM and HHsearch for recognition of transmembrane regions and homology, respectively. We used Pfam database for annotating 2985 multidomain proteins (polyproteins) composed of >1000 amino acid residues, which dominate proteomes of RNA viruses. Under strict cutoffs, LAMPA outperformed HHsearch-mediated runs using intact polyproteins as queries by three measures: number of and coverage by identified homologous regions, and number of hit Pfam profiles. Compared to HHsearch, LAMPA identified 507 extra homologous regions in 14.4% of polyproteins. This Pfam-based annotation of RNA virus polyproteins by LAMPA was also superior to RefSeq expert annotation by two measures, region number and annotated length, for 69.3% of RNA virus polyprotein entries. We rationalized the obtained results based on dependencies of HHsearch hit statistical significance for local alignment similarity score from lengths and diversities of query-target pairs in computational experiments. AVAILABILITY AND IMPLEMENTATION: LAMPA 1.0.0 R package is placed at github (https://github.com/Gorbalenya-Lab/LAMPA). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

Additional changes to taxonomy ratified in a special vote by the International Committee on Taxonomy of Viruses (October 2018).

This article reports the changes to virus taxonomy approved and ratified by the International Committee on Taxonomy of Viruses (ICTV) in October 2018. Of note, the ICTV has approved, by an absolute majority, the creation of additional taxonomical ranks above those recognized previously. A total of 15 ranks (realm, subrealm, kingdom, subkingdom, phylum, subphylum, class, subclass, order, suborder, family, subfamily, genus, subgenus, and species) are now available to encompass the entire spectrum of virus diversity. Classification at ranks above genus is not obligatory but can be used by the authors of new taxonomic proposals when scientific justification is provided.

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.

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.

Changes to taxonomy and the International Code of Virus Classification and Nomenclature ratified by the International Committee on Taxonomy of Viruses (2018).

This article lists the changes to virus taxonomy approved and ratified by the International Committee on Taxonomy of Viruses in February 2018. A total of 451 species, 69 genera, 11 subfamilies, 9 families and one new order were added to the taxonomy. The current totals at each taxonomic level now stand at 9 orders, 131 families, 46 subfamilies, 803 genera and 4853 species. A change was made to the International Code of Virus Classification and Nomenclature to allow the use of the names of people in taxon names under appropriate circumstances. An updated Master Species List incorporating the approved changes was released in March 2018 ( https://talk.ictvonline.org/taxonomy/ ).

Modulation of the Host Lipid Landscape to Promote RNA Virus Replication: The Picornavirus Encephalomyocarditis Virus Converges on the Pathway Used by Hepatitis C Virus.

Cardioviruses, including encephalomyocarditis virus (EMCV) and the human Saffold virus, are small non-enveloped viruses belonging to the Picornaviridae, a large family of positive-sense RNA [(+)RNA] viruses. All (+)RNA viruses remodel intracellular membranes into unique structures for viral genome replication. Accumulating evidence suggests that picornaviruses from different genera use different strategies to generate viral replication organelles (ROs). For instance, enteroviruses (e.g. poliovirus, coxsackievirus, rhinovirus) rely on the Golgi-localized phosphatidylinositol 4-kinase III beta (PI4KB), while cardioviruses replicate independently of the kinase. By which mechanisms cardioviruses develop their ROs is currently unknown. Here we show that cardioviruses manipulate another PI4K, namely the ER-localized phosphatidylinositol 4-kinase III alpha (PI4KA), to generate PI4P-enriched ROs. By siRNA-mediated knockdown and pharmacological inhibition, we demonstrate that PI4KA is an essential host factor for EMCV genome replication. We reveal that the EMCV nonstructural protein 3A interacts with and is responsible for PI4KA recruitment to viral ROs. The ensuing phosphatidylinositol 4-phosphate (PI4P) proved important for the recruitment of oxysterol-binding protein (OSBP), which delivers cholesterol to EMCV ROs in a PI4P-dependent manner. PI4P lipids and cholesterol are shown to be required for the global organization of the ROs and for viral genome replication. Consistently, inhibition of OSBP expression or function efficiently blocked EMCV RNA replication. In conclusion, we describe for the first time a cellular pathway involved in the biogenesis of cardiovirus ROs. Remarkably, the same pathway was reported to promote formation of the replication sites of hepatitis C virus, a member of the Flaviviridae family, but not other picornaviruses or flaviviruses. Thus, our results highlight the convergent recruitment by distantly related (+)RNA viruses of a host lipid-modifying pathway underlying formation of viral replication sites.

Changes to taxonomy and the International Code of Virus Classification and Nomenclature ratified by the International Committee on Taxonomy of Viruses (2017).

This article lists the changes to virus taxonomy approved and ratified by the International Committee on Taxonomy of Viruses (ICTV) in March 2017.

50 years of the International Committee on Taxonomy of Viruses: progress and prospects.

We mark the 50th anniversary of the International Committee on Taxonomy of Viruses (ICTV) by presenting a brief history of the organization since its foundation, showing how it has adapted to advancements in our knowledge of virus diversity and the methods used to characterize it. We also outline recent developments, supported by a grant from the Wellcome Trust (UK), that are facilitating substantial changes in the operations of the ICTV and promoting dialogue with the virology community. These developments will generate improved online resources, including a freely available and regularly updated ICTV Virus Taxonomy Report. They also include a series of meetings between the ICTV and the broader community focused on some of the major challenges facing virus taxonomy, with the outcomes helping to inform the future policy and practice of the ICTV.

Interspecific adaptation by binary choice at de novo polyomavirus T antigen site through accelerated codon-constrained Val-Ala toggling within an intrinsically disordered region.

It is common knowledge that conserved residues evolve slowly. We challenge generality of this central tenet of molecular biology by describing the fast evolution of a conserved nucleotide position that is located in the overlap of two open reading frames (ORFs) of polyomaviruses. The de novo ORF is expressed through either the ALTO protein or the Middle T antigen (MT/ALTO), while the ancestral ORF encodes the N-terminal domain of helicase-containing Large T (LT) antigen. In the latter domain the conserved Cys codon of the LXCXE pRB-binding motif constrains codon evolution in the overlapping MT/ALTO ORF to a binary choice between Val and Ala codons, termed here as codon-constrained Val-Ala (COCO-VA) toggling. We found the rate of COCO-VA toggling to approach the speciation rate and to be significantly accelerated compared to the baseline rate of chance substitution in a large monophyletic lineage including all viruses encoding MT/ALTO and three others. Importantly, the COCO-VA site is located in a short linear motif (SLiM) of an intrinsically disordered region, a typical characteristic of adaptive responders. These findings provide evidence that the COCO-VA toggling is under positive selection in many polyomaviruses, implying its critical role in interspecific adaptation, which is unprecedented for conserved residues.

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.

One severe acute respiratory syndrome coronavirus protein complex integrates processive RNA polymerase and exonuclease activities.

In addition to members causing milder human infections, the Coronaviridae family includes potentially lethal zoonotic agents causing severe acute respiratory syndrome (SARS) and the recently emerged Middle East respiratory syndrome. The ∼30-kb positive-stranded RNA genome of coronaviruses encodes a replication/transcription machinery that is unusually complex and composed of 16 nonstructural proteins (nsps). SARS-CoV nsp12, the canonical RNA-dependent RNA polymerase (RdRp), exhibits poorly processive RNA synthesis in vitro, at odds with the efficient replication of a very large RNA genome in vivo. Here, we report that SARS-CoV nsp7 and nsp8 activate and confer processivity to the RNA-synthesizing activity of nsp12. Using biochemical assays and reverse genetics, the importance of conserved nsp7 and nsp8 residues was probed. Whereas several nsp7 mutations affected virus replication to a limited extent, the replacement of two nsp8 residues (P183 and R190) essential for interaction with nsp12 and a third (K58) critical for the interaction of the polymerase complex with RNA were all lethal to the virus. Without a loss of processivity, the nsp7/nsp8/nsp12 complex can associate with nsp14, a bifunctional enzyme bearing 3'-5' exoribonuclease and RNA cap N7-guanine methyltransferase activities involved in replication fidelity and 5'-RNA capping, respectively. The identification of this tripartite polymerase complex that in turn associates with the nsp14 proofreading enzyme sheds light on how coronaviruses assemble an RNA-synthesizing machinery to replicate the largest known RNA genomes. This protein complex is a fascinating example of the functional integration of RNA polymerase, capping, and proofreading activities.