Discovery and evolution of bunyavirids in arctic phantom midges and ancient bunyavirid-like sequences in insect genomes

Promotion of order Bunyavirales to class Bunyaviricetes to accommodate a rapidly increasing number of related polyploviricotine viruses

Prior to 2017, the family Bunyaviridae included five genera of arthropod and rodent viruses with tri-segmented negative-sense RNA genomes related to the Bunyamwera virus. In 2017, the International Committee on Taxonomy of Viruses (ICTV) promoted the family to order Bunyavirales and subsequently greatly expanded its composition by adding multiple families for non-segmented to polysegmented viruses of animals, fungi, plants, and protists. The continued and accelerated discovery of bunyavirals highlighted that an order would not suffice to depict the evolutionary relationships of these viruses. Thus, in April 2024, the order was promoted to class Bunyaviricetes. This class currently includes two major orders, Elliovirales (Cruliviridae, Fimoviridae, Hantaviridae, Peribunyaviridae, Phasmaviridae, Tospoviridae, and Tulasviridae) and Hareavirales (Arenaviridae, Discoviridae, Konkoviridae, Leishbuviridae, Mypoviridae, Nairoviridae, Phenuiviridae, and Wupedeviridae), for hundreds of viruses, many of which are pathogenic for humans and other animals, plants, and fungi.

Endogenous nege-like viral elements in arthropod genomes reveal virus-host coevolution and ancient history of two plant virus families

Negevirus is a recently proposed taxon of arthropod-infecting virus, which is associated with plant viruses of two families (Virgaviridae and Kitaviridae). Nevertheless, the evolutionary history of negevirus-host and its relationship with plant viruses remain poorly understood. Endogenous nege-like viral elements (ENVEs) are ancient nege-like viral sequences integrated into the arthropod genomes, which can serve as the molecular fossil records of previous viral infection. In this study, 292 ENVEs were identified in 150 published arthropod genomes, revealing the evolutionary history of nege-like viruses and two related plant virus families. We discovered three novel and eight strains of nege-like viruses in 11 aphid species. Further analysis indicated that 10 ENVEs were detected in six aphid genomes, and they were divided into four types (ENVE1-ENVE4). Orthologous integration and phylogenetic analyses revealed that nege-like viruses had a history of infection of over 60 My and coexisted with aphid ancestors throughout the Cenozoic Era. Moreover, two nege-like viral proteins (CP and SP24) were highly homologous to those of plant viruses in the families Virgaviridae and Kitaviridae. CP- and SP24-derived ENVEs were widely integrated into numerous arthropod genomes. These results demonstrate that nege-like viruses have a long-term coexistence with arthropod hosts and plant viruses of the two families, Virgaviridae and Kitaviridae, which may have evolved from the nege-like virus ancestor through horizontal virus transfer events. These findings broaden our perspective on the history of viral infection in arthropods and the origins of plant viruses.

IMPORTANCE

Although negevirus is phylogenetically related to plant virus, the evolutionary history of negevirus-host and its relationship with plant virus remain largely unknown. In this study, we used endogenous nege-like viral elements (ENVEs) as the molecular fossil records to investigate the history of nege-like viral infection in arthropod hosts and the evolution of two related plant virus families (Virgaviridae and Kitaviridae). Our results showed the infection of nege-like viruses for over 60 My during the arthropod evolution. ENVEs highly homologous to viral sequences in Virgaviridae and Kitaviridae were present in a wide range of arthropod genomes but were absent in plant genomes, indicating that plant viruses in these two families possibly evolved from the nege-like virus ancestor through cross-species horizontal virus transmission. Our findings provide a new perspective on the virus-host coevolution and the origins of plant viruses.

New lineages of RNA viruses from clinical isolates of Rhizopus microsporus revealed by fragmented and primer-ligated dsRNA sequencing (FLDS) analysis

Rhizopus microsporus is a species in the order Mucorales that is known to cause mucormycosis, but it is poorly understood as a host of viruses. Here, we examined 25 clinical strains of R. microsporus for viral infection with a conventional double-stranded RNA (dsRNA) assay using agarose gel electrophoresis (AGE) and the recently established fragmented and primer-ligated dsRNA sequencing (FLDS) protocol. By AGE, five virus-infected strains were detected. Then, full-length genomic sequences of 12 novel RNA viruses were revealed by FLDS, which were related to the families Mitoviridae, Narnaviridae, and Endornaviridae, ill-defined groups of single-stranded RNA (ssRNA) viruses with similarity to the established families Virgaviridae and Phasmaviridae, and the proposed family "Ambiguiviridae." All the characterized viruses, except a potential phasmavirid with a negative-sense RNA genome, had positive-sense RNA genomes. One virus belonged to a previously established species within the family Mitoviridae, whereas the other 11 viruses represented new species or even new genera. These results show that the fungal pathogen R. microsporus harbors diverse RNA viruses and extend our understanding of the diversity of RNA viruses in the fungal order Mucorales, division Mucoromycota. Identifying RNA viruses from clinical isolates of R. microsporus may expand the repertoire of natural therapeutic agents for mucormycosis in the future.IMPORTANCEThe diversity of mycoviruses in fungal hosts in the division Mucoromycota has been underestimated, mainly within the species Rhizopus microsporus. Only five positive-sense RNA genomes had previously been discovered in this species. Because current sequencing methods poorly complete the termini of genomes, we used fragmented and primer-ligated double-stranded RNA sequencing to acquire the full-length genomes. Eleven novel mycoviruses were detected in this study, including the first negative-sense RNA genome reported in R. microsporus. Our findings extend the understanding of the viral diversity in clinical strains of Mucoromycota, may provide insights into the pathogenesis and ecology of this fungus, and may offer therapeutic options.

ICTV Virus Taxonomy Profile: Phasmaviridae 2024

Phasmaviridae is a family for negative-sense RNA viruses with genomes of about 9.7-15.8 kb. These viruses are maintained in and/or transmitted by insects. Phasmavirids produce enveloped virions containing three single-stranded RNA segments that encode a nucleoprotein (N), a glycoprotein precursor (GPC), and a large (L) protein containing an RNA-directed RNA polymerase (RdRP) domain. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Phasmaviridae, which is available at ictv.global/report/phasmaviridae.

A Systematic Review on the Viruses of Anopheles Mosquitoes: The Potential Importance for Public Health

Anopheles mosquitoes are the vectors of Plasmodium, the etiological agent of malaria. In addition, Anopheles funestus and Anopheles gambiae are the main vectors of the O'nyong-nyong virus. However, research on the viruses carried by Anopheles is scarce; thus, the possible transmission of viruses by Anopheles is still unexplored. This systematic review was carried out to identify studies that report viruses in natural populations of Anopheles or virus infection and transmission in laboratory-reared mosquitoes. The databases reviewed were EBSCO-Host, Google Scholar, Science Direct, Scopus and PubMed. After the identification and screening of candidate articles, a total of 203 original studies were included that reported on a variety of viruses detected in Anopheles natural populations. In total, 161 viruses in 54 species from 41 countries worldwide were registered. In laboratory studies, 28 viruses in 15 Anopheles species were evaluated for mosquito viral transmission capacity or viral infection. The viruses reported in Anopheles encompassed 25 viral families and included arboviruses, probable arboviruses and Insect-Specific Viruses (ISVs). Insights after performing this review include the need for (1) a better understanding of Anopheles-viral interactions, (2) characterizing the Anopheles virome-considering the public health importance of the viruses potentially transmitted by Anopheles and the significance of finding viruses with biological control activity-and (3) performing virological surveillance in natural populations of Anopheles, especially in the current context of environmental modifications that may potentiate the expansion of the Anopheles species distribution.

Phlebovirus diversity in ticks from livestock in arid ecologies in Kenya

Phleboviruses are emerging pathogens of public health importance. However, their association with ticks is poorly described, particularly in Africa. Here, adult ticks infesting cattle, goats and sheep were collected in two dryland pastoralist ecosystems of Kenya (Baringo and Kajiado counties) and were screened for infection with phleboviruses. Ticks mainly belonged to the species Rhipicephalus appendiculatus, Hyalomma impeltatum, and Hyalomma rufipes. A fragment of the RNA-dependent RNA polymerase (RdRp) gene was identified in thirty of 671 tick pools, of which twenty-nine were from livestock sampled in Baringo county. Phylogenetic analyses revealed that twenty-five sequences were falling in three clades within the group of tick-associated phleboviruses. The sequences of the three clades showed nucleotide distances 8%, 19% and 22%, respectively, to previously known viruses suggesting that these sequence fragments may belong to three distinct viruses. Viruses of the group of tick-associated phleboviruses have been found in several countries and continents but so far have not been associated with disease in humans or animals. In addition, five sequences were found to group with the sandfly-associated phleboviruses Bogoria virus, Perkerra virus and Ntepes virus recently detected in the same region. Further studies are needed to investigate the transmission and maintenance cycles of these viruses, as well as to assess their potential to infect vertebrates.

The diversity of endogenous viral elements in insects

We provide an overview of the currently known diversity of viral sequences integrated into insect genomes. Such endogenous viral elements (EVE) have so far been annotated in at least eight insect orders and can be assigned to at least three families of large double-stranded (ds) DNA viruses, at least 22 families of RNA viruses, and three families of single-stranded DNA viruses. The study of these EVE has already produced important insights into insect-virus interactions, including the discovery of a new form of adaptive antiviral immunity. Insect EVE diversity will continue to increase as new insect genomes and exogenous viruses are sequenced, which will continue to make paleovirology a vibrant research field in this group of animals in the years to come.

The Bunyavirales: The Plant-Infecting Counterparts

Negative-strand (-) RNA viruses (NSVs) comprise a large and diverse group of viruses that are generally divided in those with non-segmented and those with segmented genomes. Whereas most NSVs infect animals and humans, the smaller group of the plant-infecting counterparts is expanding, with many causing devastating diseases worldwide, affecting a large number of major bulk and high-value food crops. In 2018, the taxonomy of segmented NSVs faced a major reorganization with the establishment of the order Bunyavirales. This article overviews the major plant viruses that are part of the order, i.e., orthospoviruses (Tospoviridae), tenuiviruses (Phenuiviridae), and emaraviruses (Fimoviridae), and provides updates on the more recent ongoing research. Features shared with the animal-infecting counterparts are mentioned, however, special attention is given to their adaptation to plant hosts and vector transmission, including intra/intercellular trafficking and viral counter defense to antiviral RNAi.

Recent Advances in Bunyavirus Glycoprotein Research: Precursor Processing, Receptor Binding and Structure

The Bunyavirales order accommodates related viruses (bunyaviruses) with segmented, linear, single-stranded, negative- or ambi-sense RNA genomes. Their glycoproteins form capsomeric projections or spikes on the virion surface and play a crucial role in virus entry, assembly, morphogenesis. Bunyavirus glycoproteins are encoded by a single RNA segment as a polyprotein precursor that is co- and post-translationally cleaved by host cell enzymes to yield two mature glycoproteins, Gn and Gc (or GP1 and GP2 in arenaviruses). These glycoproteins undergo extensive N-linked glycosylation and despite their cleavage, remain associated to the virion to form an integral transmembrane glycoprotein complex. This review summarizes recent advances in our understanding of the molecular biology of bunyavirus glycoproteins, including their processing, structure, and known interactions with host factors that facilitate cell entry.

Viruses in Extreme Environments, Current Overview, and Biotechnological Potential

Virus research has advanced significantly since the discovery of the tobacco mosaic virus (TMV), the characterization of its infection mechanisms and the factors that determine their pathogenicity. However, most viral research has focused on pathogenic viruses to humans, animals and plants, which represent only a small fraction in the virosphere. As a result, the role of most viral genes, and the mechanisms of coevolution between mutualistic viruses, their host and their environment, beyond pathogenicity, remain poorly understood. This review focuses on general aspects of viruses that interact with extremophile organisms, characteristics and examples of mechanisms of adaptation. Finally, this review provides an overview on how knowledge of extremophile viruses sheds light on the application of new tools of relevant use in modern molecular biology, discussing their value in a biotechnological context.

Stability of the Virome in Lab- and Field-Collected Aedes albopictus Mosquitoes across Different Developmental Stages and Possible Core Viruses in the Publicly Available Virome Data of Aedes Mosquitoes

Aedes mosquitoes can efficiently transmit many pathogenic arboviruses, placing a great burden on public health worldwide. In addition, they also carry a number of insect-specific viruses (ISVs), and it was recently suggested that some of these ISVs might form a stable species-specific "core virome" in mosquito populations. However, little is known about such a core virome in laboratory colonies and if it is present across different developmental stages. In this study, we compared the viromes in eggs, larvae, pupae, and adults of Aedes albopictus mosquitoes collected from a lab colony and compared each to the virome of different developmental stages collected in the field. The virome in lab-derived A. albopictus was very stable across all stages, consistent with a vertical transmission route of these viruses, and formed a possible "vertically transmitted core virome." The different stages of field-collected A. albopictus mosquitoes also contained this stable vertically transmitted core virome, as well as another set of viruses (e.g., viruses distantly related to Guadeloupe mosquito virus, Hubei virga-like virus 2, and Sarawak virus) shared by mosquitoes across different stages, which might represent an "environment-derived core virome." To further study this core set of ISVs, we screened 48 publicly available SRA viral metagenomic data sets of mosquitoes belonging to the genus Aedes, showing that some of the identified ISVs were identified in the majority of SRAs and providing further evidence supporting the core-virome concept.IMPORTANCE Our study revealed that the virome was very stable across all developmental stages of both lab-derived and field-collected Aedes albopictus The data representing the core virome in lab A. albopictus proved the vertical transmission route of these viruses, forming a "vertically transmitted core virome." Field mosquitoes also contained this stable vertically transmitted core virome as well as additional viruses, which probably represented "environment-derived core virome" and which therefore were less stable over time and geography. By further screening publicly available SRA viral metagenomic data sets from mosquitoes belonging to the genus Aedes, some of the identified core ISVs were shown to be present in the majority of SRAs, such as Phasi Charoen-like phasivirus and Guadeloupe mosquito virus. How these core ISVs influence the biology of the mosquito host and arbovirus infection and evolution deserves to be further explored.

High prevalence of Phasi Charoen-like virus from wild-caught Aedes aegypti in Grenada, W.I. as revealed by metagenomic analysis

Arboviruses cause diseases of significant global health concerns. Interactions between mosquitoes and their microbiota as well as the important role of this interaction in the mosquito's capacity to harbor and transmit pathogens have emerged as important fields of research. Aedes aegypti is one of the most abundant mosquitoes in many geographic locations, a vector capable of transmitting a number of arboviruses such as dengue and Zika. Currently, there are few studies on the metavirome of this mosquito particularly in the Americas. This study analyzes the metavirome of A. aegypti from Grenada, a Caribbean nation with tropical weather, abundant A. aegypti, and both endemic and arboviral pathogens transmitted by this mosquito. Between January and December 2018, 1152 mosquitoes were collected from six semi-rural locations near houses in St. George Parish, Grenada, by weekly trapping using BG-Sentinel traps. From these, 300 A. aegypti were selected for analysis. The metavirome was analyzed using the Illumina HiSeq 1500 for deep sequencing. The generation sequencing library construction protocol used was NuGEN Universal RNA with an average read length of 125 bp. Reads were mapped to the A. aegypti assembly. Non-mosquito reads were analyzed using the tools FastViromeExplorer. The NCBI total virus, RNA virus, and eukaryotic virus databases were used as references. The metagenomic comparison analysis showed that the most abundant virus-related reads among all databases and assemblies was Phasi Charoen-like virus. The Phasi Charoen-like virus results are in agreement to other studies in America, Asia and Australia. Further studies using wild-caught mosquitoes is needed to assess the impact of this insect-specific virus on the A. aegypti lifecycle and vector capacity.

Arboviruses and the Challenge to Establish Systemic and Persistent Infections in Competent Mosquito Vectors: The Interaction With the RNAi Mechanism

Arboviruses are capable to establish long-term persistent infections in mosquitoes that do not affect significantly the physiology of the insect vectors. Arbovirus infections are controlled by the RNAi machinery via the production of viral siRNAs and the formation of RISC complexes targeting viral genomes and mRNAs. Engineered arboviruses that contain cellular gene sequences can therefore be transformed to "viral silencing vectors" for studies of gene function in reverse genetics approaches. More specifically, "ideal" viral silencing vectors must be competent to induce robust RNAi effects while other interactions with the host immune system should be kept at a minimum to reduce non-specific effects. Because of their inconspicuous nature, arboviruses may approach the "ideal" viral silencing vectors in insects and it is therefore worthwhile to study the mechanisms by which the interactions with the RNAi machinery occur. In this review, an analysis is presented of the antiviral RNAi response in mosquito vectors with respect to the major types of arboviruses (alphaviruses, flaviviruses, bunyaviruses, and others). With respect to antiviral defense, the exo-RNAi pathway constitutes the major mechanism while the contribution of both miRNAs and viral piRNAs remains a contentious issue. However, additional mechanisms exist in mosquitoes that are capable to enhance or restrict the efficiency of viral silencing vectors such as the amplification of RNAi effects by DNA forms, the existence of incorporated viral elements in the genome and the induction of a non-specific systemic response by Dicer-2. Of significance is the observation that no major "viral suppressors of RNAi" (VSRs) seem to be encoded by arboviral genomes, indicating that relatively tight control of the activity of the RNA-dependent RNA polymerase (RdRp) may be sufficient to maintain the persistent character of arbovirus infections. Major strategies for improvement of viral silencing vectors therefore are proposed to involve engineering of VSRs and modifying of the properties of the RdRp. Because of safety issues (pathogen status), however, arbovirus-based silencing vectors are not well suited for practical applications, such as RNAi-based mosquito control. In that case, related mosquito-specific viruses that also establish persistent infections and may cause similar RNAi responses may represent a valuable alternative solution.

Evolutionary persistence of insect bunyavirus infection despite host acquisition and expression of the viral nucleoprotein gene

How insects combat RNA virus infection is a subject of intensive research owing to its importance in insect health, virus evolution, and disease transmission. In recent years, a pair of potentially linked phenomena have come to light as a result of this work-first, the pervasive production of viral DNA from exogenous nonretroviral RNA in infected individuals, and second, the widespread distribution of nonretroviral integrated RNA virus sequences (NIRVs) in the genomes of diverse eukaryotes. The evolutionary consequences of NIRVs for viruses are unclear and the field would benefit from studies of natural virus infections co-occurring with recent integrations, an exceedingly rare circumstance in the literature. Here, we provide evidence that a novel insect-infecting phasmavirus (Order Bunyavirales) has been persisting in a phantom midge host, Chaoborus americanus, for millions of years. Interestingly, the infection persists despite the host's acquisition (during the Pliocene), fixation, and expression of the viral nucleoprotein gene. We show that virus prevalence and geographic distribution are high and broad, comparable to the host-specific infections reported in other phantom midges. Short-read mapping analyses identified a lower abundance of the nucleoprotein-encoding genome segment in this virus relative to related viruses. Finally, the novel virus has facilitated the first substitution rate estimation for insect-infecting phasmaviruses. Over a period of approximately 16 million years, we find rates of (0.6 - 1.6) × 10^-7 substitutions per site per year in protein coding genes, extraordinarily low for negative-sense RNA viruses, but consistent with the few estimates produced over comparable evolutionary timescales.

A Simple Method to Detect Candidate Overlapping Genes in Viruses Using Single Genome Sequences

Overlapping genes in viruses maximize the coding capacity of their genomes and allow the generation of new genes without major increases in genome size. Despite their importance, the evolution and function of overlapping genes are often not well understood, in part due to difficulties in their detection. In addition, most bioinformatic approaches for the detection of overlapping genes require the comparison of multiple genome sequences that may not be available in metagenomic surveys of virus biodiversity. We introduce a simple new method for identifying candidate functional overlapping genes using single virus genome sequences. Our method uses randomization tests to estimate the expected length of open reading frames and then identifies overlapping open reading frames that significantly exceed this length and are thus predicted to be functional. We applied this method to 2548 reference RNA virus genomes and find that it has both high sensitivity and low false discovery for genes that overlap by at least 50 nucleotides. Notably, this analysis provided evidence for 29 previously undiscovered functional overlapping genes, some of which are coded in the antisense direction suggesting there are limitations in our current understanding of RNA virus replication.

An insight into the sialotranscriptome and virome of Amazonian anophelines

Background

Saliva of mosquitoes contains anti-platelet, anti-clotting, vasodilatory, anti-complement and anti-inflammatory substances that help the blood feeding process. The salivary polypeptides are at a fast pace of evolution possibly due to their relative lack of structural constraint and possibly also by positive selection on their genes leading to evasion of host immune pressure.

Results

In this study, we used deep mRNA sequence to uncover for the first time the sialomes of four Amazonian anophelines species (Anopheles braziliensis, A. marajorara, A. nuneztovari and A. triannulatus) and extend the knowledge of the A. darlingi sialome. Two libraries were generated from A. darlingi mosquitoes, sampled from two localities separated ~ 1100 km apart. A total of 60,016 sequences were submitted to GenBank, which will help discovery of novel pharmacologically active polypeptides and the design of specific immunological markers of mosquito exposure. Additionally, in these analyses we identified and characterized novel phasmaviruses and anpheviruses associated to the sialomes of A. triannulatus, A. marajorara and A. darlingi species.

Conclusions

Besides their pharmacological properties, which may be exploited for the development of new drugs (e.g. anti-thrombotics), salivary proteins of blood feeding arthropods may be turned into tools to prevent and/or better control vector borne diseases; for example, through the development of vaccines or biomarkers to evaluate human exposure to vector bites. The sialotranscriptome study reported here provided novel data on four New World anopheline species and allowed to extend our knowledge on the salivary repertoire of A. darlingi. Additionally, we discovered novel viruses following analysis of the transcriptomes, a procedure that should become standard within future RNAseq studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5545-0) contains supplementary material, which is available to authorized users.

Insect-specific virus evolution and potential effects on vector competence

The advancement in high-throughput sequencing technology and bioinformatics tools has spurred a new age of viral discovery. Arthropods is the largest group of animals and has shown to be a major reservoir of different viruses, including a group known as insect-specific viruses (ISVs). The majority of known ISVs have been isolated from mosquitoes and shown to belong to viral families associated with animal arbovirus pathogens, such as Flaviviridae, Togaviridae and Phenuiviridae. These insect-specific viruses have a strict tropism and are unable to replicate in vertebrate cells, these properties are interesting for many reasons. One is that these viruses could potentially be utilised as biocontrol agents using a similar strategy as for Wolbachia. Mosquitoes infected with the viral agent could have inferior vectorial capacity of arboviruses resulting in a decrease of circulating arboviruses of public health importance. Moreover, insect-specific viruses are thought to be ancestral to arboviruses and could be used to study the evolution of the switch from single-host to dual-host. In this review, we discuss new discoveries and hypothesis in the field of arboviruses and insect-specific viruses.

Firefly genomes illuminate parallel origins of bioluminescence in beetles

Fireflies and their luminous courtships have inspired centuries of scientific study. Today firefly luciferase is widely used in biotechnology, but the evolutionary origin of bioluminescence within beetles remains unclear. To shed light on this long-standing question, we sequenced the genomes of two firefly species that diverged over 100 million-years-ago: the North American Photinus pyralis and Japanese Aquatica lateralis. To compare bioluminescent origins, we also sequenced the genome of a related click beetle, the Caribbean Ignelater luminosus, with bioluminescent biochemistry near-identical to fireflies, but anatomically unique light organs, suggesting the intriguing hypothesis of parallel gains of bioluminescence. Our analyses support independent gains of bioluminescence in fireflies and click beetles, and provide new insights into the genes, chemical defenses, and symbionts that evolved alongside their luminous lifestyle.

Taxonomy of the family Arenaviridae and the order Bunyavirales: update 2018

In 2018, the family Arenaviridae was expanded by inclusion of 1 new genus and 5 novel species. At the same time, the recently established order Bunyavirales was expanded by 3 species. This article presents the updated taxonomy of the family Arenaviridae and the order Bunyavirales as now accepted by the International Committee on Taxonomy of Viruses (ICTV) and summarizes additional taxonomic proposals that may affect the order in the near future.

Bunyavirales ribonucleoproteins: the viral replication and transcription machinery

The Bunyavirales order is one of the largest groups of segmented negative-sense single-stranded RNA viruses, which includes many pathogenic strains that cause severe human diseases. The RNA segments of the bunyavirus genome are separately encapsidated by multiple copies of nucleoprotein (N), and both termini of each N-encapsidated genomic RNA segment bind to one copy of the viral L polymerase protein. The viral genomic RNA, N and L protein together form the ribonucleoprotein (RNP) complex that constitutes the molecular machinery for viral genome replication and transcription. Recently, breakthroughs have been achieved in understanding the architecture of bunyavirus RNPs with the determination of the atomic structures of the N and L proteins from various members of this order. In this review, we discuss the structures and functions of these bunyavirus RNP components, as well as viral genome replication and transcription mechanisms.

Insect Virus Discovery by Metagenomic and Cell Culture-Based Approaches

Insects are the most abundant and diverse group of animals on earth, but our knowledge of their viruses is biased toward insect-borne viruses that cause disease in plants, animals, or humans. Recent metagenomic studies and systematic surveys of viruses in wild-caught insects have identified an unanticipated large repertoire of novel viruses and viral sequences. These include new members of existing clades, new clades, and even entirely new virus families. These studies greatly expand the known virosphere in insects, provide opportunities to study virus-host interactions, and generate new insights into virus evolution. In this chapter, we discuss the methods used to identify novel viruses in insects and highlight some notable surprises arising from these studies.

High-Resolution Metatranscriptomics Reveals the Ecological Dynamics of Mosquito-Associated RNA Viruses in Western Australia

Mosquitoes harbor a high diversity of RNA viruses, including many that impact human health. Despite a growing effort to describe the extent and nature of the mosquito virome, little is known about how these viruses persist, spread, and interact with both their hosts and other microbes. To address this issue we performed a metatranscriptomics analysis of 12 Western Australian mosquito populations structured by species and geographic location. Our results identified the complete genomes of 24 species of RNA viruses from a diverse range of viral families and orders, among which 19 are newly described. Comparisons of viromes revealed a striking difference between the two mosquito genera, with viromes of mosquitoes of the Aedes genus exhibiting substantially less diversity and lower abundances than those of mosquitoes of the Culex genus, within which the viral abundance reached 16.87% of the total non-rRNA. In addition, there was little overlap in viral diversity between the two genera, although the viromes were very similar among the three Culex species studied, suggesting that the host taxon plays a major role in structuring virus diversity. In contrast, we found no evidence that geographic location played a major role in shaping RNA virus diversity, and several viruses discovered here exhibited high similarity (95 to 98% nucleotide identity) to those from Indonesia and China. Finally, using abundance-level and phylogenetic relationships, we were able to distinguish potential mosquito viruses from those present in coinfecting bacteria, fungi, and protists. In sum, our metatranscriptomics approach provides important insights into the ecology of mosquito RNA viruses.IMPORTANCE Studies of virus ecology have generally focused on individual viral species. However, recent advances in bulk RNA sequencing make it possible to utilize metatranscriptomic approaches to reveal both complete virus diversity and the relative abundance of these viruses. We used such a metatranscriptomic approach to determine key aspects of the ecology of mosquito viruses in Western Australia. Our results show that RNA viruses are some of the most important components of the mosquito transcriptome, and we identified 19 new virus species from a diverse set of virus families. A key result was that host genetic background plays a more important role in shaping virus diversity than sampling location, with Culex species harboring more viruses at higher abundance than those from Aedes mosquitoes.

Nonretroviral integrated RNA viruses in arthropod vectors: an occasional event or something more?

With few exceptions, all arthropod-borne viruses (arboviruses) are nonretroviral RNA viruses (NRVs). Despite NRVs do not encode reverse transcriptases and integrases, NRVs-DNA fragments are detected in mosquito cells and mosquitoes at early stages of infection as episomal DNA forms. Additionally, next generation sequencing and bioinformatics analyses have convincingly shown NRVs-vDNA integrated in vector genomes. We hypothesize vDNA role may be linked to host immunity and viral persistence. Key questions remain about nonretroviral integrated RNA virus sequences (NIRVS) in mosquitoes such as what is driving vDNA synthesis from NRVs, how does integration occur and what is their biological function. Here we review current knowledge about NIRVS highlighting connections with host immunity and virus-vector co-evolution and we suggest directions for future research.

Mosquito-specific and mosquito-borne viruses: evolution, infection, and host defense

Recent virus discovery programs have identified an extensive reservoir of viruses in arthropods. It is thought that arthropod viruses, including mosquito-specific viruses, are ancestral to vertebrate-pathogenic arboviruses. Mosquito-specific viruses are restricted in vertebrate cells at multiple levels, including entry, RNA replication, assembly, and by the inability to replicate at high temperatures. Moreover, it is likely that the vertebrate immune system suppresses replication of these viruses. The evolution from single to dual-host tropism may also require changes in the course of infection in the mosquito host. In this review we explore the adaptive changes required for a switch from a mosquito-specific to a mosquito-borne transmission cycle.

Vertebrate Reservoirs of Arboviruses: Myth, Synonym of Amplifier, or Reality?

The rapid succession of the pandemic of arbovirus diseases, such as dengue, West Nile fever, chikungunya, and Zika fever, has intensified research on these and other arbovirus diseases worldwide. Investigating the unique mode of vector-borne transmission requires a clear understanding of the roles of vertebrates. One major obstacle to this understanding is the ambiguity of the arbovirus definition originally established by the World Health Organization. The paucity of pertinent information on arbovirus transmission at the time contributed to the notion that vertebrates played the role of reservoir in the arbovirus transmission cycle. Because this notion is a salient feature of the arbovirus definition, it is important to reexamine its validity. This review addresses controversial issues concerning vertebrate reservoirs and their role in arbovirus persistence in nature, examines the genesis of the problem from a historical perspective, discusses various unresolved issues from multiple points of view, assesses the present status of the notion in light of current knowledge, and provides options for a solution to resolve the issue.

Comparative genomics shows that viral integrations are abundant and express piRNAs in the arboviral vectors Aedes aegypti and Aedes albopictus

Background

Arthropod-borne viruses (arboviruses) transmitted by mosquito vectors cause many important emerging or resurging infectious diseases in humans including dengue, chikungunya and Zika. Understanding the co-evolutionary processes among viruses and vectors is essential for the development of novel transmission-blocking strategies. Episomal viral DNA fragments are produced from arboviral RNA upon infection of mosquito cells and adults. Additionally, sequences from insect-specific viruses and arboviruses have been found integrated into mosquito genomes.

Results

We used a bioinformatic approach to analyse the presence, abundance, distribution, and transcriptional activity of integrations from 425 non-retroviral viruses, including 133 arboviruses, across the presently available 22 mosquito genome sequences. Large differences in abundance and types of viral integrations were observed in mosquito species from the same region. Viral integrations are unexpectedly abundant in the arboviral vector species Aedes aegypti and Ae. albopictus, in which they are approximately ~10-fold more abundant than in other mosquito species analysed. Additionally, viral integrations are enriched in piRNA clusters of both the Ae. aegypti and Ae. albopictus genomes and, accordingly, they express piRNAs, but not siRNAs.

Conclusions

Differences in the number of viral integrations in the genomes of mosquito species from the same geographic area support the conclusion that integrations of viral sequences is not dependent on viral exposure, but that lineage-specific interactions exist. Viral integrations are abundant in Ae. aegypti and Ae. albopictus, and represent a thus far underappreciated component of their genomes. Additionally, the genome locations of viral integrations and their production of piRNAs indicate a functional link between viral integrations and the piRNA pathway. These results greatly expand the breadth and complexity of small RNA-mediated regulation and suggest a role for viral integrations in antiviral defense in these two mosquito species.

Electronic supplementary material

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

Bunyaviridae RdRps: structure, motifs, and RNA synthesis machinery

Bunyaviridae family is the largest and most diverse family of RNA viruses. It has more than 350 members divided into five genera: Orthobunyavirus, Phlebovirus, Nairovirus, Hantavirus, and Tospovirus. They are present in the five continents, causing recurrent epidemics, epizootics, and considerable agricultural loss. The genome of bunyaviruses is divided into three segments of negative single-stranded RNA according to their relative size: L (Large), M (Medium) and S (Small) segment. Bunyaviridae RNA-dependent RNA polymerase (RdRp) is encoded by the L segment, and is in charge of the replication and transcription of the viral RNA in the cytoplasm of the infected cell. Viral RdRps share a characteristic right hand-like structure with three subdomains: finger, palm, and thumb subdomains that define the formation of the catalytic cavity. In addition to the N-terminal endonuclease domain, eight conserved motifs (A-H) have been identified in the RdRp of Bunyaviridae. In this review, we have summarized the recent insights from the structural and functional studies of RdRp to understand the roles of different motifs shared by RdRps, the mechanism of viral RNA replication, genome segment packaging by the nucleoprotein, cap-snatching, mRNA transcription, and other RNA mechanisms of bunyaviruses.

The Envelope Proteins of the Bunyavirales

The Bunyavirales Order encompasses nine families of enveloped viruses containing a single-stranded negative-sense RNA genome divided into three segments. The small (S) and large (L) segments encode proteins participating in genome replication in the infected cell cytoplasm. The middle (M) segment encodes the viral glycoproteins Gn and Gc, which are derived from a precursor polyprotein by host cell proteases. Entry studies are available only for a few viruses in the Order, and in each case they were shown to enter cells via receptor-mediated endocytosis. The acidic endosomal pH triggers the fusion of the viral envelope with the membrane of an endosome. Structural studies on two members of this Order, the phleboviruses and the hantaviruses, have shown that the membrane fusion protein Gc displays a class II fusion protein fold and is homologous to its counterparts in flaviviruses and alphaviruses, which are positive-sense, single-stranded RNA viruses. We analyze here recent data on the structure and function of the structure of the phlebovirus Gc and hantavirus Gn and Gc glycoproteins, and extrapolate common features identified in the amino acid sequences to understand also the structure and function of their counterparts in other families of the Bunyavirales Order. Our analysis also identified clear structural homology between the hantavirus Gn and alphavirus E2 glycoproteins, which make a heterodimer with the corresponding fusion proteins Gc and E1, respectively, revealing that not only the fusion protein has been conserved across viral families.

Unexpected differences in the population genetics of phasmavirids (Bunyavirales) from subarctic ponds

Little is known of the evolution of RNA viruses in aquatic systems. Here, we assess the genetic connectivity of two bunyaviruses (Kigluaik phantom orthophasmavirus or KIGV and Nome phantom orthophasmavirus or NOMV) with zooplanktonic hosts from subarctic ponds. We expected weak genetic structure among populations as the hosts (phantom midges) have a terrestrial winged dispersal stage. To test whether their respective viruses mirror this structure, we collected and analyzed population datasets from 21 subarctic freshwater ponds and obtained sequences from all four genes in the viral genomes. Prevalence averaged 66 per cent for 514 host specimens and was not significantly different between recently formed thaw ponds and glacial ponds. Unexpectedly, KIGV from older ponds showed pronounced haplotype divergence with little evidence of genetic connectivity. However, KIGV populations from recent thaw ponds appeared to be represented by a closely related haplotype group, perhaps indicating a genotypic dispersal bias. Unlike KIGV, NOMV had modest structure and diversity in recently formed thaw ponds. For each virus, we found elevated genetic diversity relative to the host, but similar population structures to the host. Our results suggest that non-random processes such as virus-host interactions, genotypic bias, and habitat effects differ among polar aquatic RNA viruses.

Unbiased RNA Shotgun Metagenomics in Social and Solitary Wild Bees Detects Associations with Eukaryote Parasites and New Viruses

The diversity of eukaryote organisms and viruses associated with wild bees remains poorly characterized in contrast to the well-documented pathosphere of the western honey bee, Apis mellifera. Using a deliberate RNA shotgun metagenomic sequencing strategy in combination with a dedicated bioinformatics workflow, we identified the (micro-)organisms and viruses associated with two bumble bee hosts, Bombus terrestris and Bombus pascuorum, and two solitary bee hosts, Osmia cornuta and Andrena vaga. Ion Torrent semiconductor sequencing generated approximately 3.8 million high quality reads. The most significant eukaryote associations were two protozoan, Apicystis bombi and Crithidia bombi, and one nematode parasite Sphaerularia bombi in bumble bees. The trypanosome protozoan C. bombi was also found in the solitary bee O. cornuta. Next to the identification of three honey bee viruses Black queen cell virus, Sacbrood virus and Varroa destructor virus-1 and four plant viruses, we describe two novel RNA viruses Scaldis River bee virus (SRBV) and Ganda bee virus (GABV) based on their partial genomic sequences. The novel viruses belong to the class of negative-sense RNA viruses, SRBV is related to the order Mononegavirales whereas GABV is related to the family Bunyaviridae. The potential biological role of both viruses in bees is discussed in the context of recent advances in the field of arthropod viruses. Further, fragmentary sequence evidence for other undescribed viruses is presented, among which a nudivirus in O. cornuta and an unclassified virus related to Chronic bee paralysis virus in B. terrestris. Our findings extend the current knowledge of wild bee parasites in general and addsto the growing evidence of unexplored arthropod viruses in valuable insects.

Redefining the invertebrate RNA virosphere

Current knowledge of RNA virus biodiversity is both biased and fragmentary, reflecting a focus on culturable or disease-causing agents. Here we profile the transcriptomes of over 220 invertebrate species sampled across nine animal phyla and report the discovery of 1,445 RNA viruses, including some that are sufficiently divergent to comprise new families. The identified viruses fill major gaps in the RNA virus phylogeny and reveal an evolutionary history that is characterized by both host switching and co-divergence. The invertebrate virome also reveals remarkable genomic flexibility that includes frequent recombination, lateral gene transfer among viruses and hosts, gene gain and loss, and complex genomic rearrangements. Together, these data present a view of the RNA virosphere that is more phylogenetically and genomically diverse than that depicted in current classification schemes and provide a more solid foundation for studies in virus ecology and evolution.

Novel Viruses Isolated from Mosquitoes in Pantanal, Brazil

Genomic sequences are described from five novel viruses and divergent strains of Brejeira and Guaico Culex viruses from mosquitoes collected in Pantanal, Brazil, in 2010.

Resistance to Tospoviruses in Vegetable Crops: Epidemiological and Molecular Aspects

During the past three decades, the economic impact of tospoviruses has increased, causing high yield losses in a variety of crops and ornamentals. Owing to the difficulty in combating thrips vectors with insecticides, the best way to limit/prevent tospovirus-induced diseases involves a management strategy that includes virus resistance. This review briefly presents current tospovirus taxonomy, diversity, molecular biology, and cytopathology as an introduction to a more extensive description of the two main resistance genes employed against tospoviruses: the Sw5 gene in tomato and the Tsw in pepper. Natural and experimental resistance-breaking (RB) isolates allowed the identification of the viral avirulence protein triggering each of the two resistance gene products; epidemiology of RB isolates is discussed to reinforce the need for allelic variants and the need to search for new/alternative resistance genes. Ongoing efforts for alternative resistance strategies are described not only for Tomato spotted wilt virus (TSWV) in pepper and tomato but also for other vegetable crops heavily impacted by tospoviruses.

Evolutionary origin of pathogenic arthropod-borne viruses-a case study in the family Bunyaviridae

Arthropod-borne viruses have a dual-host tropism and their transmission requires the infection of two disparate hosts, arthropods and vertebrates. Arboviruses occur in several RNA families that also contain viruses with a monotropism for either arthropods or vertebrates. The evolutionary origin of the dual-host tropism of arboviruses was recently identified for the family Bunyaviridae. Bunyaviruses were suggested to have evolved from viruses that are restricted to arthropods as hosts (arthropod-specific viruses). Additional findings of an immense genetic diversity of bunyaviruses in non-blood feeding arthropods support the hypothesis of an arthropod origin of vertebrate-pathogenic bunyaviruses.

Twenty-Five New Viruses Associated with the Drosophilidae (Diptera)

Drosophila melanogaster is an important laboratory model for studies of antiviral immunity in invertebrates, and Drosophila species provide a valuable system to study virus host range and host switching. Here, we use metagenomic RNA sequencing of about 1600 adult flies to discover 25 new RNA viruses associated with six different drosophilid hosts in the wild. We also provide a comprehensive listing of viruses previously reported from the Drosophilidae. The new viruses include Iflaviruses, Rhabdoviruses, Nodaviruses, and Reoviruses, and members of unclassified lineages distantly related to Negeviruses, Sobemoviruses, Poleroviruses, Flaviviridae, and Tombusviridae. Among these are close relatives of Drosophila X virus and Flock House virus, which we find in association with wild Drosophila immigrans. These two viruses are widely used in experimental studies but have not been previously reported to naturally infect Drosophila. Although we detect no new DNA viruses, in D. immigrans and Drosophila obscura, we identify sequences very closely related to Armadillidium vulgare iridescent virus (Invertebrate iridescent virus 31), bringing the total number of DNA viruses found in the Drosophilidae to three.

Genomic Characterization of the Genus Nairovirus (Family Bunyaviridae)

Nairovirus, one of five bunyaviral genera, includes seven species. Genomic sequence information is limited for members of the Dera Ghazi Khan, Hughes, Qalyub, Sakhalin, and Thiafora nairovirus species. We used next-generation sequencing and historical virus-culture samples to determine 14 complete and nine coding-complete nairoviral genome sequences to further characterize these species. Previously unsequenced viruses include Abu Mina, Clo Mor, Great Saltee, Hughes, Raza, Sakhalin, Soldado, and Tillamook viruses. In addition, we present genomic sequence information on additional isolates of previously sequenced Avalon, Dugbe, Sapphire II, and Zirqa viruses. Finally, we identify Tunis virus, previously thought to be a phlebovirus, as an isolate of Abu Hammad virus. Phylogenetic analyses indicate the need for reassignment of Sapphire II virus to Dera Ghazi Khan nairovirus and reassignment of Hazara, Tofla, and Nairobi sheep disease viruses to novel species. We also propose new species for the Kasokero group (Kasokero, Leopards Hill, Yogue viruses), the Ketarah group (Gossas, Issyk-kul, Keterah/soft tick viruses) and the Burana group (Wēnzhōu tick virus, Huángpí tick virus 1, Tǎchéng tick virus 1). Our analyses emphasize the sister relationship of nairoviruses and arenaviruses, and indicate that several nairo-like viruses (Shāyáng spider virus 1, Xīnzhōu spider virus, Sānxiá water strider virus 1, South Bay virus, Wǔhàn millipede virus 2) require establishment of novel genera in a larger nairovirus-arenavirus supergroup.

Characterization of Viral Communities of Biting Midges and Identification of Novel Thogotovirus Species and Rhabdovirus Genus

More than two thirds of emerging viruses are of zoonotic origin, and among them RNA viruses represent the majority. Ceratopogonidae (genus Culicoides) are well-known vectors of several viruses responsible for epizooties (bluetongue, epizootic haemorrhagic disease, etc.). They are also vectors of the only known virus infecting humans: the Oropouche virus. Female midges usually feed on a variety of hosts, leading to possible transmission of emerging viruses from animals to humans. In this context, we report here the analysis of RNA viral communities of Senegalese biting midges using next-generation sequencing techniques as a preliminary step toward the identification of potential viral biohazards. Sequencing of the RNA virome of three pools of Culicoides revealed the presence of a significant diversity of viruses infecting plants, insects and mammals. Several novel viruses were detected, including a novel Thogotovirus species, related but genetically distant from previously described tick-borne thogotoviruses. Novel rhabdoviruses were also detected, possibly constituting a novel Rhabdoviridae genus, and putatively restricted to insects. Sequences related to the major viruses transmitted by Culicoides, i.e., African horse sickness, bluetongue and epizootic haemorrhagic disease viruses were also detected. This study highlights the interest in monitoring the emergence and circulation of zoonoses and epizooties using their arthropod vectors.

Diverse Array of New Viral Sequences Identified in Worldwide Populations of the Asian Citrus Psyllid (Diaphorina citri) Using Viral Metagenomics

The Asian citrus psyllid, Diaphorina citri, is the natural vector of the causal agent of Huanglongbing (HLB), or citrus greening disease. Together; HLB and D. citri represent a major threat to world citrus production. As there is no cure for HLB, insect vector management is considered one strategy to help control the disease, and D. citri viruses might be useful. In this study, we used a metagenomic approach to analyze viral sequences associated with the global population of D. citri. By sequencing small RNAs and the transcriptome coupled with bioinformatics analysis, we showed that the virus-like sequences of D. citri are diverse. We identified novel viral sequences belonging to the picornavirus superfamily, the Reoviridae, Parvoviridae, and Bunyaviridae families, and an unclassified positive-sense single-stranded RNA virus. Moreover, a Wolbachia prophage-related sequence was identified. This is the first comprehensive survey to assess the viral community from worldwide populations of an agricultural insect pest. Our results provide valuable information on new putative viruses, some of which may have the potential to be used as biocontrol agents.

IMPORTANCE

Insects have the most species of all animals, and are hosts to, and vectors of, a great variety of known and unknown viruses. Some of these most likely have the potential to be important fundamental and/or practical resources. In this study, we used high-throughput next-generation sequencing (NGS) technology and bioinformatics analysis to identify putative viruses associated with Diaphorina citri, the Asian citrus psyllid. D. citri is the vector of the bacterium causing Huanglongbing (HLB), currently the most serious threat to citrus worldwide. Here, we report several novel viral sequences associated with D. citri.

Arbovirus-mosquito interactions: RNAi pathway

Arthropod-borne (arbo) viruses infect hematophagous arthropods (vectors) to maintain virus transmission between vertebrate hosts. The mosquito vector actively controls arbovirus infection to minimize its fitness costs. The RNA interference (RNAi) pathway is the major antiviral response vectors use to restrict arbovirus infections. We know this because depleting RNAi gene products profoundly impacts arbovirus replication, the antiviral RNAi pathway genes undergo positive, diversifying selection and arboviruses have evolved strategies to evade the vector's RNAi responses. The vector's RNAi defense and arbovirus countermeasures lead to an arms race that prevents potential virus-induced fitness costs yet maintains arbovirus infections needed for transmission. This review will discuss the latest findings in RNAi-arbovirus interactions in the model insect (Drosophila melanogaster) and in specific mosquito vectors.

Genetic and Phylogenetic Characterization of Tataguine and Witwatersrand Viruses and Other Orthobunyaviruses of the Anopheles A, Capim, Guamá, Koongol, Mapputta, Tete, and Turlock Serogroups

The family Bunyaviridae has more than 530 members that are distributed among five genera or remain to be classified. The genus Orthobunyavirus is the most diverse bunyaviral genus with more than 220 viruses that have been assigned to more than 18 serogroups based on serological cross-reactions and limited molecular-biological characterization. Sequence information for all three orthobunyaviral genome segments is only available for viruses belonging to the Bunyamwera, Bwamba/Pongola, California encephalitis, Gamboa, Group C, Mapputta, Nyando, and Simbu serogroups. Here we present coding-complete sequences for all three genome segments of 15 orthobunyaviruses belonging to the Anopheles A, Capim, Guamá, Kongool, Tete, and Turlock serogroups, and of two unclassified bunyaviruses previously not known to be orthobunyaviruses (Tataguine and Witwatersrand viruses). Using those sequence data, we established the most comprehensive phylogeny of the Orthobunyavirus genus to date, now covering 15 serogroups. Our results emphasize the high genetic diversity of orthobunyaviruses and reveal that the presence of the small nonstructural protein (NSs)-encoding open reading frame is not as common in orthobunyavirus genomes as previously thought.

The Discovery, Distribution, and Evolution of Viruses Associated with Drosophila melanogaster

Drosophila melanogaster is a valuable invertebrate model for viral infection and antiviral immunity, and is a focus for studies of insect-virus coevolution. Here we use a metagenomic approach to identify more than 20 previously undetected RNA viruses and a DNA virus associated with wild D. melanogaster. These viruses not only include distant relatives of known insect pathogens but also novel groups of insect-infecting viruses. By sequencing virus-derived small RNAs, we show that the viruses represent active infections of Drosophila. We find that the RNA viruses differ in the number and properties of their small RNAs, and we detect both siRNAs and a novel miRNA from the DNA virus. Analysis of small RNAs also allows us to identify putative viral sequences that lack detectable sequence similarity to known viruses. By surveying >2,000 individually collected wild adult Drosophila we show that more than 30% of D. melanogaster carry a detectable virus, and more than 6% carry multiple viruses. However, despite a high prevalence of the Wolbachia endosymbiont--which is known to be protective against virus infections in Drosophila--we were unable to detect any relationship between the presence of Wolbachia and the presence of any virus. Using publicly available RNA-seq datasets, we show that the community of viruses in Drosophila laboratories is very different from that seen in the wild, but that some of the newly discovered viruses are nevertheless widespread in laboratory lines and are ubiquitous in cell culture. By sequencing viruses from individual wild-collected flies we show that some viruses are shared between D. melanogaster and D. simulans. Our results provide an essential evolutionary and ecological context for host-virus interaction in Drosophila, and the newly reported viral sequences will help develop D. melanogaster further as a model for molecular and evolutionary virus research.

Evolutionary and phenotypic analysis of live virus isolates suggests arthropod origin of a pathogenic RNA virus family

The evolutionary origins of arboviruses are unknown because their typical dual host tropism is paraphyletic within viral families. Here we studied one of the most diversified and medically relevant RNA virus families, the Bunyaviridae, in which four of five established genera are transmitted by arthropods. We define two cardinally novel bunyavirus groups based on live isolation of 26 viral strains from mosquitoes (Jonchet virus [JONV], eight strains; Ferak virus [FERV], 18 strains). Both viruses were incapable of replicating at vertebrate-typical temperatures but replicated efficiently in insect cells. Replication involved formation of virion-sense RNA (vRNA) and mRNA, including cap-snatching activity. SDS/PAGE, mass spectrometry, and Edman degradation identified translation products corresponding to virion-associated RNA-dependent RNA polymerase protein (RdRp), glycoprotein precursor protein, glycoproteins Gn and Gc, as well as putative nonstructural proteins NSs and NSm. Distinct virion morphologies suggested ancient evolutionary divergence, with bunyavirus-typical morphology for FERV (spheres of 60-120 nm) as opposed to an unusual bimorphology for JONV (tubular virions of 60 × 600 nm and spheres of 80 nm). Both viruses were genetically equidistant from all other bunyaviruses, showing <15% amino acid identity in the RdRp palm domain. Both had different and unique conserved genome termini, as in separate bunyavirus genera. JONV and FERV define two novel sister taxons to the superclade of orthobunyaviruses, tospoviruses, and hantaviruses. Phylogenetic ancestral state reconstruction with probabilistic hypothesis testing suggested ancestral associations with arthropods at deep nodes throughout the bunyavirus tree. Our findings suggest an arthropod origin of bunyaviruses.

Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses

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.

DOI:http://dx.doi.org/10.7554/eLife.05378.001

Many illnesses, including influenza, hemorrhagic fever, and rabies, are caused by a group of viruses called negative-sense RNA viruses. The genetic information—or genome—of these viruses is encoded in strands of RNA that must be copied before they can be translated into the proteins needed to build new viruses. It is currently known that there are at least eight different families of these viruses, which have a wide range of shapes and sizes and arrange their RNA in different ways.

Insects, spiders, and other arthropods carry many different RNA viruses. Many of these viruses have not previously been studied, and those that have been studied so far are mainly those that cause diseases in humans and other vertebrates. Researchers therefore only know a limited amount about the diversity of the negative-sense RNA viruses that arthropods harbor and how these viruses evolved. Studying how viruses evolve helps scientists to understand what makes some viruses deadly and others harmless and can also help develop treatments or vaccines for the diseases caused by the viruses.

Li, Shi, Tian, Lin, Kang et al. collected 70 species of insects, spiders, centipedes, and other arthropods in China and sequenced all the negative-sense RNA viruses in the creatures. This revealed an enormous number of negative-sense RNA viruses, including 112 new viruses. Many of the newly discovered arthropod viruses appear to be the ancestors of disease-causing viruses, including influenza viruses and the filoviruses—the group that includes the Ebola virus. Indeed, it appears that arthropods host many—if not all—of the negative-sense RNA viruses that cause disease in vertebrates and plants.

While documenting the new RNA viruses and how they are related to each other, Li et al. found many different genome structures. Some genomes were segmented, which may play an important role in evolution as segments can be easily swapped to create new genetic combinations. Non-segmented and circular genomes were also found. This genetic diversity suggests that arthropods are likely to have played a key role in the evolution of new viruses by acting as a site where many different viruses can interact and exchange genetic information.

DOI:http://dx.doi.org/10.7554/eLife.05378.002

Replication of boid inclusion body disease-associated arenaviruses is temperature sensitive in both boid and mammalian cells

Boid inclusion body disease (BIDB) is a fatal disease of boid snakes, the etiology of which has only recently been revealed following the identification of several novel arenaviruses in diseased snakes. BIBD-associated arenaviruses (BIBDAV) are genetically divergent from the classical Old and New World arenaviruses and also differ substantially from each other. Even though there is convincing evidence that BIBDAV are indeed the etiological agent of BIBD, the BIBDAV reservoir hosts--if any exist besides boid snakes themselves--are not yet known. In this report, we use University of Helsinki virus (UHV; a virus that we isolated from a Boa constrictor with BIBD) to show that BIBDAV can also replicate effectively in mammalian cells, including human cells, provided they are cultured at 30°C. The infection induces the formation of cytoplasmic inclusion bodies (IB), comprised mainly of viral nucleoprotein (NP), similar to those observed in BIBD and in boid cell cultures. Transferring infected cells from 30°C to 37°C ambient temperature resulted in progressive declines in IB formation and in the amounts of viral NP and RNA, suggesting that BIBDAV growth is limited at 37°C. These observations indirectly indicate that IB formation is linked to viral replication. In addition to mammalian and reptilian cells, UHV infected arthropod (tick) cells when grown at 30°C. Even though our findings suggest that BIBDAV have a high potential to cross the species barrier, their inefficient growth at mammalian body temperatures indicates that the reservoir hosts of BIBDAV are likely species with a lower body temperature, such as snakes.

IMPORTANCE

The newly discovered boid inclusion body disease-associated arenaviruses (BIBDAV) of reptiles have drastically altered the phylogeny of the family Arenavirus. Prior to their discovery, known arenaviruses were considered mainly rodent-borne viruses, with each arenavirus species having its own reservoir host. BIBDAV have so far been demonstrated in captive boid snakes, but their possible reservoir host(s) have not yet been identified. Here we show, using University of Helsinki virus as a model, that these viruses are able to infect mammalian (including human) and arthropod cells. Our results provide in vitro proof of the considerable ability of arenaviruses to cross species barriers. However, our data indicate that BIBDAV growth occurs at 30°C but is inhibited at 37°C, implying that crossing of the species barrier would be hindered by the body temperature of mammalian species.

Evidence that ebolaviruses and cuevaviruses have been diverging from marburgviruses since the Miocene

An understanding of the timescale of evolution is critical for comparative virology but remains elusive for many RNA viruses. Age estimates based on mutation rates can severely underestimate divergences for ancient viral genes that are evolving under strong purifying selection. Paleoviral dating, however, can provide minimum age estimates for ancient divergence, but few orthologous paleoviruses are known within clades of extant viruses. For example, ebolaviruses and marburgviruses are well-studied mammalian pathogens, but their comparative biology is difficult to interpret because the existing estimates of divergence are controversial. Here we provide evidence that paleoviral elements of two genes (ebolavirus-like VP35 and NP) in cricetid rodent genomes originated after the divergence of ebolaviruses and cuevaviruses from marburgviruses. We provide evidence of orthology by identifying common paleoviral insertion sites among the rodent genomes. Our findings indicate that ebolaviruses and cuevaviruses have been diverging from marburgviruses since the early Miocene.