A new word, phylodynamics, was coined to emphasize the interconnection between phylogenetic properties, as observed for instance in a phylogenetic tree, and the epidemic dynamics of viruses, where selection, mediated by the host immune response, and transmission play a crucial role. The challenges faced when investigating the evolution of RNA viruses call for a virtuous loop of data collection, data analysis and modeling. This already resulted both in the collection of massive sequences databases and in the formulation of hypotheses on the main mechanisms driving qualitative differences observed in the (reconstructed) evolutionary patterns of different RNA viruses. Qualitatively, it has been observed that selection driven by the host immune response induces an uneven survival ability among co-existing strains. As a consequence, the imbalance level of the phylogenetic tree is manifestly more pronounced if compared to the case when the interaction with the host immune system does not play a central role in the evolutive dynamics. While many imbalance metrics have been introduced, reliable methods to discriminate in a quantitative way different level of imbalance are still lacking. In our work, we reconstruct and analyze the phylogenetic trees of six RNA viruses, with a special emphasis on the human Influenza A virus, due to its relevance for vaccine preparation as well as for the theoretical challenges it poses due to its peculiar evolutionary dynamics. We focus in particular on topological properties. We point out the limitation featured by standard imbalance metrics, and we introduce a new methodology with which we assign the correct imbalance level of the phylogenetic trees, in agreement with the phylodynamics of the viruses. Our thorough quantitative analysis allows for a deeper understanding of the evolutionary dynamics of the considered RNA viruses, which is crucial in order to provide a valuable framework for a quantitative assessment of theoretical predictions.

Although arthropods are important viral vectors, the biodiversity of arthropod viruses, as well as the role that arthropods have played in viral origins and evolution, is unclear. Through RNA sequencing of 70 arthropod species we discovered 112 novel viruses that appear to be ancestral to much of the documented genetic diversity of negative-sense RNA viruses, a number of which are also present as endogenous genomic copies. With this greatly enriched diversity we revealed that arthropods contain viruses that fall basal to major virus groups, including the vertebrate-specific arenaviruses, filoviruses, hantaviruses, influenza viruses, lyssaviruses, and paramyxoviruses. We similarly documented a remarkable diversity of genome structures in arthropod viruses, including a putative circular form, that sheds new light on the evolution of genome organization. Hence, arthropods are a major reservoir of viral genetic diversity and have likely been central to viral evolution.

CRISPR RNA-guided endonucleases have enabled precise editing of DNA. However, options for editing RNA remain limited. Here, we combine sequence-specific RNA cleavage by CRISPR ribonucleases with programmable RNA repair to make precise deletions and insertions in RNA. This work establishes a new recombinant RNA technology with immediate applications for the facile engineering of RNA viruses.

CRISPR RNA-guided endonucleases have enabled precise editing of DNA. However, options for editing RNA remain limited. Here, we combine sequence-specific RNA cleavage by CRISPR ribonucleases with programmable RNA repair to make precise deletions and insertions in RNA. This work establishes a recombinant RNA technology with immediate applications for the facile engineering of RNA viruses.

Nudibranchia is an under-studied taxonomic group of gastropods, including more than 3,000 species with colourful and extravagant body shapes and peculiar predatory and defensive strategies. Although symbiosis with bacteria has been reported, no data are available for the nudibranch microbiome nor regarding viruses possibly associated with these geographically widespread species.

Viruses are obligatory intracellular pathogens and completely depend on their hosts for survival and reproduction. The strategies adopted by viruses to exploit host cell processes and to evade host immune systems during infections may differ largely with the type of the viral genetic material. An improved understanding of these viral infection mechanisms is only possible through a better understanding of the pathogen-host interactions (PHIs) that enable viruses to enter into the host cells and manipulate the cellular mechanisms to their own advantage. Experimentally-verified protein-protein interaction (PPI) data of pathogen-host systems only became available at large scale within the last decade. In this study, we comparatively analyzed the current PHI networks belonging to DNA and RNA viruses and their human host, to get insights into the infection strategies used by these viral groups. We investigated the functional properties of human proteins in the PHI networks, to observe and compare the attack strategies of DNA and RNA viruses. We observed that DNA viruses are able to attack both human cellular and metabolic processes simultaneously during infections. On the other hand, RNA viruses preferentially interact with human proteins functioning in specific cellular processes as well as in intracellular transport and localization within the cell. Observing virus-targeted human proteins, we propose heterogeneous nuclear ribonucleoproteins and transporter proteins as potential antiviral therapeutic targets. The observed common and specific infection mechanisms in terms of viral strategies to attack human proteins may provide crucial information for further design of broad and specific next-generation antiviral therapeutics.

Mutational robustness is defined as the constancy of a phenotype in the face of deleterious mutations. Whether robustness can be directly favored by natural selection remains controversial. Theory and in silico experiments predict that, at high mutation rates, slow-replicating genotypes can potentially outcompete faster counterparts if they benefit from a higher robustness. Here, we experimentally validate this hypothesis, dubbed the "survival of the flattest," using two populations of the vesicular stomatitis RNA virus. Characterization of fitness distributions and genetic variability indicated that one population showed a higher replication rate, whereas the other was more robust to mutation. The faster replicator outgrew its robust counterpart in standard competition assays, but the outcome was reversed in the presence of chemical mutagens. These results show that selection can directly favor mutational robustness and reveal a novel viral resistance mechanism against treatment by lethal mutagenesis.

RNA viruses are a major threat to human health. Here, based on extensive literature searches carried out over a period of 18 years, we provide a catalogue of all 214 known human-infective RNA virus species. We link these viruses to metadata for a number of traits that influence their epidemiology, including the date of the first report of human infection, transmissibility in human populations, transmission route(s) and host range. This database can be used in comparative studies of human-infective RNA viruses to identify the characteristics of viruses most likely to pose the greatest public health threat, both now and in the future.

The study of invertebrate-and particularly nematode-viruses is emerging with the advancement of transcriptome sequencing. Five single-stranded RNA viruses have now been confirmed within the economically important soybean cyst nematode (SCN; Heterodera glycines). From previous research, we know these viruses to be widespread in greenhouse and field populations of SCN. Several of the SCN viruses were also confirmed within clover (H. trifolii) and beet (H. schachtii) cyst nematodes. In the presented study, we sequenced the transcriptomes of several inbred SCN populations and identified two previously undiscovered viral-like genomes. Both of these proposed viruses are negative-sense RNA viruses and have been named SCN nyami-like virus (NLV) and SCN bunya-like virus (BLV). Finally, we analyzed publicly available transcriptome data of two potato cyst nematode (PCN) species, Globodera pallida and G. rostochiensis. From these data, a third potential virus was discovered and called PCN picorna-like virus (PLV). PCN PLV is a positive-sense RNA virus, and to the best of our knowledge, is the first virus described within PCN. The presence of these novel viruses was confirmed via qRT-PCR, endpoint PCR, and Sanger sequencing with the exception of PCN PLV due to quarantine restrictions on the nematode host. While much work needs to be done to understand the biological and evolutionary significance of these viruses, they offer insight into nematode ecology and the possibility of novel nematode management strategies.

Avian haemosporidian parasites can cause malaria-like symptoms in their hosts and have been implicated in the demise of some bird species. The newly described Matryoshka RNA viruses (MaRNAV1 and MaRNAV2) infect haemosporidian parasites that in turn infect their vertebrate hosts. MaRNAV2 was the first RNA virus discovered associated with parasites of the genus Leucocytozoon. By analyzing metatranscriptomes from the NCBI SRA database with local sequence alignment tools, we detected two novel RNA viruses; we describe them as MaRNAV3 associated with Leucocytozoon and MaRNAV4 associated with Parahaemoproteus. MaRNAV3 had ~59% amino acid identity to the RNA-dependent RNA-polymerase (RdRp) of MaRNAV1 and ~63% amino acid identity to MaRNAV2. MaRNAV4 had ~44% amino acid identity to MaRNAV1 and ~47% amino acid identity to MaRNAV2. These findings lead us to hypothesize that MaRNAV_like viruses are widespread and tightly associated with the order Haemosporida since they have been described in human Plasmodium vivax, and now two genera of avian haemosporidians.

CpG repression in RNA viruses has been known for decades, but a reasonable explanation has not yet been proposed to explain this phenomenon. In this study, we calculated the CpG odds ratio of all RNA viruses that have available genome sequences and analyzed the correlation with their genome polarity, base composition, synonymous codon usage, phylogenetic relationship, and host. The results indicated that the viral base composition, synonymous codon usage and host selection were the dominant factors that determined the CpG bias in RNA viruses. CpG usage variation between the different viral groups was caused by different combinations of these pressures, which also differed from each other in strength. The consistent under-representation of CpG usage in -ssRNA viruses is determined predominantly by base composition, which may be a consequence of the U/A preferred mutation bias of -ssRNA viruses, whereas the CpG usage of +ssRNA viruses is affected greatly by their hosts. As a result, most +ssRNA viruses mimic their hosts' CpG usage. Unbiased CpG usage in dsRNA viruses is most likely a result of their dsRNA genome, which allows the viruses to escape from the host-driven CpG elimination pressure. CpG was under-represented in all reverse-transcribing viruses (RT viruses), suggesting that DNA methylation is an important factor affecting the CpG usage of retroviruses. However, vertebrate-infecting RT viruses may also suffer host' CpG elimination pressure that also acts on +ssRNA viruses, which results in further under-representation of CpG in the vertebrate-infecting RT viruses.

RNA viruses are among the most prevalent pathogens and are a major burden on society. Although RNA viruses have been studied extensively, little is known about the processes that occur during the first several hours of infection because of a lack of sensitive assays. Here we develop a single-molecule imaging assay, virus infection real-time imaging (VIRIM), to study translation and replication of individual RNA viruses in live cells. VIRIM uncovered a striking heterogeneity in replication dynamics between cells and revealed extensive coordination between translation and replication of single viral RNAs. Furthermore, using VIRIM, we identify the replication step of the incoming viral RNA as a major bottleneck of successful infection and identify host genes that are responsible for inhibition of early virus replication. Single-molecule imaging of virus infection is a powerful tool to study virus replication and virus-host interactions that may be broadly applicable to RNA viruses.

RNA 3' uridylation occurs pervasively in eukaryotes, but is poorly characterized in viruses. In this study, we demonstrate that a broad array of RNA viruses, including mycoviruses, plant viruses and animal viruses, possess a novel population of RNA species bearing nontemplated oligo(U) or (U)-rich tails, suggesting widespread 3' uridylation in eukaryotic viruses. Given the biological relevance of 3' uridylation to eukaryotic RNA degradation, we propose a conserved but as-yet-unknown mechanism in virus-host interaction.

Understanding the microbiome of ticks in Australia is of considerable interest given the ongoing debate over whether Lyme disease and its causative agent, the bacterium Borrelia burgdorferisensu lato, are present in Australia. The diversity of bacteria infecting Australian ticks has been studied using both culture- and metagenomics-based techniques. However, little is known about the virome of Australian ticks, including whether this includes viruses with the potential to infect mammals. We used a meta-transcriptomics approach to reveal the diversity and evolution of viruses from Australian ticks collected from two locations on the central east coast of Australia, including metropolitan Sydney. From this we identified 19 novel RNA viruses belonging to 12 families, as well as 1 previously described RNA virus. The majority of these viruses were related to arthropod-associated viruses, suggesting that they do not utilize mammalian hosts. However, two novel viruses discovered in ticks feeding on bandicoot marsupials clustered closely within the mammal-associated hepacivirus and pestivirus groups (family Flaviviridae). Another bandicoot tick yielded a novel coltivirus (family Reoviridae), a group of largely tick-associated viruses containing the known human pathogen Colorado tick fever virus and its relative, Eyach virus. Importantly, our transcriptomic data provided no evidence for the presence of B. burgdorferisensu lato in any tick sample, providing further evidence against the presence of Lyme disease in Australia. In sum, this study reveals that Australian ticks harbor a diverse virome, including some viruses that merit additional screening in the context of emerging infectious disease.IMPORTANCE Each year a growing number of individuals along the east coast of Australia experience debilitating disease following tick bites. As there is no evidence for the presence of the causative agent of Lyme disease, Borrelia burgdorferisensu lato, in Australian ticks, the etiological basis of this disease syndrome remains controversial. To characterize the viruses associated with Australian ticks, particularly those that might be associated with mammalian infection, we performed unbiased RNA sequencing on 146 ticks collected across two locations along the coast of New South Wales, Australia. This revealed 19 novel RNA viruses from a diverse set of families. Notably, three of these viruses clustered with known mammalian viruses, including a novel coltivirus that was related to the human pathogen Colorado tick fever virus.

Infection of human cells with poliovirus induces the proliferation of double-membraned cytoplasmic vesicles whose surfaces are used as the sites of viral RNA replication and whose origin is unknown. Here, we show that several hallmarks of cellular autophagosomes can be identified in poliovirus-induced vesicles, including colocalization of LAMP1 and LC3, the human homolog of Saccharomyces cerevisiae Atg8p, and staining with the fluorophore monodansylcadaverine followed by fixation. Colocalization of LC3 and LAMP1 was observed early in the poliovirus replicative cycle, in cells infected with rhinoviruses 2 and 14, and in cells that express poliovirus proteins 2BC and 3A, known to be sufficient to induce double-membraned vesicles. Stimulation of autophagy increased poliovirus yield, and inhibition of the autophagosomal pathway by 3-methyladenine or by RNA interference against mRNAs that encode two different proteins known to be required for autophagy decreased poliovirus yield. We propose that, for poliovirus and rhinovirus, components of the cellular machinery of autophagosome formation are subverted to promote viral replication. Although autophagy can serve in the innate immune response to microorganisms, our findings are inconsistent with a role for the induced autophagosome-like structures in clearance of poliovirus. Instead, we argue that these double-membraned structures provide membranous supports for viral RNA replication complexes, possibly enabling the nonlytic release of cytoplasmic contents, including progeny virions, from infected cells.

RNA viruses are the etiological agents of many infectious diseases. Since RNA viruses are error-prone during genome replication, rapid, accurate and economical whole RNA viral genome sequence determination is highly demanded. Next-generation sequencing (NGS) techniques perform whole viral genome sequencing due to their high-throughput sequencing capacity. However, the NGS techniques involve a significant burden for sample preparation. Since to generate complete viral genome coverage, genomic nucleic acid enrichment is required by reverse transcription PCR using virus-specific primers or by viral particle concentration. Furthermore, conventional NGS techniques cannot determine the 5' and 3' terminal sequences of the RNA viral genome. Therefore, the terminal sequences are determined one by one using rapid amplification of cDNA ends (RACE). However, since some RNA viruses have segmented genomes, the burden of the determination using RACE is proportional to the number of segments. To date, there is only one study attempting whole genome sequencing of multiple RNA viruses without using above mentioned methods, but the generated sequences' accuracy compared to the reference sequences was up to 97% and did not reach 100% due to the low read depth. Hence, we established novel methods, named PCR-NGS and RCA-NGS, that were optimized for an NGS machine, MinION. These methods do not require nucleic acid amplification with virus-specific PCR primers, physical viral particle enrichment, and RACE. These methods enable whole RNA viral genome sequencing by combining the following techniques: (1) removal of unwanted DNA and RNA other than the RNA viral genome by nuclease treatment; (2) the terminal of viral genome sequence determination by barcoded linkers ligation; (3) amplification of the viral genomic cDNA using ligated linker sequences-specific PCR or an isothermal DNA amplification technique, such as rolling circle amplification (RCA). The established method was evaluated using isolated RNA viruses with single-stranded, double-stranded, positive-stranded, negative-stranded, non-segmented or multi-segmented genomes. As a result, all the viral genome sequences could be determined with 100% accuracy, and these mean read depths were greater than 2,500×, at least using either of the methods. This method should allow for easy and economical determination of accurate RNA viral genomes.

RNA-dependent RNA polymerase (RdRp) is a unique and highly conserved enzyme across all members of the RNA virus superfamilies. Besides, humans do not have a homolog of this protein. Therefore, the RdRp is an attractive target for a broadly effective therapeutic agent against RNA viruses. In this study, a formerly generated cell-penetrating human single-chain antibody variable fragment (superantibody) to a conformational epitope of hepatitis C virus (HCV) RdRp, which inhibited the polymerase activity leading to the HCV replication inhibition and the host innate immunity restoration, was tested against emerging/reemerging RNA viruses. The superantibody could inhibit the replication of the other members of the Flaviviridae (DENV serotypes 1-4, ZIKV, and JEV), Picornaviridae (genus Enterovirus: EV71, CVA16), and Coronaviridae (genus Alphacoronavirus: PEDV, and genus Betacoronavirus: SARS-CoV-2 (Wuhan wild-type and the variants of concern), in a dose-dependent manner, as demonstrated by the reduction of intracellular viral RNAs and numbers of the released infectious particles. Computerized simulation indicated that the superantibody formed contact interfaces with many residues at the back of the thumb domain (thumb II site, T2) of DENV, ZIKV, JEV, EV71, and CVA16 and fingers and thumb domains of the HCV and coronaviruses (PEDV and SARS-CoV-2). The superantibody binding may cause allosteric change in the spatial conformation of the enzyme and disrupt the catalytic activity, leading to replication inhibition. Although the speculated molecular mechanism of the superantibody needs experimental support, existing data indicate that the superantibody has high potential as a non-chemical broadly effective anti-positive sense-RNA virus agent.

Weeds often grow alongside crop plants. In addition to competing with crops for nutrients, water and space, weeds host insect vectors or act as reservoirs for viral diversity. However, little is known about viruses infecting rice weeds. In this work, we used metatranscriptomic deep sequencing to identify RNA viruses from 29 weed samples representing 23 weed species. A total of 224 RNA viruses were identified: 39 newly identified viruses are sufficiently divergent to comprise new families and genera. The newly identified RNA viruses clustered within 18 viral families. Of the identified viruses, 196 are positive-sense single-stranded RNA viruses, 24 are negative-sense single-stranded RNA viruses and 4 are double-stranded RNA viruses. We found that some novel RNA viruses clustered within the families or genera of several plant virus species and have the potential to infect plants. Collectively, these results expand our understanding of viral diversity in rice weeds. Our work will contribute to developing effective strategies with which to manage the spread and epidemiology of plant viruses.

Arctic permafrost is thawing due to global warming, with unknown consequences on the microbial inhabitants or associated viruses. DNA viruses have previously been shown to be abundant and active in thawing permafrost, but little is known about RNA viruses in these systems. To address this knowledge gap, we assessed the composition of RNA viruses in thawed permafrost samples that were incubated for 97 days at 4°C to simulate thaw conditions. A diverse RNA viral community was assembled from metatranscriptome data including double-stranded RNA viruses, dominated by Reoviridae and Hypoviridae, and negative and positive single-stranded RNA viruses, with relatively high representations of Rhabdoviridae and Leviviridae, respectively. Sequences corresponding to potential plant and human pathogens were also detected. The detected RNA viruses primarily targeted dominant eukaryotic taxa in the samples (e.g., fungi, Metazoa and Viridiplantae) and the viral community structures were significantly associated with predicted host populations. These results indicate that RNA viruses are linked to eukaryotic host dynamics. Several of the RNA viral sequences contained auxiliary metabolic genes encoding proteins involved in carbon utilization (e.g., polygalacturosase), implying their potential roles in carbon cycling in thawed permafrost. IMPORTANCE Permafrost is thawing at a rapid pace in the Arctic with largely unknown consequences on ecological processes that are fundamental to Arctic ecosystems. This is the first study to determine the composition of RNA viruses in thawed permafrost. Other recent studies have characterized DNA viruses in thawing permafrost, but the majority of DNA viruses are bacteriophages that target bacterial hosts. By contrast RNA viruses primarily target eukaryotic hosts and thus represent potential pathogenic threats to humans, animals, and plants. Here, we find that RNA viruses in permafrost are novel and distinct from those in other habitats studied to date. The COVID-19 pandemic has heightened awareness of the importance of potential environmental reservoirs of emerging RNA viral pathogens. We demonstrate that some potential pathogens were detected after an experimental thawing regime. These results are important for understanding critical viral-host interactions and provide a better understanding of the ecological roles that RNA viruses play as permafrost thaws.

Acyclovir (ACV) is an effective antiviral agent for treating lytic Herpes Simplex virus, type 1 (HSV-1) infections, and it has dramatically reduced the mortality rate of herpes simplex encephalitis. However, HSV-1 resistance to ACV and its derivatives is being increasingly documented, particularly among immunocompromised individuals. The burgeoning drug resistance compels the search for a new generation of more efficacious anti-herpetic drugs. We have previously shown that trans-dihydrolycoricidine (R430), a lycorane-type alkaloid derivative, effectively inhibits HSV-1 infections in cultured cells. We now report that R430 also inhibits ACV-resistant HSV-1 strains, accompanied by global inhibition of viral gene transcription and enrichment of H3K27me3 methylation on viral gene promoters. Furthermore, we demonstrate that R430 prevents HSV-1 reactivation from latency in an ex vivo rodent model. Finally, among a panel of DNA viruses and RNA viruses, R430 inhibited Zika virus with high therapeutic index. Its therapeutic index is comparable to standard antiviral drugs, though it has greater toxicity in non-neuronal cells than in neuronal cells. Synthesis of additional derivatives could enable more efficacious antivirals and the identification of active pharmacophores.