Reply to Holmes and Duchêne, "Can Sequence Phylogenies Safely Infer the Origin of the Global Virome?": Deep Phylogenetic Analysis of RNA Viruses Is Highly Challenging but Not Meaningless

Over two decades of research on the marine RNA virosphere

RNA viruses (realm: Riboviria), including RNA phages and eukaryote-infecting RNA viruses, are essential components of marine ecosystems. A large number of marine RNA viruses have been discovered in the last two decades because of the rapid development of next-generation sequencing (NGS) technology. Indeed, the combination of NGS and state-of-the-art meta-omics methods (viromics, the study of all viruses in a specific environment) has led to a fundamental understanding of the taxonomy and genetic diversity of RNA viruses in the sea, suggesting the complex ecological roles played by RNA viruses in this complex ecosystem. Furthermore, comparisons of viromes in the context of highly variable marine niches reveal the biogeographic patterns and ecological impact of marine RNA viruses, whose role in global ecology is becoming increasingly clearer. In this review, we summarize the characteristics of the global marine RNA virosphere and outline the taxonomic hierarchy of RNA viruses with a specific focus on their ancient evolutionary history. We also review the development of methodology and the major progress resulting from its applications in RNA viromics. The aim of this review is not only to provide an in-depth understanding of multifaceted aspects of marine RNA viruses, but to offer future perspectives on developing a better methodology for discovery, and exploring the evolutionary origin and major ecological significance of marine RNA virosphere.

The origins of COVID-19 pandemic: A brief overview

The novel coronavirus disease (COVID-19) outbreak that emerged at the end of 2019 has now swept the world for more than 2 years, causing immeasurable damage to the lives and economies of the world. It has drawn so much attention to discovering how the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) originated and entered the human body. The current argument revolves around two contradictory theories: a scenario of laboratory spillover events and human contact with zoonotic diseases. Here, we reviewed the transmission, pathogenesis, possible hosts, as well as the genome and protein structure of SARS-CoV-2, which play key roles in the COVID-19 pandemic. We believe the coronavirus was originally transmitted to human by animals rather than by a laboratory leak. However, there still needs more investigations to determine the source of the pandemic. Understanding how COVID-19 emerged is vital to developing global strategies for mitigating future outbreaks.

A discussion of RNA virus taxonomy based on the 2020 International Committee on Taxonomy of Viruses report

RNA viruses have a higher mutation rate than DNA viruses; however, RNA viruses are insufficiently studied outside disease settings. The International Committee on Taxonomy of Viruses (ICTV) is an organization set up by virologists to standardize virus classification. To better understand ICTV taxonomy and the characteristics and rules of different RNA virus families, we analyzed the 3,529 RNA viruses included in the 2020 ICTV report using five widely used metrics: length, host, GC content, number of predicted ORFs, and sequence similarity. The results show that host type has a significant influence on viral genome length and GC content. The genome lengths of virus members within the same genus are quite similar: 98.28% of the genome length differences within any particular genus are less than 20%. The species within those genera containing segmented viruses also have a similar length and number of segments. The number of predicted ORFs in the RNA viral genomes also shows a strong, statistically significant correlation with genome length. We suggest that due to the high mutation rate of RNA virus genomes, current RNA virus classification should mainly rely on protein similarities rather than nucleic acid similarities.

NeoRdRp: A Comprehensive Dataset for Identifying RNA-dependent RNA Polymerases of Various RNA Viruses from Metatranscriptomic Data

RNA viruses are distributed throughout various environments, and most have recently been identified by metatranscriptome sequencing. However, due to the high nucleotide diversity of RNA viruses, it is still challenging to identify novel RNA viruses from metatranscriptome data. To overcome this issue, we created a dataset of RNA-dependent RNA polymerase (RdRp) domains that are essential for all RNA viruses belonging to Orthornavirae. Genes with RdRp domains from various RNA viruses were clustered based on amino acid sequence similarities. A multiple sequence alignment was generated for each cluster, and a hidden Markov model (HMM) profile was created when the number of sequences was greater than three. We further refined 426 HMM profiles by detecting RefSeq RNA virus sequences and subsequently combined the hit sequences with the RdRp domains. As a result, 1,182 HMM profiles were generated from 12,502 RdRp domain sequences, and the dataset was named NeoRdRp. The majority of NeoRdRp HMM profiles successfully detected RdRp domains, specifically in the UniProt dataset. Furthermore, we compared the NeoRdRp dataset with two previously reported methods for RNA virus detection using metatranscriptome sequencing data. Our methods successfully identified the majority of RNA viruses in the datasets; however, some RNA viruses were not detected, similar to the other two methods. NeoRdRp may be repeatedly improved by the addition of new RdRp sequences and is applicable as a system for detecting various RNA viruses from diverse metatranscriptome data.

Revised Taxonomy of Rhabdoviruses Infecting Fish and Marine Mammals

Simple Summary

The Rhabdoviridae is a family of viruses that includes some important pathogens of fish and marine mammals. Aspects of the taxonomic classification of fish viruses assigned to this family have recently been reviewed by the International Committee on Taxonomy of Viruses (ICTV). This paper describes the newly approved taxonomy, including the assignment of new subfamilies and new virus species. The paper also considers a taxonomic conundrum presented by viruses assigned to one group of fish rhabdoviruses (genus Novirhabdovirus) for which assignment to the family Rhabdoviridae may not be appropriate.

Abstract

The Rhabdoviridae is a large family of negative-sense (-) RNA viruses that includes important pathogens of ray-finned fish and marine mammals. As for all viruses, the taxonomic assignment of rhabdoviruses occurs through a process implemented by the International Committee on Taxonomy of Viruses (ICTV). A recent revision of taxonomy conducted in conjunction with the ICTV Rhabdoviridae Study Group has resulted in the establishment of three new subfamilies (Alpharhabdovirinae, Betarhabdovirinae, and Gammarhabdovirinae) within the Rhabdoviridae, as well as three new genera (Cetarhavirus, Siniperhavirus, and Scophrhavirus) and seven new species for viruses infecting fish or marine mammals. All rhabdovirus species have also now been named or renamed to comply with the binomial format adopted by the ICTV in 2021, comprising the genus name followed by a species epithet. Phylogenetic analyses of L protein (RNA-dependent RNA polymerase) sequences of (-) RNA viruses indicate that members of the genus Novirhabdovirus (subfamily Gammarhabdovirinae) do not cluster within the Rhabdoviridae, suggesting the need for a review of their current classification.

Cryptic and abundant marine viruses at the evolutionary origins of Earth's RNA virome

Whereas DNA viruses are known to be abundant, diverse, and commonly key ecosystem players, RNA viruses are insufficiently studied outside disease settings. In this study, we analyzed ≈28 terabases of Global Ocean RNA sequences to expand Earth's RNA virus catalogs and their taxonomy, investigate their evolutionary origins, and assess their marine biogeography from pole to pole. Using new approaches to optimize discovery and classification, we identified RNA viruses that necessitate substantive revisions of taxonomy (doubling phyla and adding >50% new classes) and evolutionary understanding. "Species"-rank abundance determination revealed that viruses of the new phyla "Taraviricota," a missing link in early RNA virus evolution, and "Arctiviricota" are widespread and dominant in the oceans. These efforts provide foundational knowledge critical to integrating RNA viruses into ecological and epidemiological models.

Analysis of the virome associated to grapevine downy mildew lesions reveals new mycovirus lineages

The obligatory biotrophic oomycetes Plasmopara viticola is the causal agent of downy mildew, a destructive disease of grapevine worldwide. So far, chemical fungicides are widely employed to limit this pathogen, but their adverse effects are stimulating the quest for environmentally friendly alternative approaches. Here, we report on the search for mycoviruses that might be later developed as biocontrol agents for this pathogen. Symptomatic leaves were collected from various regions in Spain and Italy and mycelia associated to leaf lesions was harvested. Total RNA extractions were depleted of rRNA and metatranscriptomes were generated using a high-throughput sequencing approach. The virome associated to leaf lesions was then characterized through a bioinformatic pipeline relying on blast searches against current viral databases. Here, we present an inventory of 283 new RNA viruses: 222 positive strand RNA viruses, 29 negative strand RNA viruses, 27 double-stranded RNA viruses and 5 ORFan virus RdRP segments, which could not be reliably assigned to any existing group in the Riboviria. In addition to ORFan viruses, we found other surprising new evolutionary trajectories in this wide inventory of viruses. The most represented viruses in our collection are those in phylum Lenarviricota, and, among them, a group of mycovirus segments distantly related to narnaviruses, but characterized by a polymerase palm domain lacking subdomain C, with the putative GDD catalytic triad. We also provided evidence of a strict association between two RNA segments that form a new mycovirus clade of positive strand RNA in the phylum Kitrinoviricota, order Martellivirales. In the phylum Negarnaviricota, we report for the first time in the order Mononegavirales a clade of viruses that is ambisense, a feature that so far was present only in the order Bunyavirales. Furthermore, in the same phylum we detected the widespread occurrence and abundant accumulation in our libraries of a distinct mycovirus clade distantly related to the Muvirales and Goujanvirales orders, which so far include only viruses infecting invertebrates. Possible new oomycetes-specific virus clades are also described in the phylum Duplornaviricota. These data greatly expand the evolutionary history of mycoviruses adding new layers of diversity to the realm Riboviria.

Changes to virus taxonomy and the Statutes ratified by the International Committee on Taxonomy of Viruses (2020)

This article reports the changes to virus classification and taxonomy approved and ratified by the International Committee on Taxonomy of Viruses (ICTV) in March 2020. The entire ICTV was invited to vote on 206 taxonomic proposals approved by the ICTV Executive Committee at its meeting in July 2019, as well as on the proposed revision of the ICTV Statutes. All proposals and the revision of the Statutes were approved by an absolute majority of the ICTV voting membership. Of note, ICTV has approved a proposal that extends the previously established realm Riboviria to encompass nearly all RNA viruses and reverse-transcribing viruses, and approved three separate proposals to establish three realms for viruses with DNA genomes.

Understanding evolution of SARS-CoV-2: A perspective from analysis of genetic diversity of RdRp gene

Coronavirus disease 2019 emerged as the first example of “Disease X”, a hypothetical disease of humans caused by an unknown infectious agent that was named as novel coronavirus and subsequently designated as severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2). The origin of the outbreak at the animal market in Wuhan, China implies it as a case of zoonotic spillover. The study was designed to understand evolution of Betacoronaviruses and in particular diversification of SARS‐CoV‐2 using RNA dependent RNA polymerase (RdRp) gene, a stable genetic marker. Phylogenetic and population stratification analyses were carried out using maximum likelihood and Bayesian methods, respectively. Molecular phylogeny using RdRp showed that SARS‐CoV‐2 isolates cluster together. Bat‐CoV isolate RaTG13 and Pangolin‐CoVs are observed to branch off prior to SARS‐CoV‐2 cluster. While SARS‐CoV form a single cluster, Bat‐CoVs form multiple clusters. Population‐based analyses revealed that both SARS‐CoV‐2 and SARS‐CoV form separate clusters with no admixture. Bat‐CoVs were found to have single and mixed ancestry and clustered as four sub‐populations. Population‐based analyses of Betacoronaviruses using RdRp revealed that SARS‐CoV‐2 is a homogeneous population. SARS‐CoV‐2 appears to have evolved from Bat‐CoV isolate RaTG13, which diversified from a common ancestor from which Pangolin‐CoVs have also evolved. The admixed Bat‐CoV sub‐populations indicate that bats serve as reservoirs harboring virus ensembles that are responsible for zoonotic spillovers such as SARS‐CoV and SARS‐CoV‐2. The extent of admixed isolates of Bat‐CoVs observed in population diversification studies underline the need for periodic surveillance of bats and other animal reservoirs for potential spillovers as a measure towards preparedness for emergence of zoonosis.

SARS‐CoV‐2 shares a common ancestor with Bat‐CoV (isolate RaTG13) and Pangolin‐CoVs.

Analysis based on RdRp gene revealed that SARS‐CoV‐2 form a distinct homogeneous population.

Bat‐CoVs are found to be admixed ensembles indicative of their potential for zoonotic spillovers.

Our study recommends use of RdRp gene for survelliance to track emerging RNA viruses.

Hadaka Virus 1: a Capsidless Eleven-Segmented Positive-Sense Single-Stranded RNA Virus from a Phytopathogenic Fungus, Fusarium oxysporum

The search for viruses infecting fungi, or mycoviruses, has extended our knowledge about the diversity of RNA viruses, as exemplified by the discovery of polymycoviruses, a phylogenetic group of multisegmented RNA viruses with unusual forms. The genomic RNAs of known polymycoviruses, which show a phylogenetic affinity for animal positive-sense single-stranded RNA [(+)RNA] viruses such as caliciviruses, are comprised of four conserved segments with an additional zero to four segments. The double-stranded form of polymycovirus genomic RNA is assumed to be associated with a virally encoded protein (proline-alanine-serine-rich protein [PASrp]) in either of two manners: a capsidless colloidal form or a filamentous encapsidated form. Detailed molecular characterizations of polymycoviruses, however, have been conducted for only a few strains. Here, a novel polymyco-related virus named Hadaka virus 1 (HadV1), from the phytopathogenic fungus Fusarium oxysporum, was characterized. The genomic RNA of HadV1 consisted of an 11-segmented positive-sense RNA with highly conserved terminal nucleotide sequences. HadV1 shared the three conserved segments with known polymycoviruses but lacked the PASrp-encoding segment. Unlike the known polymycoviruses and encapsidated viruses, HadV1 was not pelleted by conventional ultracentrifugation, possibly due to the lack of PASrp. This result implied that HadV1 exists only as a soluble form with naked RNA. Nevertheless, the 11 genomic segments of HadV1 have been stably maintained through host subculturing and conidiation. Taken together, the results of this study revealed a virus with a potential novel virus lifestyle, carrying many genomic segments without typical capsids or PASrp-associated forms.IMPORTANCE Fungi collectively host various RNA viruses. Examples include encapsidated double-stranded RNA (dsRNA) viruses with diverse numbers of genomic segments (from 1 to 12) and capsidless viruses with nonsegmented (+)RNA genomes. Recently, viruses with unusual intermediate features of an infectious entity between encapsidated dsRNA viruses and capsidless (+)RNA viruses were found. They are called polymycoviruses, which typically have four to eight dsRNA genomic segments associated with one of the virus-encoded proteins and are phylogenetically distantly related to animal (+)RNA caliciviruses. Here, we identified a novel virus phylogenetically related to polymycoviruses, from the phytopathogenic fungus Fusarium oxysporum The virus, termed Hadaka virus 1 (HadV1), has 11 (+)RNA genomic segments, the largest number in known (+)RNA viruses. Nevertheless, HadV1 lacked a typical structural protein of polymycoviruses and was not pelleted by standard ultracentrifugation, implying an unusual capsidless nature of HadV1. This study reveals a potential novel lifestyle of multisegmented RNA viruses.

Global Organization and Proposed Megataxonomy of the Virus World

Viruses and mobile genetic elements are molecular parasites or symbionts that coevolve with nearly all forms of cellular life. The route of virus replication and protein expression is determined by the viral genome type. Comparison of these routes led to the classification of viruses into seven "Baltimore classes" (BCs) that define the major features of virus reproduction. However, recent phylogenomic studies identified multiple evolutionary connections among viruses within each of the BCs as well as between different classes. Due to the modular organization of virus genomes, these relationships defy simple representation as lines of descent but rather form complex networks. Phylogenetic analyses of virus hallmark genes combined with analyses of gene-sharing networks show that replication modules of five BCs (three classes of RNA viruses and two classes of reverse-transcribing viruses) evolved from a common ancestor that encoded an RNA-directed RNA polymerase or a reverse transcriptase. Bona fide viruses evolved from this ancestor on multiple, independent occasions via the recruitment of distinct cellular proteins as capsid subunits and other structural components of virions. The single-stranded DNA (ssDNA) viruses are a polyphyletic class, with different groups evolving by recombination between rolling-circle-replicating plasmids, which contributed the replication protein, and positive-sense RNA viruses, which contributed the capsid protein. The double-stranded DNA (dsDNA) viruses are distributed among several large monophyletic groups and arose via the combination of distinct structural modules with equally diverse replication modules. Phylogenomic analyses reveal the finer structure of evolutionary connections among RNA viruses and reverse-transcribing viruses, ssDNA viruses, and large subsets of dsDNA viruses. Taken together, these analyses allow us to outline the global organization of the virus world. Here, we describe the key aspects of this organization and propose a comprehensive hierarchical taxonomy of viruses.

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

This article reports the changes to virus taxonomy approved and ratified by the International Committee on Taxonomy of Viruses (ICTV) in February 2019. Of note, in addition to seven new virus families, the ICTV has approved, by an absolute majority, the creation of the realm Riboviria, a likely monophyletic group encompassing all viruses with positive-strand, negative-strand and double-strand genomic RNA that use cognate RNA-directed RNA polymerases for replication.