An exploration of ambigrammatic sequences in narnaviruses

Viromes of Tabanids from Russia

Advances in sequencing technologies and bioinformatics have greatly enhanced our knowledge of virus biodiversity. Currently, the viromes of hematophagous invertebrates, such as mosquitoes and ixodid ticks, are being actively studied. Tabanidae (Diptera) are a widespread family, with members mostly known for their persistent hematophagous behavior. They transmit viral, bacterial, and other pathogens, both biologically and mechanically. However, tabanid viromes remain severely understudied. In this study, we used high-throughput sequencing to describe the viromes of several species in the Hybomitra, Tabanus, Chrysops, and Haematopota genera, which were collected in two distant parts of Russia: the Primorye Territory and Ryazan Region. We assembled fourteen full coding genomes of novel viruses, four partial coding genomes, as well as several fragmented viral sequences, which presumably belong to another twelve new viruses. All the discovered viruses were tested for their ability to replicate in mammalian porcine embryo kidney (PEK), tick HAE/CTVM8, and mosquito C6/36 cell lines. In total, 16 viruses were detected in at least one cell culture after three passages (for PEK and C6/36) or 3 weeks of persistence in HAE/CTVM8. However, in the majority of cases, qPCR showed a decline in virus load over time.

The virome of the invasive Asian bush mosquito Aedes japonicus in Europe

The Asian bush mosquito Aedes japonicus is rapidly invading North America and Europe. Due to its potential to transmit multiple pathogenic arthropod-borne (arbo)viruses including Zika virus, West Nile virus, and chikungunya virus, it is important to understand the biology of this vector mosquito in more detail. In addition to arboviruses, mosquitoes can also carry insect-specific viruses that are receiving increasing attention due to their potential effects on host physiology and arbovirus transmission. In this study, we characterized the collection of viruses, referred to as the virome, circulating in Ae. japonicus populations in the Netherlands and France. Applying a small RNA-based metagenomic approach to Ae. japonicus, we uncovered a distinct group of viruses present in samples from both the Netherlands and France. These included one known virus, Ae. japonicus narnavirus 1 (AejapNV1), and three new virus species that we named Ae. japonicus totivirus 1 (AejapTV1), Ae. japonicus anphevirus 1 (AejapAV1) and Ae. japonicus bunyavirus 1 (AejapBV1). We also discovered sequences that were presumably derived from two additional novel viruses: Ae. japonicus bunyavirus 2 (AejapBV2) and Ae. japonicus rhabdovirus 1 (AejapRV1). All six viruses induced strong RNA interference responses, including the production of twenty-one nucleotide-sized small interfering RNAs, a signature of active replication in the host. Notably, AejapBV1 and AejapBV2 belong to different viral families; however, no RNA-dependent RNA polymerase sequence has been found for AejapBV2. Intriguingly, our small RNA-based approach identified an ∼1-kb long ambigrammatic RNA that is associated with AejapNV1 as a secondary segment but showed no similarity to any sequence in public databases. We confirmed the presence of AejapNV1 primary and secondary segments, AejapTV1, AejapAV1, and AejapBV1 by reverse transcriptase polymerase chain reaction (PCR) in wild-caught Ae. japonicus mosquitoes. AejapNV1 and AejapTV1 were found at high prevalence (87-100 per cent) in adult females, adult males, and larvae. Using a small RNA-based, sequence-independent metagenomic strategy, we uncovered a conserved and prevalent virome among Ae. japonicus mosquito populations. The high prevalence of AejapNV1 and AejapTV1 across all tested mosquito life stages suggests that these viruses are intimately associated with Ae. japonicus.

A novel ambigrammatic mycovirus, PsV5, works hand in glove with wheat stripe rust fungus to facilitate infection

Here we describe a novel narnavirus, Puccinia striiformis virus 5 (PsV5), from the devastating wheat stripe rust fungus P. striiformis f. sp. tritici (Pst). The genome of PsV5 contains two predicted open reading frames (ORFs) that largely overlap on reverse strands: an RNA-dependent RNA polymerase (RdRp) and a reverse-frame ORF (rORF) with unknown function. Protein translations of both ORFs were demonstrated by immune technology. Transgenic wheat lines overexpressing PsV5 (RdRp-rORF), RdRp ORF, or rORF were more susceptible to Pst infection, whereas PsV5-RNA interference (RNAi) lines were more resistant. Overexpression of PsV5 (RdRp-rORF), RdRp ORF, or rORF in Fusarium graminearum also boosted fungal virulence. We thus report a novel ambigrammatic mycovirus that promotes the virulence of its fungal host. The results are a significant addition to our understanding of virosphere diversity and offer insights for sustainable wheat rust disease control.

Putative Mitoviruses without In-Frame UGA(W) Codons: Evolutionary Implications

Mitoviruses are small vertically transmitted RNA viruses found in fungi, plants and animals. Taxonomically, a total of 105 species and 4 genera have been formally recognized by ICTV, and recently, 18 new putative species have been included in a new proposed genus. Transcriptomic and metatranscriptomic studies are a major source of countless new virus-like sequences that are continually being added to open databases and these may be good sources for identifying new putative mitoviruses. The search for mitovirus-like sequences in the NCBI databases resulted in the discovery of more than one hundred new putative mitoviruses, with important implications for taxonomy and also for the evolutionary scenario. Here, we propose the inclusion of four new putative members to the genus Kvaramitovirus, and the existence of a new large basally divergent lineage composed of 144 members that lack internal UGA codons (subfamily "Arkeomitovirinae"), a feature not shared by the vast majority of mitoviruses. Finally, a taxonomic categorization proposal and a detailed description of the evolutionary history of mitoviruses were carried out. This in silico study supports the hypothesis of the existence of a basally divergent lineage that could have had an impact on the early evolutionary history of mitoviruses.

Mycoviruses as a part of the global virome: Diversity, evolutionary links and lifestyle

Knowledge of mycovirus diversity, evolution, horizontal gene transfer and shared ancestry with viruses infecting distantly related hosts, such as plants and arthropods, has increased vastly during the last few years due to advances in the high throughput sequencing methodologies. This also has enabled the discovery of novel mycoviruses with previously unknown genome types, mainly new positive and negative single-stranded RNA mycoviruses ((+) ssRNA and (-) ssRNA) and single-stranded DNA mycoviruses (ssDNA), and has increased our knowledge of double-stranded RNA mycoviruses (dsRNA), which in the past were thought to be the most common viruses infecting fungi. Fungi and oomycetes (Stramenopila) share similar lifestyles and also have similar viromes. Hypothesis about the origin and cross-kingdom transmission events of viruses have been raised and are supported by phylogenetic analysis and by the discovery of natural exchange of viruses between different hosts during virus-fungus coinfection in planta. In this review we make a compilation of the current information on the genome organization, diversity and taxonomy of mycoviruses, discussing their possible origins. Our focus is in recent findings suggesting the expansion of the host range of many viral taxa previously considered to be exclusively fungal, but we also address factors affecting virus transmissibility and coexistence in single fungal or oomycete isolates, as well as the development of synthetic mycoviruses and their use in investigating mycovirus replication cycles and pathogenicity.

Identification of Three Viruses Infecting Mulberry Varieties

Viruses-mediated genome editing in plants is a powerful strategy to develop plant cultivars with important and novel agricultural traits. Mulberry alba is an important economic tree species that has been cultivated in China for more than 5000 years. So far, only a few viruses have been identified from mulberry trees, and their application potential is largely unknown. Therefore, mining more virus resources from the mulberry tree can pave the way for the establishment of useful engineering tools. In this study, eight old mulberry plants were gathered in seven geographic areas for virome analysis. Based on transcriptome analysis, we discovered three viruses associated with mulberries: Citrus leaf blotch virus isolate mulberry alba 2 (CLBV-ML2), Mulberry-associated virga-like virus (MaVLV), and Mulberry-associated narna-like virus (MaNLV). The genome of CLBV-ML2 was completely sequenced and exhibited high homology with Citriviruses, considered to be members of the genus Citrivirus, while the genomes of MaVLV and MaNLV were nearly completed lacking the 5' and 3' termini sequences. We tentatively consider MaVLV to be members of the family Virgaviridae and MaNLV to be members of the genus Narnavirus based on the results of phylogenetic trees. The infection experiments showed that CLBV-ML2 could be detected in the inoculated seedlings of both N. benthamiana and Morus alba, while MaVLV could only be detected in N. benthamiana. All of the infected seedlings did not show obvious symptoms.

Discovery and Evolution of Six Positive-Sense RNA Viruses Co-infecting the Hypovirulent Strain SCH733 of Sclerotinia sclerotiorum

Sclerotinia sclerotiorum is a well-known phytopathogenic fungus with a wide host range. Identifying novel mycoviruses in phytopathogenic fungi is necessary to develop novel strategies for plant health protection and contribute to understanding the origin of viruses. Six new mycoviruses with positive single-stranded RNA genomes co-infecting the hypovirulent strain SCH733 of S. sclerotiorum were identified using a metatranscriptomic approach, and their complete genome sequences were molecularly determined. These mycoviruses belong to the following five families: Narnaviridae, Mitoviridae, Deltaflexviridae, Botourmiaviridae, and Ambiguiviridae. Three of these mycoviruses belong to existing International Committee on Taxonomy of Viruses (ICTV)-recognized species. Two of these newly identified mycoviruses have unique genomic features that are significantly different from those of all known mycoviruses. Phylogenetic analysis revealed that these six mycoviruses included close as well as distant relatives of known mycoviruses, thereby providing new insight into virus evolution and classification. Mycovirus horizontal transmission and elimination experiments revealed that Sclerotinia sclerotiorum narnavirus 5 is associated with hypovirulence of S. sclerotiorum, although we have not shown that it is independently responsible for the hypovirulence phenotype. This study broadens the diversity of known mycoviruses infecting S. sclerotiorum and provides a clue toward limiting hypovirulence in S. sclerotiorum.

Mycovirus Diversity and Evolution Revealed/Inferred from Recent Studies

High-throughput virome analyses with various fungi, from cultured or uncultured sources, have led to the discovery of diverse viruses with unique genome structures and even neo-lifestyles. Examples in the former category include splipalmiviruses and ambiviruses. Splipalmiviruses, related to yeast narnaviruses, have multiple positive-sense (+) single-stranded (ss) RNA genomic segments that separately encode the RNA-dependent RNA polymerase motifs, the hallmark of RNA viruses (members of the kingdom Orthornavirae). Ambiviruses appear to have an undivided ssRNA genome of 3∼5 kb with two large open reading frames (ORFs) separated by intergenic regions. Another narna-like virus group has two fully overlapping ORFs on both strands of a genomic segment that span more than 90% of the genome size. New virus lifestyles exhibited by mycoviruses include the yado-kari/yado-nushi nature characterized by the partnership between the (+)ssRNA yadokarivirus and an unrelated dsRNA virus (donor of the capsid for the former) and the hadaka nature of capsidless 10-11 segmented (+)ssRNA accessible by RNase in infected mycelial homogenates. Furthermore, dsRNA polymycoviruses with phylogenetic affinity to (+)ssRNA animal caliciviruses have been shown to be infectious as dsRNA-protein complexes or deproteinized naked dsRNA. Many previous phylogenetic gaps have been filled by recently discovered fungal and other viruses, which haveprovided interesting evolutionary insights. Phylogenetic analyses and the discovery of natural and experimental cross-kingdom infections suggest that horizontal virus transfer may have occurred and continue to occur between fungi and other kingdoms.

High Diversity of Novel Viruses in the Tree Pathogen Phytophthora castaneae Revealed by High-Throughput Sequencing of Total and Small RNA

Phytophthora castaneae, an oomycete pathogen causing root and trunk rot of different tree species in Asia, was shown to harbor a rich diversity of novel viruses from different families. Four P. castaneae isolates collected from Chamaecyparis hodginsii in a semi-natural montane forest site in Vietnam were investigated for viral presence by traditional and next-generation sequencing (NGS) techniques, i.e., double-stranded RNA (dsRNA) extraction and high-throughput sequencing (HTS) of small RNAs (sRNAs) and total RNA. Genome organization, sequence similarity, and phylogenetic analyses indicated that the viruses were related to members of the order Bunyavirales and families Endornaviridae, Megabirnaviridae, Narnaviridae, Totiviridae, and the proposed family “Fusagraviridae.” The study describes six novel viruses: Phytophthora castaneae RNA virus 1–5 (PcaRV1-5) and Phytophthora castaneae negative-stranded RNA virus 1 (PcaNSRV1). All six viruses were detected by sRNA sequencing, which demonstrates an active RNA interference (RNAi) system targeting viruses in P. castaneae. To our knowledge, this is the first report of viruses in P. castaneae and the whole Phytophthora major Clade 5, as well as of the activity of an RNAi mechanism targeting viral genomes among Clade 5 species. PcaRV1 is the first megabirnavirus described in oomycetes and the genus Phytophthora.

Resolving Deep Evolutionary Relationships within the RNA Virus Phylum Lenarviricota

The RNA virus phylum Lenarviricota is composed of the fungi-associated families Narnaviridae and Mitoviridae, the RNA bacteriophage Leviviridae, and the plant and fungi-associated Botourmiaviridae. Members of the Lenarviricota are abundant in most environments and boast remarkable phylogenetic and genomic diversity. As this phylum includes both RNA bacteriophage and fungi- and plant-associated species, the Lenarviricota likely mark a major evolutionary transition between those RNA viruses associated with prokaryotes and eukaryotes. Despite the remarkable expansion of this phylum following metagenomic studies, the phylogenetic relationships among the families within the Lenarviricota remain uncertain. Utilising a large data set of relevant viral sequences, we performed phylogenetic and genomic analyses to resolve the complex evolutionary history within this phylum and identify patterns in the evolution of virus genome organisation. Despite limitations reflecting very high levels of sequence diversity, our phylogenetic analyses suggest that the Leviviridae comprise the basal lineage within the Lenarviricota. Our phylogenetic results also support the construction of a new virus family—the Narliviridae—comprising a set of diverse and phylogenetically distinct species, including a number of uniquely encapsidated viruses. We propose a taxonomic restructuring within the Lenarviricota to better reflect the phylogenetic relationships documented here, with the Botourmiaviridae and Narliviridae combined into the order Ourlivirales, the Narnaviridae remaining in the order Wolframvirales, and these orders combined into the single class, the Amabiliviricetes. In sum, this study provides insights into the complex evolutionary relationships among the diverse families that make up the Lenarviricota.

Spotlight on alternative frame coding: Two long overlapping genes in Pseudomonas aeruginosa are translated and under purifying selection

The existence of overlapping genes (OLGs) with significant coding overlaps revolutionizes our understanding of genomic complexity. We report two exceptionally long (957 nt and 1536 nt), evolutionarily novel, translated antisense open reading frames (ORFs) embedded within annotated genes in the pathogenic Gram-negative bacterium Pseudomonas aeruginosa. Both OLG pairs show sequence features consistent with being genes and transcriptional signals in RNA sequencing. Translation of both OLGs was confirmed by ribosome profiling and mass spectrometry. Quantitative proteomics of samples taken during different phases of growth revealed regulation of protein abundances, implying biological functionality. Both OLGs are taxonomically restricted, and likely arose by overprinting within the genus. Evidence for purifying selection further supports functionality. The OLGs reported here, designated olg1 and olg2, are the longest yet proposed in prokaryotes and are among the best attested in terms of translation and evolutionary constraint. These results highlight a potentially large unexplored dimension of prokaryotic genomes.

•

Two novel, very long, overlapping genes were found in Pseudomonas aeruginosa

•

Both overlapping genes, olg1 and olg2, are transcribed, translated, and regulated

•

Mass spectrometry verifies translation of the overlapping and their mother genes

•

Both overlapping genes are taxonomically restricted, but under purifying selection

Revealing RNA virus diversity and evolution in unicellular algae transcriptomes

Remarkably little is known about the diversity and evolution of RNA viruses in unicellular eukaryotes. We screened a total of 570 transcriptomes from the Marine Microbial Eukaryote Transcriptome Sequencing Project that encompasses a wide diversity of microbial eukaryotes, including most major photosynthetic lineages (i.e. the microalgae). From this, we identified thirty new and divergent RNA virus species, occupying a range of phylogenetic positions within the overall diversity of RNA viruses. Approximately one-third of the newly described viruses comprised single-stranded positive-sense RNA viruses from the order Lenarviricota associated with fungi, plants, and protists, while another third were related to the order Ghabrivirales, including members of the protist and fungi-associated Totiviridae. Other viral species showed sequence similarity to positive-sense RNA viruses from the algae-associated Marnaviridae, the double-stranded RNA (ds-RNA) Partitiviridae, as well as tentative evidence for one negative-sense RNA virus related to the Qinviridae. Importantly, we were able to identify divergent RNA viruses from distant host taxa, revealing the ancestry of these viral families and greatly extending our knowledge of the RNA viromes of microalgal cultures. Both the limited number of viruses detected per sample and the low sequence identity to known RNA viruses imply that additional microalgal viruses exist that could not be detected at the current sequencing depth or were too divergent to be identified using sequence similarity. Together, these results highlight the need for further investigation of algal-associated RNA viruses as well as the development of new tools to identify RNA viruses that exhibit very high levels of sequence divergence.

Global genome analysis reveals a vast and dynamic anellovirus landscape within the human virome

Anelloviruses are a ubiquitous component of healthy human viromes and remain highly prevalent after being acquired early in life. The full extent of "anellome" diversity and its evolutionary dynamics remain unexplored. We employed in-depth sequencing of blood-transfusion donor(s)-recipient pairs coupled with public genomic resources for a large-scale assembly of anellovirus genomes and used the data to characterize global and personal anellovirus diversity through time. The breadth of the anellome is much greater than previously appreciated, and individuals harbor unique anellomes and transmit lineages that can persist for several months within a diverse milieu of endemic host lineages. Anellovirus sequence diversity is shaped by extensive recombination at all levels of divergence, hindering traditional phylogenetic analyses. Our findings illuminate the transmission dynamics and vast diversity of anelloviruses and set the foundation for future studies to characterize their biology.

Polysomally protected viruses

It is conceivable that an RNA virus could use a polysome, that is, a string of ribosomes covering the RNA strand, to protect the genetic material from degradation inside a host cell. This paper discusses how such a virus might operate, and how its presence might be detected by ribosome profiling. There are two possible forms for such apolysomally protected virus, depending upon whether just the forward strand or both the forward and complementary strands can be encased by ribosomes (these will be termed type 1 and type 2, respectively). It is argued that in the type 2 case the viral RNA would evolve anambigrammaticproperty, whereby the viral genes are free of stop codons in a reverse reading frame (with forward and reverse codons aligned). Recent observations of ribosome profiles of ambigrammatic narnavirus sequences are consistent with our predictions for the type 2 case.

Persistence of Ambigrammatic Narnaviruses Requires Translation of the Reverse Open Reading Frame

Narnaviruses are RNA viruses detected in diverse fungi, plants, protists, arthropods, and nematodes. Though initially described as simple single-gene nonsegmented viruses encoding RNA-dependent RNA polymerase (RdRp), a subset of narnaviruses referred to as "ambigrammatic" harbor a unique genomic configuration consisting of overlapping open reading frames (ORFs) encoded on opposite strands. Phylogenetic analysis supports selection to maintain this unusual genome organization, but functional investigations are lacking. Here, we establish the mosquito-infecting Culex narnavirus 1 (CxNV1) as a model to investigate the functional role of overlapping ORFs in narnavirus replication. In CxNV1, a reverse ORF without homology to known proteins covers nearly the entire 3.2-kb segment encoding the RdRp. Additionally, two opposing and nearly completely overlapping novel ORFs are found on the second putative CxNV1 segment, the 0.8-kb "Robin" RNA. We developed a system to launch CxNV1 in a naive mosquito cell line and then showed that functional RdRp is required for persistence of both segments, and an intact reverse ORF is required on the RdRp segment for persistence. Mass spectrometry of persistently CxNV1-infected cells provided evidence for translation of this reverse ORF. Finally, ribosome profiling yielded a striking pattern of footprints for all four CxNV1 RNA strands that was distinct from actively translating ribosomes on host mRNA or coinfecting RNA viruses. Taken together, these data raise the possibility that the process of translation itself is important for persistence of ambigrammatic narnaviruses, potentially by protecting viral RNA with ribosomes, thus suggesting a heretofore undescribed viral tactic for replication and transmission. IMPORTANCE Fundamental to our understanding of RNA viruses is a description of which strand(s) of RNA are transmitted as the viral genome relative to which encode the viral proteins. Ambigrammatic narnaviruses break the mold. These viruses, found broadly in fungi, plants, and insects, have the unique feature of two overlapping genes encoded on opposite strands, comprising nearly the full length of the viral genome. Such extensive overlap is not seen in other RNA viruses and comes at the cost of reduced evolutionary flexibility in the sequence. The present study is motivated by investigating the benefits which balance that cost. We show for the first time a functional requirement for the ambigrammatic genome configuration in Culex narnavirus 1, which suggests a model for how translation of both strands might benefit this virus. Our work highlights a new blueprint for viral persistence, distinct from strategies defined by canonical definitions of the coding strand.

Novel Mycoviruses Discovered in the Mycovirome of a Necrotrophic Fungus

Botrytis cinerea is one of the most important plant-pathogenic fungus. Products based on microorganisms can be used in biocontrol strategies alternative to chemical control, and mycoviruses have been explored as putative biological agents in such approaches. Here, we have explored the mycovirome of B. cinerea isolates from grapevine of Italy and Spain to increase the knowledge about mycoviral diversity and evolution, and to search for new widely distributed mycoviruses that could be active ingredients in biological products to control this hazardous fungus. A total of 248 B. cinerea field isolates were used for our metatranscriptomic study. Ninety-two mycoviruses were identified: 62 new mycoviral species constituting putative novel viral genera and families. Of these mycoviruses, 57 had a positive-sense single-stranded RNA (ssRNA) genome, 19 contained a double-stranded RNA (dsRNA) genome, 15 had a negative-sense ssRNA genome, and 1 contained a single-stranded DNA (ssDNA) genome. In general, ssRNA mycoviruses were widely distributed in all sampled regions, the ssDNA mycovirus was more frequently found in Spain, and dsRNA mycoviruses were scattered in some pools of both countries. Some of the identified mycoviruses belong to clades that have never been found associated with Botrytis species: Botrytis-infecting narnaviruses; alpha-like, umbra-like, and tymo-like ssRNA+ mycoviruses; trisegmented ssRNA- mycovirus; bisegmented and tetrasegmented dsRNA mycoviruses; and finally, an ssDNA mycovirus. Among the results obtained in this massive mycovirus screening, the discovery of novel bisegmented viruses, phylogenetically related to narnaviruses, is remarkable.IMPORTANCE The results obtained here have expanded our knowledge of mycoviral diversity, horizontal transfers, and putative cross-kingdom events. To date, this study presents the most extensive and wide diversity collection of mycoviruses infecting the necrotrophic fungus B. cinerea The collection included all types of mycoviruses, with dsRNA, ssRNA+, ssRNA-, and ssDNA genomes, most of which were discovered here, and some of which were previously reported as infecting B. cinerea or other plant-pathogenic fungi. Some of these mycoviruses are reported for the first time here associated with B. cinerea, as a trisegmented ssRNA- mycovirus and as an ssDNA mycovirus, but even more remarkablly, we also describe here four novel bisegmented viruses (binarnaviruses) not previously described in nature. The present findings significantly contribute to general knowledge in virology and more particularly in the field of mycovirology.

Simultaneous codon usage, the origin of the proteome, and the emergence of de-novo proteins

Genetic coding generally uses only one of a gene's two strands; its complement serving as template for replication. Aminoacyl-tRNA synthetases, aaRS, apparently first emerged as pairs on bidirectional genes, in which anticodons in the template strand served as codons for an entirely different protein. Interpreting both strands in frame constrained such genes sufficiently that it was rapidly superseded, leaving only traces in the elevated pairing between codon middle bases in antiparallel alignments. Codon assignments actually promote using information from both strands in multiple reading frames. Related phenomena, known as overprinting, are widely associated with viruses. In-frame bidirectional coding and overprinting nevertheless imply different structural and functional relationships, and different roles in generating folded proteins throughout the evolution of the proteome.

Polymorphism of genetic ambigrams

Double synonyms in the genetic code can be used as a tool to test competing hypotheses regarding ambigrammatic narnavirus genomes. Applying the analysis to recent observations of Culex narnavirus 1 and Zhejiang mosquito virus 3 ambigrammatic viruses indicates that the open reading frame on the complementary strand of the segment coding for RNA-dependent RNA polymerase does not code for a functional protein. Culex narnavirus 1 has been shown to possess a second segment, also ambigrammatic, termed 'Robin'. We find a comparable segment for Zhejiang mosquito virus 3, a moderately diverged relative of Culex narnavirus 1. Our analysis of Robin polymorphisms suggests that its reverse open reading frame also does not code for a functional protein. We make a hypothesis about its role.

Characterization of a Novel Mitovirus of the Sand Fly Lutzomyia longipalpis Using Genomic and Virus-Host Interaction Signatures

Hematophagous insects act as the major reservoirs of infectious agents due to their intimate contact with a large variety of vertebrate hosts. Lutzomyia longipalpis is the main vector of Leishmania chagasi in the New World, but its role as a host of viruses is poorly understood. In this work, Lu. longipalpis RNA libraries were subjected to progressive assembly using viral profile HMMs as seeds. A sequence phylogenetically related to fungal viruses of the genus Mitovirus was identified and this novel virus was named Lul-MV-1. The 2697-base genome presents a single gene coding for an RNA-directed RNA polymerase with an organellar genetic code. To determine the possible host of Lul-MV-1, we analyzed the molecular characteristics of the viral genome. Dinucleotide composition and codon usage showed profiles similar to mitochondrial DNA of invertebrate hosts. Also, the virus-derived small RNA profile was consistent with the activation of the siRNA pathway, with size distribution and 5′ base enrichment analogous to those observed in viruses of sand flies, reinforcing Lu. longipalpis as a putative host. Finally, RT-PCR of different insect pools and sequences of public Lu. longipalpis RNA libraries confirmed the high prevalence of Lul-MV-1. This is the first report of a mitovirus infecting an insect host.

Soybean Thrips (Thysanoptera: Thripidae) Harbor Highly Diverse Populations of Arthropod, Fungal and Plant Viruses

Soybean thrips (Neohydatothrips variabilis) are one of the most efficient vectors of soybean vein necrosis virus, which can cause severe necrotic symptoms in sensitive soybean plants. To determine which other viruses are associated with soybean thrips, the metatranscriptome of soybean thrips, collected by the Midwest Suction Trap Network during 2018, was analyzed. Contigs assembled from the data revealed a remarkable diversity of virus-like sequences. Of the 181 virus-like sequences identified, 155 were novel and associated primarily with taxa of arthropod-infecting viruses, but sequences similar to plant and fungus-infecting viruses were also identified. The novel viruses were predicted to have positive-sense RNA, negative-stranded RNA, double-stranded RNA, and single-stranded DNA genomes. The assembled sequences included 100 contigs that represented at least 95% coverage of a virus genome or genome segment. Sequences represented 12 previously described arthropod viruses including eight viruses reported from Hubei Province in China, and 12 plant virus sequences of which six have been previously described. The presence of diverse populations of plant viruses within soybean thrips suggests they feed on and acquire viruses from multiple host plant species that could be transmitted to soybean. Assessment of the virome of soybean thrips provides, for the first time, information on the diversity of viruses present in thrips.

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.

The virome from a collection of endomycorrhizal fungi reveals new viral taxa with unprecedented genome organization

Mutualistic plant-associated fungi are recognized as important drivers in plant evolution, diversity, and health. The discovery that mycoviruses can take part and play important roles in symbiotic tripartite interactions has prompted us to study the viromes associated with a collection of ericoid and orchid mycorrhizal (ERM and ORM, respectively) fungi. Our study, based on high-throughput sequencing of transcriptomes (RNAseq) from fungal isolates grown in axenic cultures, revealed in both ERM and ORM fungi the presence of new mycoviruses closely related to already classified virus taxa, but also new viruses that expand the boundaries of characterized RNA virus diversity to previously undescribed evolutionary trajectories. In ERM fungi, we provide first evidence of a bipartite virus, distantly related to narnaviruses, that splits the RNA-dependent RNA polymerase (RdRP) palm domain into two distinct proteins, encoded by each of the two segments. Furthermore, in one isolate of the ORM fungus Tulasnella spp. we detected a 12 kb genomic fragment coding for an RdRP with features of bunyavirus-like RdRPs. However, this 12 kb genomic RNA has the unique features, for Bunyavirales members, of being tri-cistronic and carrying ORFs for the putative RdRP and putative nucleocapsid in ambisense orientation on the same genomic RNA. Finally, a number of ORM fungal isolates harbored a group of ambisense bicistronic viruses with a genomic size of around 5 kb, where we could identify a putative RdRP palm domain that has some features of plus strand RNA viruses; these new viruses may represent a new lineage in the Riboviria, as they could not be reliably assigned to any of the branches in the recently derived monophyletic tree that includes most viruses with an RNA genome.

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

Despite its overwhelming clinical importance, the SARS-CoV-2 gene set remains unresolved, hindering dissection of COVID-19 biology. Here, we use comparative genomics to provide a high-confidence protein-coding gene set, characterize protein-level and nucleotide-level evolutionary constraint, and prioritize functional mutations from the ongoing COVID-19 pandemic. We select 44 complete Sarbecovirus genomes at evolutionary distances ideally-suited for protein-coding and non-coding element identification, create whole-genome alignments, and quantify protein-coding evolutionary signatures and overlapping constraint. We find strong protein-coding signatures for all named genes and for 3a, 6, 7a, 7b, 8, 9b, and also ORF3c, a novel alternate-frame gene. By contrast, ORF10, and overlapping-ORFs 9c, 3b, and 3d lack protein-coding signatures or convincing experimental evidence and are not protein-coding. Furthermore, we show no other protein-coding genes remain to be discovered. Cross-strain and within-strain evolutionary pressures largely agree at the gene, amino-acid, and nucleotide levels, with some notable exceptions, including fewer-than-expected mutations in nsp3 and Spike subunit S1, and more-than-expected mutations in Nucleocapsid. The latter also shows a cluster of amino-acid-changing variants in otherwise-conserved residues in a predicted B-cell epitope, which may indicate positive selection for immune avoidance. Several Spike-protein mutations, including D614G, which has been associated with increased transmission, disrupt otherwise-perfectly-conserved amino acids, and could be novel adaptations to human hosts. The resulting high-confidence gene set and evolutionary-history annotations provide valuable resources and insights on COVID-19 biology, mutations, and evolution.

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

Despite its overwhelming clinical importance, the SARS-CoV-2 gene set remains unresolved, hindering dissection of COVID-19 biology. Here, we use comparative genomics to provide a high-confidence protein-coding gene set, characterize protein-level and nucleotide-level evolutionary constraint, and prioritize functional mutations from the ongoing COVID-19 pandemic. We select 44 complete Sarbecovirus genomes at evolutionary distances ideally-suited for protein-coding and non-coding element identification, create whole-genome alignments, and quantify protein-coding evolutionary signatures and overlapping constraint. We find strong protein-coding signatures for all named genes and for 3a, 6, 7a, 7b, 8, 9b, and also ORF3c, a novel alternate-frame gene. By contrast, ORF10, and overlapping-ORFs 9c, 3b, and 3d lack protein-coding signatures or convincing experimental evidence and are not protein-coding. Furthermore, we show no other protein-coding genes remain to be discovered. Cross-strain and within-strain evolutionary pressures largely agree at the gene, amino-acid, and nucleotide levels, with some notable exceptions, including fewer-than-expected mutations in nsp3 and Spike subunit S1, and more-than-expected mutations in Nucleocapsid. The latter also shows a cluster of amino-acid-changing variants in otherwise-conserved residues in a predicted B-cell epitope, which may indicate positive selection for immune avoidance. Several Spike-protein mutations, including D614G, which has been associated with increased transmission, disrupt otherwise-perfectly-conserved amino acids, and could be novel adaptations to human hosts. The resulting high-confidence gene set and evolutionary-history annotations provide valuable resources and insights on COVID-19 biology, mutations, and evolution.

Virome Analysis of Aphid Populations That Infest the Barley Field: The Discovery of Two Novel Groups of Nege/Kita-Like Viruses and Other Novel RNA Viruses

Aphids (order Hemiptera) are important insect pests of crops and are also vectors of many plant viruses. However, little is known about aphid-infecting viruses, particularly their diversity and relationship to plant viruses. To investigate the aphid viromes, we performed deep sequencing analyses of the aphid transcriptomes from infested barley plants in a field in Japan. We discovered virus-like sequences related to nege/kita-, flavi-, tombus-, phenui-, mononega-, narna-, chryso-, partiti-, and luteoviruses. Using RT-PCR and sequence analyses, we determined almost complete sequences of seven nege/kitavirus-like virus genomes; one of which was a variant of the Wuhan house centipede virus (WHCV-1). The other six seem to belong to four novel viruses distantly related to Wuhan insect virus 9 (WhIV-9) or Hubei nege-like virus 4 (HVLV-4). We designated the four viruses as barley aphid RNA virus 1 to 4 (BARV-1 to -4). Moreover, some nege/kitavirus-like sequences were found by searches on the transcriptome shotgun assembly (TSA) libraries of arthropods and plants. Phylogenetic analyses showed that BARV-1 forms a clade with WHCV-1 and HVLV-4, whereas BARV-2 to -4 clustered with WhIV-9 and an aphid virus, Aphis glycines virus 3. Both virus groups (tentatively designated as Centivirus and Aphiglyvirus, respectively), together with arthropod virus-like TSAs, fill the phylogenetic gaps between the negeviruses and kitaviruses lineages. We also characterized the flavi/jingmen-like and tombus-like virus sequences as well as other RNA viruses, including six putative novel viruses, designated as barley aphid RNA viruses 5 to 10. Interestingly, we also discovered that some aphid-associated viruses, including nege/kita-like viruses, were present in different aphid species, raising a speculation that these viruses might be distributed across different aphid species with plants being the reservoirs. This study provides novel information on the diversity and spread of nege/kitavirus-related viruses and other RNA viruses that are associated with aphids.

A case for a negative-strand coding sequence in a group of positive-sense RNA viruses

Positive-sense single-stranded RNA viruses form the largest and most diverse group of eukaryote-infecting viruses. Their genomes comprise one or more segments of coding-sense RNA that function directly as messenger RNAs upon release into the cytoplasm of infected cells. Positive-sense RNA viruses are generally accepted to encode proteins solely on the positive strand. However, we previously identified a surprisingly long (∼1,000-codon) open reading frame (ORF) on the negative strand of some members of the family Narnaviridae which, together with RNA bacteriophages of the family Leviviridae, form a sister group to all other positive-sense RNA viruses. Here, we completed the genomes of three mosquito-associated narnaviruses, all of which have the long reverse-frame ORF. We systematically identified narnaviral sequences in public data sets from a wide range of sources, including arthropod, fungal, and plant transcriptomic data sets. Long reverse-frame ORFs are widespread in one clade of narnaviruses, where they frequently occupy >95 per cent of the genome. The reverse-frame ORFs correspond to a specific avoidance of CUA, UUA, and UCA codons (i.e. stop codon reverse complements) in the forward-frame RNA-dependent RNA polymerase ORF. However, absence of these codons cannot be explained by other factors such as inability to decode these codons or GC3 bias. Together with other analyses, we provide the strongest evidence yet of coding capacity on the negative strand of a positive-sense RNA virus. As these ORFs comprise some of the longest known overlapping genes, their study may be of broad relevance to understanding overlapping gene evolution and de novo origin of genes.