Novel virus discovery and genome reconstruction from field RNA samples reveals highly divergent viruses in dipteran hosts

Tips and tools to obtain and assess mosquito viromes

Due to their vectorial capacity, mosquitoes (Diptera: Culicidae) receive special attention from health authorities and entomologists. These cosmopolitan insects are responsible for the transmission of many viral diseases, such as dengue and yellow fever, causing huge impacts on human health and justifying the intensification of research focused on mosquito-borne diseases. In this context, the study of the virome of mosquitoes can contribute to anticipate the emergence and/or the reemergence of infectious diseases. The assessment of mosquito viromes also contributes to the surveillance of a wide variety of viruses found in these insects, allowing the early detection of pathogens with public health importance. However, the study of mosquito viromes can be challenging due to the number and complexities of steps involved in this type of research. Therefore, this article aims to describe, in a straightforward and simplified way, the steps necessary for obtention and assessment of mosquito viromes. In brief, this article explores: the capture and preservation of specimens; sampling strategies; treatment of samples before DNA/RNA extraction; extraction methodologies; enrichment and purification processes; sequencing choices; and bioinformatics analysis.

Novel strains of Culex flavivirus and Hubei chryso-like virus 1 from the Anopheles mosquito in western Kenya

Surveillance of mosquito vectors is critical for early detection, prevention and control of vector borne diseases. In this study we used advanced molecular tools, such as DNA barcoding in combination with novel sequencing technologies to discover new and already known viruses in genetically identified mosquito species. Mosquitoes were captured using BG sentinel traps in Western Kenya during May and July 2019, and homogenized individually before pooled into groups of ten mosquitoes. The pools and individual samples were then used for molecular analysis and to infect cell cultures. Of a total of fifty-four (54) 10-pools, thirteen (13) showed cytopathic effect (CPE) on VeroB4 cells, eighteen (18) showed CPE on C6/36 cells. Eight (8) 10-pools out of the 31 CPE positive pools showed CPE on both VeroB4 and C6/36 cells. When using reverse transcriptase polymerase chain reaction (RT-PCR), Sanger sequencing and Twist Comprehensive Viral Research Panel (CVRP) (Twist Biosciences), all pools were found negative by RT-PCR when using genus specific primers targeting alphaviruses, orthobunyaviruses and virus specific primers towards o'nyong-nyong virus, chikungunya virus and Sindbis virus (previously reported to circulate in the region). Interestingly, five pools were RT-PCR positive for flavivirus. Two of the RT-PCR positive pools showed CPE on both VeroB4 and C6/36 cells, two pools showed CPE on C6/36 cells alone and one pool on VeroB4 cells only. Fifty individual mosquito homogenates from the five RT-PCR positive 10-pools were analyzed further for flavivirus RNA. Of these, 19 out of the 50 individual mosquito homogenates indicated the presence of flavivirus RNA. Barcoding of the flavivirus positive mosquitoes revealed the mosquito species as Aedes aegypti (1), Mansonia uniformis (6), Anopheles spp (3), Culex pipiens (5), Culex spp (1), Coquilletidia metallica (2) and Culex quinquefasciatus (1). Of the 19 flavivirus positive individual mosquitoes, five (5) virus positive homogenates were sequenced. Genome sequences of two viruses were completed. One was identified as the single-stranded RNA Culex flavivirus and the other as the double-stranded RNA Hubei chryso-like virus 1. Both viruses were found in the same Anopheles spp. homogenate extracted from a sample that showed CPE on both VeroB4 and C6/36 cells. The detection of both viruses in a single mosquito homogenate indicated coinfection. Phylogenetic analyses suggested that the Culex flavivirus sequence detected was closely related to a Culex flavivirus isolated from Uganda in 2008. All four Hubei chryso-like virus 1 segments clusters closely to Hubei chryso-like virus 1 strains isolated in Australia, China and USA. Two novel strains of insect-specific viruses in Anopheles mosquitoes were detected and characterized.

The minus strand of positive-sense RNA viruses encodes small proteins

It is widely accepted that the minus strands of positive single-strand RNA (+ssRNA) viruses function as replication templates only. Gong et al. revealed that the minus strand of two unrelated +ssRNA viruses encodes proteins. This textbook-changing discovery calls for the reconsideration of the molecular mechanisms underlying the infection cycle of +ssRNA viruses.

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.

Identification of RNA Virus-Derived RdRp Sequences in Publicly Available Transcriptomic Data Sets

RNA viruses are abundant and highly diverse and infect all or most eukaryotic organisms. However, only a tiny fraction of the number and diversity of RNA virus species have been catalogued. To cost-effectively expand the diversity of known RNA virus sequences, we mined publicly available transcriptomic data sets. We developed 77 family-level Hidden Markov Model profiles for the viral RNA-dependent RNA polymerase (RdRp)-the only universal "hallmark" gene of RNA viruses. By using these to search the National Center for Biotechnology Information Transcriptome Shotgun Assembly database, we identified 5,867 contigs encoding RNA virus RdRps or fragments thereof and analyzed their diversity, taxonomic classification, phylogeny, and host associations. Our study expands the known diversity of RNA viruses, and the 77 curated RdRp Profile Hidden Markov Models provide a useful resource for the virus discovery community.

Composition and global distribution of the mosquito virome - A comprehensive database of insect-specific viruses

Mosquitoes are vectors for emerging and re-emerging infectious viral diseases of humans, livestock and other animals. In addition to these arthropod-borne (arbo)viruses, mosquitoes are host to an array of insect-specific viruses, collectively referred to as the mosquito virome. Mapping the mosquito virome and understanding if and how its composition modulates arbovirus transmission is critical to understand arboviral disease emergence and outbreak dynamics. In recent years, next-generation sequencing as well as PCR and culture-based methods have been extensively used to identify mosquito-associated viruses, providing insights into virus ecology and evolution. Until now, the large amount of mosquito virome data, specifically those acquired by metagenomic sequencing, has not been comprehensively integrated. We have constructed a searchable database of insect-specific viruses associated with vector mosquitoes from 175 studies, published between October 2000 and February 2022. We identify the most frequently detected and widespread viruses of the Culex, Aedes and Anopheles mosquito genera and report their global distribution. In addition, we highlight the challenges of extracting and integrating published virome data and we propose that a standardized reporting format will facilitate data interpretation and re-use by other scientists. We expect our comprehensive database, summarizing mosquito virome data collected over 20 years, to be a useful resource for future studies.

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.

An in silico analysis revealed a novel evolutionary lineage of putative mitoviruses

Mitoviruses (family Mitoviridae) are small capsid-less RNA viruses that replicate in the mitochondria of fungi and plants. However, to date, the only authentic animal mitovirus infecting an insect was identified as Lutzomyia longipalpis mitovirus 1 (LulMV1). Public databases of transcriptomic studies from several animals may be a good source for identifying the often missed mitoviruses. Consequently, a search of mitovirus-like transcripts at the NCBI transcriptome shotgun assembly (TSA) library, and a search for the mitoviruses previously recorded at the NCBI non-redundant (nr) protein sequences library, were performed in order to identify new mitovirus-like sequences associated with animals. In total, 10 new putative mitoviruses were identified in the TSA database and 8 putative mitoviruses in the nr protein database. To our knowledge, these results represent the first evidence of putative mitoviruses associated with poriferan, cnidarians, echinoderms, crustaceans, myriapods and arachnids. According to different phylogenetic inferences using the maximum likelihood method, these 18 putative mitoviruses form a robust monophyletic lineage with LulMV1, the only known animal-infecting mitovirus. These findings based on in silico procedures provide strong evidence for the existence of a clade of putative mitoviruses associated with animals, which has been provisionally named 'kvinmitovirus'.

Isolation of a novel rhabdovirus and detection of multiple novel viral sequences in Culex species mosquitoes in the United States

To increase our understanding of the diversity of the mosquito virome, 6956 mosquitoes of five species (Culex erraticus, Culex pipiens, Culex restuans, Culex tarsalis, and Culex territans) collected in Iowa in the United States in 2017 and 2020 were assayed for novel viruses by performing polyethylene glycol precipitation, virus isolation in cell culture, and unbiased high-throughput sequencing. A novel virus, provisionally named "Walnut Creek virus", was isolated from Cx. tarsalis, and its genomic sequence and organization are characteristic of viruses in the genus Hapavirus (family Rhabdoviridae). Replication of Walnut Creek virus occurred in avian, mammalian, and mosquito, but not tick, cell lines. A novel virus was also isolated from Cx. restuans, and partial genome sequencing revealed that it is distantly related to an unclassified virus of the genus Phytoreovirus (family Sedoreoviridae). Two recognized viruses were also isolated: Culex Y virus (family Birnaviridae) and Houston virus (family Mesoniviridae). We also identified sequences of eight novel viruses from six families (Amalgaviridae, Birnaviridae, Partitiviridae, Sedoreoviridae, Tombusviridae, and Totiviridae), two viruses that do not belong to any established families, and many previously recognized viruses. In summary, we provide evidence of multiple novel and recognized viruses in Culex spp. mosquitoes in the United States.

Interactions of host miRNAs in the flavivirus 3´UTR genome: From bioinformatics predictions to practical approaches

The genus Flavivirus of the Flaviviridae family includes important viruses, such as Dengue, Zika, West Nile, Japanese encephalitis, Murray Valley encephalitis, tick-borne encephalitis, Yellow fever, Saint Louis encephalitis, and Usutu viruses. They are transmitted by mosquitoes or ticks, and they can infect humans, causing fever, encephalitis, or haemorrhagic fever. The treatment resources for these diseases and the number of vaccines available are limited. It has been discovered that eukaryotic cells synthesize small RNA molecules that can bind specifically to sequences present in messenger RNAs to inhibit the translation process, thus regulating gene expression. These small RNAs have been named microRNAs, and they have an important impact on viral infections. In this review, we compiled the available information on miRNAs that can interact with the 3’ untranslated region (3’UTR) of the flavivirus genome, a conserved region that is important for viral replication and translation.

Potential origins of fish toti-like viruses in invertebrates

The virus family Totiviridae had originally been considered to include only viruses which infected fungal and protist hosts, but since 2006 a growing number of viruses found in invertebrates and fish have been shown to cluster phylogenetically within this family. These Totiviridae-like, or toti-like, viruses do not appear to belong within any existing genera of Totiviridae, and whilst a number of new genus names have been suggested, none has yet been universally accepted. Within this growing number of toti-like viruses from animal hosts, there exists emerging viral threats particularly to aquaculture, namely Infectious myonecrosis virus in whiteleg shrimp and Piscine myocarditis virus (PMCV) in Atlantic salmon (Salmo salar). PMCV in particular continues to be an issue in salmon aquaculture as a number of questions remain unanswered about how the virus is transmitted and the route of entry into host fish. Using a phylogenetic approach, this study shows how PMCV and the other fish toti-like viruses probably have deeper origins in an arthropod host. Based on this, it is hypothesized that sea lice could be acting as a vector for PMCV, as seen with other RNA viruses in Atlantic salmon aquaculture and in the toti-like Cucurbit yellows-associated virus which is spread by the greenhouse whitefly Trialeurodes vaporariorum.

Complete genome analysis of a novel chuvirus from a southern green stink bug (Nezara viridula)

A novel chuvirus from a southern green stink bug (Nezara viridula) was identified by RNA sequencing in this study and was tentatively named "Ningbo southern green stink bug chuvirus 1" (NBSGSBV-1). The complete genome sequence of NBSGSBV-1 consists of 11,375 nucleotides, and the genome was found to be circular by 'around-the-genome' reverse transcription polymerase chain reaction (RT-PCR) and Sanger sequencing. Three open reading frames (ORFs) were predicted in the NBSGSBV-1 genome, encoding a large polymerase protein (L protein), a glycoprotein (G protein), and a nucleocapsid protein (N protein). A phylogenetic tree was constructed based on all of the currently available RNA-dependent RNA polymerase amino acid sequences of viruses of the family Chuviridae, and NBSGSBV-1 was found to cluster together with Sanya chuvirus 2 and Hubei odonate virus 11, indicating that NBSGSBV-1 might belong to the genus Odonatavirus. Five conserved sites were identified in the L proteins of NBSGSBV-1 and other chuviruses. The abundance and characteristics of the NBSGSBV-1-derived small interfering RNAs suggested that NBSGSBV-1 actively replicates in the host insect. To the best of our knowledge, this is the first report of a chuvirus identified in a member of the insect family Pentatomidae. The discovery and characterization of NBSGSBV-1 will help us to understand the diversity of chuviruses in insects.

Fusarium Graminearum Virus-1 Strain FgV1-SD4 Infection Eliminates Mycotoxin Deoxynivalenol Synthesis by Fusarium graminearum in FHB

Deoxynivalenol (DON) toxin production during the infection of F. graminearum in small grain crops is one of the most harmful virulence factors associated with economic losses. Metatranscriptome sequencing and RT-qPCR traced back that the only mycovirus infecting an F. graminearum isolate, designated as Fg-4-2, was a novel strain of Fusarium graminearum virus 1 (FgV1), designated as FgV1-SD4. The isolate Fg-4-2 showed significantly reduced virulence against wheat compared to the virus-free culture, designated as isolate Fg-4-1, which was obtained by deep freezing and single conidial germination. Notably, no DON accumulation was detected in the harvested wheat seeds infected by Fg-4-2, whereas ~18 ppm DON was detected in seeds infected by Fg-4-1. Comparison of the genome sequence of FgV1-SD4 with other identified strains of FgV1, i.e., FgV1-DK21 and FgV1-ch, indicates mutations on ORF-2 and the 3′-UTR in the genome that might be associated with hypovirulence. This mycovirus strain alone and specific genetic components of FgV1-SD4 can be further optimized to be developed as a biocontrol agent to reduce Fusarium head blight and to lower the DON accumulation levels in small grain crops due to this fungal disease.

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.

Metagenomic sequencing characterizes a wide diversity of viruses in field mosquito samples in Nigeria

Mosquito vectors are a tremendous public health threat. One in six diseases worldwide is vector-borne transmitted mainly by mosquitoes. In the last couple of years, there have been active Yellow fever virus (YFV) outbreaks in many settings in Nigeria, and nationwide, entomological surveillance has been a significant effort geared towards understanding these outbreaks. In this study, we used a metagenomic sequencing approach to characterize viruses present in vector samples collected during various outbreaks of Yellow fever (YF) in Nigeria between 2017 and 2020. Mosquito samples were grouped into pools of 1 to 50 mosquitoes, each based on species, sex and location. Twenty-five pools of Aedes spp and one pool of Anopheles spp collected from nine states were sequenced and metagenomic analysis was carried out. We identified a wide diversity of viruses belonging to various families in this sample set. Seven different viruses detected included: Fako virus, Phasi Charoen-like virus, Verdadero virus, Chaq like-virus, Aedes aegypti totivirus, cell fusing agent virus and Tesano Aedes virus. Although there are no reports of these viruses being pathogenic, they are an understudied group in the same families and closely related to known pathogenic arboviruses. Our study highlights the power of next generation sequencing in identifying Insect specific viruses (ISVs), and provide insight into mosquito vectors virome in Nigeria.

Dissecting the Species-Specific Virome in Culicoides of Thrace

Biting midges (Culicoides) are vectors of arboviruses of both veterinary and medical importance. The surge of emerging and reemerging vector-borne diseases and their expansion in geographical areas affected by climate change has increased the importance of understanding their capacity to contribute to novel and emerging infectious diseases. The study of Culicoides virome is the first step in the assessment of this potential. In this study, we analyzed the RNA virome of 10 Culicoides species within the geographical area of Thrace in the southeastern part of Europe, a crossing point between Asia and Europe and important path of various arboviruses, utilizing the Ion Torrent next-generation sequencing (NGS) platform and a custom bioinformatics pipeline based on TRINITY assembler and alignment algorithms. The analysis of the RNA virome of 10 Culicoides species resulted in the identification of the genomic signatures of 14 novel RNA viruses, including three fully assembled viruses and four segmented viruses with at least one segment fully assembled, most of which were significantly divergent from previously identified related viruses from the Solemoviridae, Phasmaviridae, Phenuiviridae, Reoviridae, Chuviridae, Partitiviridae, Orthomyxoviridae, Rhabdoviridae, and Flaviviridae families. Each Culicoides species carried a species-specific set of viruses, some of which are related to viruses from other insect vectors in the same area, contributing to the idea of a virus-carrier web within the ecosystem. The identified viruses not only expand our current knowledge on the virome of Culicoides but also set the basis of the genetic diversity of such viruses in the area of southeastern Europe. Furthermore, our study highlights that such metagenomic approaches should include as many species as possible of the local virus-carrier web that interact and share the virome of a geographical area.

Updated occurrence and bionomics of potential malaria vectors in Europe: a systematic review (2000-2021)

Despite the eradication of malaria across most European countries in the 1960s and 1970s, the anopheline vectors are still present. Most of the malaria cases that have been reported in Europe up to the present time have been infections acquired in endemic areas by travelers. However, the possibility of acquiring malaria by locally infected mosquitoes has been poorly investigated in Europe, despite autochthonous malaria cases having been occasionally reported in several European countries. Here we present an update on the occurrence of potential malaria vector species in Europe. Adopting a systematic review approach, we selected 288 papers published between 2000 and 2021 for inclusion in the review based on retrieval of accurate information on the following Anopheles species: An. atroparvus, An. hyrcanus sensu lato (s.l.), An. labranchiae, An. maculipennis sensu stricto (s.s.), An. messeae/daciae, An. sacharovi, An. superpictus and An. plumbeus. The distribution of these potential vector species across Europe is critically reviewed in relation to areas of major presence and principal bionomic features, including vector competence to Plasmodium. Additional information, such as geographical details, sampling approaches and species identification methods, are also reported. We compare the information on each species extracted from the most recent studies to comparable information reported from studies published in the early 2000s, with particular reference to the role of each species in malaria transmission before eradication. The picture that emerges from this review is that potential vector species are still widespread in Europe, with the largest diversity in the Mediterranean area, Italy in particular. Despite information on their vectorial capacity being fragmentary, the information retrieved suggests a re-definition of the relative importance of potential vector species, indicating An. hyrcanus s.l., An. labranchiae, An. plumbeus and An. sacharovi as potential vectors of higher importance, while An. messeae/daciae and An. maculipennis s.s. can be considered to be moderately important species. In contrast, An. atroparvus and An. superpictus should be considered as vectors of lower importance, particularly in relation to their low anthropophily. The presence of gaps in current knowledge of vectorial systems in Europe becomes evident in this review, not only in terms of vector competence but also in the definition of sampling approaches, highlighting the need for further research to adopt the appropriate surveillance system for each species.

Analysis of avian Usutu virus infections in Germany from 2011 to 2018 with focus on dsRNA detection to demonstrate viral infections

Usutu virus (USUV) is a zoonotic arbovirus causing avian mass mortalities. The first outbreak in North-Western Germany occurred in 2018. This retrospective analysis focused on combining virological and pathological findings in birds and immunohistochemistry. 25 common blackbirds, one great grey owl, and one kingfisher collected from 2011 to 2018 and positive for USUV by qRT-PCR were investigated. Macroscopically, most USUV infected birds showed splenomegaly and hepatomegaly. Histopathological lesions included necrosis and lymphohistiocytic inflammation within spleen, Bursa fabricii, liver, heart, brain, lung and intestine. Immunohistochemistry revealed USUV antigen positive cells in heart, spleen, pancreas, lung, brain, proventriculus/gizzard, Bursa fabricii, kidney, intestine, skeletal muscle, and liver. Analysis of viral genome allocated the virus to Europe 3 or Africa 2 lineage. This study investigated whether immunohistochemical detection of double-stranded ribonucleic acid (dsRNA) serves as an alternative tool to detect viral intermediates. Tissue samples of six animals with confirmed USUV infection by qRT-PCR but lacking viral antigen in liver and spleen, were further examined immunohistochemically. Two animals exhibited a positive signal for dsRNA. This could indicate either an early state of infection without sufficient formation of virus translation products, occurrence of another concurrent virus infection or endogenous dsRNA not related to infectious pathogens and should be investigated in more detail in future studies.

Novel Picornavirus Detected in Wild Deer: Identification, Genomic Characterisation, and Prevalence in Australia

The use of high-throughput sequencing has facilitated virus discovery in wild animals and helped determine their potential threat to humans and other animals. We report the complete genome sequence of a novel picornavirus identified by next-generation sequencing in faeces from Australian fallow deer. Genomic analysis revealed that this virus possesses a typical picornavirus-like genomic organisation of 7554 nt with a single open reading frame (ORF) encoding a polyprotein of 2225 amino acids. Based on the amino acid identity comparison and phylogenetic analysis of the P1, 2C, 3CD, and VP1 regions, this novel picornavirus was closely related to but distinct from known bopiviruses detected to date. This finding suggests that deer/bopivirus could belong to a novel species within the genus Bopivirus, tentatively designated as "Bopivirus C". Epidemiological investigation of 91 deer (71 fallow, 14 sambar and 6 red deer) and 23 cattle faecal samples showed that six fallow deer and one red deer (overall prevalence 7.7%, 95% confidence interval [CI] 3.8-15.0%) tested positive, but deer/bopivirus was undetectable in sambar deer and cattle. In addition, phylogenetic and sequence analyses indicate that the same genotype is circulating in south-eastern Australia. To our knowledge, this study reports for the first time a deer-origin bopivirus and the presence of a member of genus Bopivirus in Australia. Further epidemiological and molecular studies are needed to investigate the geographic distribution and pathogenic potential of this novel Bopivirus species in other domestic and wild animal species.

Identification of potential new mosquito-associated viruses of adult Aedes aegypti mosquitoes from Tocantins state, Brazil

Medically important arboviruses such as dengue virus (DENV), Zika virus (ZIKV), and chikungunya virus (CHIKV) are primarily transmitted by the globally distributed mosquito Aedes aegypti. Increasing evidence suggests that the transmission of some viruses can be influenced by mosquito-specific and mosquito-borne viruses. Advancements in high-throughput sequencing (HTS) and bioinformatics have expanded our knowledge on the richness of viruses harbored by mosquitoes. HTS was used to characterize the presence of virus sequences in wild-caught adult Ae. aegypti from Tocantins (TO) state, Brazil. Samples of mosquitoes were collected in four cities of Tocantins state and submitted to RNA isolation, followed by sequencing at an Illumina HiSeq platform. Our results showed initially by Krona the presence of 3% of the sequenced reads belonging to the viral database. After further analysis, the virus sequences were found to have homology to two viral families found in insects Phenuiviridae and Metaviridae. Three possible viral strains including putative new viruses were detected and named Phasi Charoen-like phasivirus isolate To-1 (PCLV To-1), Aedes aegypti To virus 1 (AAToV1), and Aedes aegypti To virus 2 (AAToV2). The results presented in this work contribute to the growing knowledge about the diversity of viruses in mosquitoes and might be useful for future studies on the interaction between insect-specific viruses and arboviruses.

Uncovering the Worldwide Diversity and Evolution of the Virome of the Mosquitoes Aedes aegypti and Aedes albopictus

Aedes aegypti, the yellow fever mosquito, and Aedes albopictus, the Asian tiger mosquito, are the most significant vectors of dengue, Zika, and Chikungunya viruses globally. Studies examining host factors that control arbovirus transmission demonstrate that insect-specific viruses (ISVs) can modulate mosquitoes’ susceptibility to arbovirus infection in both in vivo and in vitro co-infection models. While research is ongoing to implicate individual ISVs as proviral or antiviral factors, we have a limited understanding of the composition and diversity of the Aedes virome. To address this gap, we used a meta-analysis approach to uncover virome diversity by analysing ~3000 available RNA sequencing libraries representing a worldwide geographic range for both mosquitoes. We identified ten novel viruses and previously characterised viruses, including mononegaviruses, orthomyxoviruses, negeviruses, and a novel bi-segmented negev-like group. Phylogenetic analysis suggests close relatedness to mosquito viruses implying likely insect host range except for one arbovirus, the multi-segmented Jingmen tick virus (Flaviviridae) in an Italian colony of Ae. albopictus. Individual mosquito transcriptomes revealed remarkable inter-host variation of ISVs within individuals from the same colony and heterogeneity between different laboratory strains. Additionally, we identified striking virus diversity in Wolbachia infected Aedes cell lines. This study expands our understanding of the virome of these important vectors. It provides a resource for further assessing the ecology, evolution, and interaction of ISVs with their mosquito hosts and the arboviruses they transmit.

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.

Discovery and Characterization of Actively Replicating DNA and Retro-Transcribing Viruses in Lower Vertebrate Hosts Based on RNA Sequencing

In a previous study, a metatranscriptomics survey of RNA viruses in several important lower vertebrate host groups revealed huge viral diversity, transforming the understanding of the evolution of vertebrate-associated RNA virus groups. However, the diversity of the DNA and retro-transcribing viruses in these host groups was left uncharacterized. Given that RNA sequencing is capable of revealing viruses undergoing active transcription and replication, we collected previously generated datasets associated with lower vertebrate hosts, and searched them for DNA and retro-transcribing viruses. Our results revealed the complete genome, or “core gene sets”, of 18 vertebrate-associated DNA and retro-transcribing viruses in cartilaginous fishes, ray-finned fishes, and amphibians, many of which had high abundance levels, and some of which showed systemic infections in multiple organs, suggesting active transcription or acute infection within the host. Furthermore, these new findings recharacterized the evolutionary history in the families Hepadnaviridae, Papillomaviridae, and Alloherpesviridae, confirming long-term virus–host codivergence relationships for these virus groups. Collectively, our results revealed reliable and sufficient information within metatranscriptomics sequencing to characterize not only RNA viruses, but also DNA and retro-transcribing viruses, and therefore established a key methodology that will help us to understand the composition and evolution of the total “infectome” within a diverse range of vertebrate hosts.

Molecular epidemiology of mosquito-borne viruses at the China-Myanmar border: discovery of a potential epidemic focus of Japanese encephalitis

BACKGROUND

Mosquito-based arbovirus surveillance can serve as an early warning in evaluating the status of mosquito-borne virus prevalence and thus prevent local outbreaks. Although Tengchong County in Yunnan Province-which borders Myanmar-is abundant and diverse in mosquitoes, very few mosquito-based arbovirus investigations have been conducted in the recent decade. Herein, this study aims to evaluate the presence and the diffusion of mosquito-borne pathogens, currently prevalent in this region.

METHODS

We collected 9486 mosquitoes, representing eight species, with Culex tritaeniorhynchus and Anopheles sinensis as the dominant species, during high mosquito activity seasons (July-October) in Tengchong, in 2018. Samples collected from 342 pools were tested using reverse-transcription PCR to determine the species, distribution, and infection rates of virus and parasite, and further analyze their genotypes, phylogenetic relationships, infection rate, and potential pathogenicity.

RESULTS

Fifteen Japanese encephalitis virus (JEV) strains from Cx. tritaeniorhynchus pools were detected. Seven strains of insect-specific flaviviruses (ISFVs), including two Aedes flavivirus (AeFV) and Yunnan Culex flavivirus strains each, one Culex theileri flavivirus, Yamadai flavivirus (YDFV) and Anopheles-associated flavivirus (AAFV) strains each were detected in Aedes albopictus, Cx. tritaeniorhynchus, Cx. vagans, Cx. pseudovihnui, and An. sinensis pools, respectively. The whole-genome was successfully amplified in one strain of JEV and AeFV each. Phylogenetic analysis using the E gene placed all the newly detected JEV strains into the GI-b genotype. They showed highly nucleotide identities, and were most closely related to the strain detected in Tengchong in 2010. The comparison of the E protein of JEV strains and vaccine-derived strain, showed six amino residue differences. The bias-corrected maximum likelihood estimation values (and 95% confidence interval) for JEV in Cx. tritaeniorhynchus collected in Tengchong in 2018 were 2.4 (1.4-3.9).

CONCLUSIONS

A potential Japanese encephalitis epidemic focus with the abundance of host mosquitoes and high JEV infection rate was observed in Tengchong. In addition, at least five species of ISFVs co-circulate in this area. This study highlights the importance of widespread and sustained mosquito-based arbovirus surveillance in local areas to prevent the transmission of JEV, and other emerging/re-emerging mosquito-borne pathogens.

Discovery of a Novel Flavivirus (Flaviviridae) From the Horse Fly, Tabanus rufidens (Diptera: Tabanidae): The Possible Coevolutionary Relationships Between the Classical Insect-Specific Flaviviruses and Host Dipteran Insects

Tabanid flies (Tabanidae: Diptera) are common hematophagous insects known to transmit some pathogens mechanically or biologically to animals; they are widely distributed throughout the world. However, no tabanid-borne viruses, except mechanically transmitted viruses, have been reported to date. In this study, we conducted RNA virome analysis of several human-biting tabanid species in Japan, to discover and characterize viruses associated with tabanids. A novel flavivirus was encountered during the study in the Japanese horse fly, Tabanus rufidens (Bigot, 1887). The virus was detected only in T. rufidens, but not in other tabanid species, and as such was designated Tabanus rufidens flavivirus (TrFV). TrFV could not be isolated using a mammalian cell line and showed a closer phylogenetic relationship to the classical insect-specific flaviviruses (cISFs) rather than the vertebrate-infecting flaviviruses (VIFs), suggesting that it is a novel member of the cISFs. The first discovery of a cISF from Brachycera provides new insight into the evolutionary history and dynamics of flaviviruses.

Insect-Specific Viruses: An overview and their relationship to arboviruses of concern to humans and animals

The group of Insect-specific viruses (ISVs) includes viruses apparently restricted to insects based on their inability to replicate in the vertebrates. Increasing numbers of ISVs have been discovered and characterized representing a diverse number of viral families. However, most studies have focused on those ISVs belonging to the family Flaviviridae, which highlights the importance of ISV study from other viral families, which allow a better understanding for the mechanisms of transmission and evolution used for this diverse group of viruses. Some ISVs have shown the potential to modulate arboviruses replication and vector competence of mosquitoes. Based on this, ISVs may be used as an alternative tool for biological control, development of vaccines, and diagnostic platforms for arboviruses. In this review, we provide an update of the general characteristics of ISVs and their interaction with arboviruses that infect vertebrates.

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.

A Unique Relative of Rotifer Birnavirus Isolated from Australian Mosquitoes

The family Birnaviridae are a group of non-enveloped double-stranded RNA viruses which infect poultry, aquatic animals and insects. This family includes agriculturally important pathogens of poultry and fish. Recently, next-generation sequencing technologies have identified closely related birnaviruses in Culex, Aedes and Anopheles mosquitoes. Using a broad-spectrum system based on detection of long double-stranded RNA, we have discovered and isolated a birnavirus from Aedes notoscriptus mosquitoes collected in northern New South Wales, Australia. Phylogenetic analysis of Aedes birnavirus (ABV) showed that it is related to Rotifer birnavirus, a pathogen of microscopic aquatic animals. In vitro cell infection assays revealed that while ABV can replicate in Aedes-derived cell lines, the virus does not replicate in vertebrate cells and displays only limited replication in Culex- and Anopheles-derived cells. A combination of SDS-PAGE and mass spectrometry analysis suggested that the ABV capsid precursor protein (pVP2) is larger than that of other birnaviruses and is partially resistant to trypsin digestion. Reactivity patterns of ABV-specific polyclonal and monoclonal antibodies indicate that the neutralizing epitopes of ABV are SDS sensitive. Our characterization shows that ABV displays a number of properties making it a unique member of the Birnaviridae and represents the first birnavirus to be isolated from Australian mosquitoes.

Understanding the Mechanisms Underlying Host Restriction of Insect-Specific Viruses

Arthropod-borne viruses contribute significantly to global mortality and morbidity in humans and animals. These viruses are mainly transmitted between susceptible vertebrate hosts by hematophagous arthropod vectors, especially mosquitoes. Recently, there has been substantial attention for a novel group of viruses, referred to as insect-specific viruses (ISVs) which are exclusively maintained in mosquito populations. Recent discoveries of novel insect-specific viruses over the past years generated a great interest not only in their potential use as vaccine and diagnostic platforms but also as novel biological control agents due to their ability to modulate arbovirus transmission. While arboviruses infect both vertebrate and invertebrate hosts, the replication of insect-specific viruses is restricted in vertebrates at multiple stages of virus replication. The vertebrate restriction factors include the genetic elements of ISVs (structural and non-structural genes and the untranslated terminal regions), vertebrate host factors (agonists and antagonists), and the temperature-dependent microenvironment. A better understanding of these bottlenecks is thus warranted. In this review, we explore these factors and the complex interplay between ISVs and their hosts contributing to this host restriction phenomenon.

Aedes aegypti from Amazon Basin Harbor High Diversity of Novel Viral Species

Viruses are the most diverse and abundant microorganisms on earth, highly adaptive to a wide range of hosts. Viral diversity within invertebrate hosts has gained notoriety in recent years in public health as several such viruses have been of medical importance. Aedes aegypti serves as a vector for several viruses that have caused epidemics within the last year throughout Brazil; including Dengue, Zika and Chikungunya. This study aimed to identify new viral agents within Aedes aegypti mosquito in a city of the Amazonian region, where it is highly endemic. Metagenomic investigation was performed on 60 mosquito pools and viral RNA sequences present in their microbiota were characterized using genomic and phylogenetic tools. In total, we identified five putative novel virus species related to the Sobemovirus genus, Iflavirus genus and Permutatetraviridae family. These findings indicate a diverse taxonomy of viruses present in the mosquito microbiota of the Amazon, the region with the greatest invertebrate diversity in the world.

Leafhopper Psammotettix alienus hosts chuviruses with different genomic structures

The leafhopper Psammotettix alienus causes serious economic losses by directly sucking the plant sap and by transmitting plant viruses. Here, we obtained the full-genomic sequence of a novel chuvirus, Hancheng leafhopper mivirus (HCLeV), using RNA sequencing (RNA-seq) and rapid amplification of cDNA ends (RACE). The full genome of HCLeV comprises 9,852 nucleotides (nt) and includes two open reading frames (ORFs), which encode nucleocapsid and large polymerase protein, respectively. Its genomic features are similar to that of Tacheng tick virus 5, a member of chuviruses. However, HCLeV had only 36.6-44.4% identities in amino acid sequence for polymerase and 27.1-44.2% identities for nucleocapsid protein with chuviruses and other viruses in the order Mononegavirales and Jingchuvirales. Interestingly, we found this leafhopper can host two chuvirus species with different genomic structures. The discovery of this new virus potentially adds a new species to the family Chuviridae, and its new genomic structural form indicates broader genomic diversity among the chuviruses.

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.

Characterization of Port Bolivar Virus, a Novel Entomobirnavirus (Birnaviridae) Isolated from Mosquitoes Collected in East Texas, USA

This report describes and characterizes a novel entomobirnavirus, designated Port Bolivar virus (PTBV), that was isolated from a pool of Aedes sollicitans mosquitoes collected in a saltwater marsh in East Texas, USA. Full genome sequencing and phylogenetic analyses indicate that PTBV is distinct but genetically related to Drosophila X virus and mosquito X virus, which are assigned to species in the genus Entomobirnavirus, family Birnaviridae. PTBV produced cytopathic effect (CPE) in cultures of mosquito (C6/36) cells, but not in Vero cell cultures. Ultrastructural studies of PTBV in infected C6/36 cells demonstrated unenveloped virus particles about 55 nm in diameter.

New Viral Sequences Identified in the Flavescence Dorée Phytoplasma Vector Scaphoideus titanus

(1) Background: The leafhopper Scaphoideus titanus is the primary vector of Flavescence dorée phytoplasma (FDp) in European vineyards. Flavescence dorée is one of the most severely damaging diseases of Vitis vinifera and, consequently, a major threat to grape and wine production in several European countries. Control measures are compulsory, but they mainly involve large-scale insecticide treatments, with detrimental impacts on the environment. One possible solution is to exploit the largely unexplored genetic diversity of viruses infecting S. titanus as highly specific and environmentally benign tools for biological control. (2) Methods: A metatranscriptomic approach was adopted to identify viruses that may infect individuals caught in the wild in both its native (United States) and invasive (Europe) areas. Reverse transcription PCR was used to confirm their presence in RNA pools and explore their prevalence. (3) Results: We described nine new RNA viruses, including members of “Picorna-Calici”, “Permutotetra”, “Bunya-Arena”, “Reo”, “Partiti-Picobirna”, “Luteo-Sobemo” and “Toti-Chryso” clades. A marked difference in the diversity and abundance of the viral species was observed between the USA population and the European ones. (4) Conclusions: This work represents the first survey to assess the viral community of a phytoplasma insect vector. The possibility to exploit these naturally occurring viruses as specific and targeted biocontrol agents of S. titanus could be the answer to increasing demand for a more sustainable viticulture.

Characterization and Incidence of the First Member of the Genus Mitovirus Identified in the Phytopathogenic Species Fusarium oxysporum

A novel mycovirus named Fusarium oxysporum f. sp. dianthi mitovirus 1 (FodMV1) has been identified infecting a strain of Fusarium oxysporum f. sp. dianthi from Colombia. The genome of FodMV1 is 2313 nt long, and comprises a 172-nt 5’-UTR, a 2025-nt single ORF encoding an RdRp of 675 amino acid residues, and a 113-nt 3´-UTR. Homology BlastX searches identifies FodMV1 as a novel member of the genus Mitovirus in the family Narnaviridae. As the rest of mitoviruses, the genome of FodMV1 presents a high percentage of A+U (58.8%) and contains a number of UGA codons that encode the amino acid tryptophan rather than acting as stop codons as in the universal genetic code. Another common feature with other mitoviruses is that the 5′- and 3′-UTR regions of FodMV1 can be folded into potentially stable stem-loop structures. Result from phylogenetic analysis place FodMV1 in a different clade than the rest of mitoviruses described in other Fusarium spp. Incidence of FodMV1-infections in the collection of F. oxysporum f. sp. dianthi isolates analyzed is relatively high. Of particular interest is the fact that FodMV1 has been detected infecting isolates from two geographical areas as distant as Spain and Colombia.

Viromics on Honey-Baited FTA Cards as a New Tool for the Detection of Circulating Viruses in Mosquitoes

Worldwide, emerging and re-emerging infectious diseases (EIDs) are a major burden on public and animal health. Arthropod vectors, with mosquitoes being the main contributors of global disease, transmit more than 70% of the recognized EIDs. To assess new alternatives for arthropod-borne viral diseases surveillance, and for the detection of new viruses, honey-baited Flinders Technology Associates (FTA) cards were used as sugar bait in mosquito traps during entomological surveys at the Llobregat River Delta (Catalonia, Spain). Next generation sequencing (NGS) metagenomics analysis was applied on honey-baited FTA cards, which had been exposed to field-captured mosquitoes to characterize their associated virome. Arthropod- and plant-infecting viruses governed the virome profile on FTA cards. Twelve near-complete viral genomes were successfully obtained, suggesting good quality preservation of viral RNAs. Mosquito pools linked to the FTA cards were screened for the detection of mosquito-associated viruses by specific RT-PCRs to confirm the presence of these viruses. The circulation of viruses related to Alphamesonivirus, Quaranjavirus and unclassified Bunyavirales was detected in mosquitoes, and phylogenetic analyses revealed their similarities to viruses previously reported in other continents. To the best our knowledge, our findings constitute the first distribution record of these viruses in European mosquitoes and the first hint of insect-specific viruses in mosquitoes' saliva in field conditions, demonstrating the feasibility of this approach to monitor the transmissible fraction of the mosquitoes' virome. In conclusion, this pilot viromics study on honey-baited FTA cards was shown to be a valid approach for the detection of viruses circulating in mosquitoes, thereby setting up an alternative tool for arbovirus surveillance and control programs.

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.

The Detection and Surveillance of Asian Citrus Psyllid (Diaphorina citri)-Associated Viruses in Florida Citrus Groves

The plant pathogenic bacterium Candidatus Liberibacter asiaticus (CLas), the causal agent of the citrus disease Huanglongbing (HLB), and its insect vector, the Asian citrus psyllid (ACP; Diaphorina citri), have been devastating the Florida citrus industry. To restore the competitive production presence of Florida in the worldwide citrus market, effective and sustainable control of HLB and the ACP needs to be identified. As alternatives for resistance-inducing insecticides, viruses are currently being considered for biological control of the ACP. To identify possible biological control candidates, we conducted one of the most comprehensive surveys of natural ACP populations in major citrus production regions spanning 21 counties in Florida. By optimizing PCRs and RT-PCRs, we were able to successfully detect and monitor the prevalence of five previously identified ACP-associated RNA and DNA viruses throughout Florida citrus groves, which include: Diaphorina citri-associated C virus (DcACV), Diaphorina citri flavi-like virus (DcFLV), Diaphorina citri densovirus (DcDNV), Diaphorina citri reovirus (DcRV), and Diaphorina citri picorna-like virus (DcPLV). Adult and nymph ACP populations from 21 of Florida's major citrus-producing counties were collected each month during approximately 18 consecutive months. RNA extracts used for these viral screens were also regionally combined and subjected to High Throughput Sequencing (HTS) to reveal a more comprehensive picture of known and unknown viruses in Florida ACP populations. We discovered that DcACV was the most prevalent ACP-associated virus throughout nymph and adult ACP populations in Florida, detected in more than 60% of all samples tested, followed by DcPLV and DcFLV. HTS allowed us to identify a novel ACP-associated reo-like virus and a picorna-like virus. The putative reo-like virus, tentatively named Diaphorina citri cimodo-like virus, was later surveyed and detected back in seasonal adult and nymph ACP samples collected in Florida during this study. HTS generated data also revealed that the most abundant virus in Florida ACP populations was Citrus tristeza virus (CTV), which is not an ACP-associated virus, suggesting persistent presence of CTV infection in citrus throughout Florida groves. Collectively, information obtained from our study may be able to help guide the direction of biotechnological pest control efforts involving a number of viruses that were detected for the first time in Florida ACP populations, including two newly identified ACP-associated viruses.

Sensitivity and specificity of metatranscriptomics as an arbovirus surveillance tool

The ability to identify all the viruses within a sample makes metatranscriptomic sequencing an attractive tool to screen mosquitoes for arboviruses. Practical application of this technique, however, requires a clear understanding of its analytical sensitivity and specificity. To assess this, five dilutions (1:1, 1:20, 1:400, 1:8,000 and 1:160,000) of Ross River virus (RRV) and Umatilla virus (UMAV) isolates were spiked into subsamples of a pool of 100 Culex australicus mosquitoes. The 1:1 dilution represented the viral load of one RRV-infected mosquito in a pool of 100 mosquitoes. The subsamples underwent nucleic acid extraction, mosquito-specific ribosomal RNA depletion, and Illumina HiSeq sequencing. The viral load of the subsamples was also measured using reverse transcription droplet digital PCR (RT-ddPCR) and quantitative PCR (RT-qPCR). Metatranscriptomic sequencing detected both RRV and UMAV in the 1:1, 1:20 and 1:400 subsamples. A high specificity was achieved, with 100% of RRV and 99.6% of UMAV assembled contigs correctly identified. Metatranscriptomic sequencing was not as sensitive as RT-qPCR or RT-ddPCR; however, it recovered whole genome information and detected 19 other viruses, including four first detections for Australia. These findings will assist arbovirus surveillance programs in utilising metatranscriptomics in routine surveillance activities to enhance arbovirus detection.

An exploration of ambigrammatic sequences in narnaviruses

Narnaviruses have been described as positive-sense RNA viruses with a remarkably simple genome of ~3 kb, encoding only a highly conserved RNA-dependent RNA polymerase (RdRp). Many narnaviruses, however, are 'ambigrammatic' and harbour an additional uninterrupted open reading frame (ORF) covering almost the entire length of the reverse complement strand. No function has been described for this ORF, yet the absence of stops is conserved across diverse narnaviruses, and in every case the codons in the reverse ORF and the RdRp are aligned. The >3 kb ORF overlap on opposite strands, unprecedented among RNA viruses, motivates an exploration of the constraints imposed or alleviated by the codon alignment. Here, we show that only when the codon frames are aligned can all stop codons be eliminated from the reverse strand by synonymous single-nucleotide substitutions in the RdRp gene, suggesting a mechanism for de novo gene creation within a strongly conserved amino-acid sequence. It will be fascinating to explore what implications this coding strategy has for other aspects of narnavirus biology. Beyond narnaviruses, our rapidly expanding catalogue of viral diversity may yet reveal additional examples of this broadly-extensible principle for ambigrammatic-sequence development.

A molecular clone of Chronic Bee Paralysis Virus (CBPV) causes mortality in honey bee pupae (Apis mellifera)

Among the many diseases compromising the well-being of the honey bee (Apis mellifera) the chronic paralysis syndrome of adult honey bees is one of the best described. The causative agent, chronic bee paralysis virus (CBPV), is a positive sense, single-stranded RNA virus with a segmented genome. Segment 1 encodes three putative open reading frames (ORFs), including the RNA-dependent RNA polymerase and other non-structural protein coding regions. Segment 2 encodes four putative ORFs, which contain the genes of supposed structural proteins. In this study, we established a reverse genetic system for CBPV by molecular cloning of DNA copies of both genome segments. CBPV rescue was studied in imago and honey bee pupae infection models. Virus replication and progeny virus production was only initiated when capped RNAs of both genome segments were injected in honey bees. As injection of these clonal RNAs caused clinical symptoms similar to wild-type CBPV infection, we conclude that the novel molecular clone fulfilled Koch’s postulates. Our virus clone will enable in-depth analysis of CBPV pathogenesis and help to increase knowledge about this important honey bee disease.

A Novel Totivirus Naturally Occurring in Two Different Fungal Genera

Mycoviruses are widely distributed across different phyla of the fungal kingdom. Viruses that share significant sequence similarities have been reported in different fungi, suggesting descent from a common ancestor. In this study, two fungal genera isolated from the same sample, Trichoderma koningiopsis isolate Mg10 and Clonostachys rosea isolate Mg06, were reported to have identical double-stranded RNA (dsRNA) profiles that consist of two virus-like, dsRNA elements (dsRNA-L and dsRNA-S). The complete sequence and genome organization of dsRNA-L from isolate Mg10 was determined. It is 4712 nucleotides (nt) long and contains two non-overlapping open reading frames (ORFs) that code for proteins with similarities to totiviruses. Consequently the virus was given the proposed name Trichoderma koningiopsis totivirus 1 (TkTV1/Mg10). The TkTV1/Mg10 genome structure resembles that of yeast totiviruses in which the region preceding the stop codon of ORF1 contains the elements required for -1 ribosomal frameshifting which may induce the expression of an ORF1–ORF2 (CP-RdRp) fusion protein. Sequence analyses of viral dsRNA-L from C. rosea isolate Mg06 revealed that it is nearly identical with that of TkTV1/Mg10. This relatedness was confirmed by northern blot hybridization and indicates very recent natural horizontal transmission of this virus between unrelated fungi. TkTV1 purified isometric virions were ∼38–40 nm in diameter and were able to transfect T. koningiopsis and C. rosea protoplasts. This is another report of a mycovirus present naturally in two taxonomically distinct fungi.

Aedes spp. and Their Microbiota: A Review

Aedes spp. are a major public health concern due to their ability to be efficient vectors of dengue, Chikungunya, Zika, and other arboviruses. With limited vaccines available and no effective therapeutic treatments against arboviruses, the control of Aedes spp. populations is currently the only strategy to prevent disease transmission. Host-associated microbes (i.e., microbiota) recently emerged as a promising field to be explored for novel environmentally friendly vector control strategies. In particular, gut microbiota is revealing its impact on multiple aspects of Aedes spp. biology, including vector competence, thus being a promising target for manipulation. Here we describe the technological advances, which are currently expanding our understanding of microbiota composition, abundance, variability, and function in the two main arboviral vectors, the mosquitoes Aedes aegypti and Aedes albopictus. Aedes spp. microbiota is described in light of its tight connections with the environment, with which mosquitoes interact during their various developmental stages. Unraveling the dynamic interactions among the ecology of the habitat, the mosquito and the microbiota have the potential to uncover novel physiological interdependencies and provide a novel perspective for mosquito control.

De novo profiling of RNA viruses in Anopheles malaria vector mosquitoes from forest ecological zones in Senegal and Cambodia

BACKGROUND

Mosquitoes are colonized by a large but mostly uncharacterized natural virome of RNA viruses, and the composition and distribution of the natural RNA virome may influence the biology and immunity of Anopheles malaria vector populations.

RESULTS

Anopheles mosquitoes were sampled in malaria endemic forest village sites in Senegal and Cambodia, including Anopheles funestus, Anopheles gambiae group sp., and Anopheles coustani in Senegal, and Anopheles hyrcanus group sp., Anopheles maculatus group sp., and Anopheles dirus in Cambodia. The most frequent mosquito species sampled at both study sites are human malaria vectors. Small and long RNA sequences were depleted of mosquito host sequences, de novo assembled and clustered to yield non-redundant contigs longer than 500 nucleotides. Analysis of the assemblies by sequence similarity to known virus families yielded 115 novel virus sequences, and evidence supports a functional status for at least 86 of the novel viral contigs. Important monophyletic virus clades in the Bunyavirales and Mononegavirales orders were found in these Anopheles from Africa and Asia. The remaining non-host RNA assemblies that were unclassified by sequence similarity to known viruses were clustered by small RNA profiles, and 39 high-quality independent contigs strongly matched a pattern of classic RNAi processing of viral replication intermediates, suggesting they are entirely undescribed viruses. One thousand five hundred sixty-six additional high-quality unclassified contigs matched a pattern consistent with Piwi-interacting RNAs (piRNAs), suggesting that strand-biased piRNAs are generated from the natural virome in Anopheles. To functionally query piRNA effect, we analyzed piRNA expression in Anopheles coluzzii after infection with O'nyong nyong virus (family Togaviridae), and identified two piRNAs that appear to display specifically altered abundance upon arbovirus infection.

CONCLUSIONS

Anopheles vectors of human malaria in Africa and Asia are ubiquitously colonized by RNA viruses, some of which are monophyletic but clearly diverged from other arthropod viruses. The interplay between small RNA pathways, immunity, and the virome may represent part of the homeostatic mechanism maintaining virome members in a commensal or nonpathogenic state, and could potentially influence vector competence.

The discovery and global distribution of novel mosquito-associated viruses in the last decade (2007-2017)

In the last decade, virus hunting and discovery has gained pace. This achievement has been driven by three major factors: (a) advancements in sequencing technologies, (b) scaled-up routine arbovirus surveillance strategies, and (c) the "hunt" for emerging pathogens and novel viruses. Many novel viruses have been discovered from a myriad of hosts, vectors, and environmental samples. To help promote understanding of the global diversity and distribution of mosquito-associated viruses and facilitate future studies, we review mosquito-associated viruses discovered between years 2007 and 2017, across the world. In the analyzed period, novel mosquito-associated viruses belonging to 25 families and a general group of unclassified viruses were categorized. The top three discovered novel mosquito-associated viruses belonged to families Flaviviridae (n=32), Rhabdoviridae (n=16), and Peribunyaviridae (n=14). Also, 67 unclassified viruses were reported. Majority of these novel viruses were identified from Culex spp, Anopheles spp, Aedes spp, and Mansonia spp mosquitoes, respectively. Notably, the number of these discovered novels is not representative of intercontinental virus diversity but rather is influenced by the number of studies done in the study period. Some of these newly discovered mosquito-associated viruses have medical significance, either directly or indirectly. For instance, in the study period, 14 novel mosquito-borne viruses that infect mammalian cells in vitro were reported. These viruses pose a danger to the global health security on emerging viral diseases. On the other hand, some of the newly discovered insect specific viruses described herein have potential application as future biocontrol and vaccine agents against known pathogenic arboviruses. Overall, this review outlines the crucial role played by mosquitoes as viral vectors in the global virosphere.

Novel insights into endogenous RNA viral elements in Ixodes scapularis and other arbovirus vector genomes

Many emerging arboviruses are not transmitted by traditional mosquito vectors, but by lesser-studied arthropods such as ticks, midges, and sand flies. Small RNA (sRNA) silencing pathways are the main antiviral defence mechanism for arthropods, which lack adaptive immunity. Non-retroviral integrated RNA virus sequences (NIRVS) are one potential source of sRNAs which comprise these pathways. NIRVS are remnants of past germline RNA viral infections, where viral cDNA integrates into the host genome and is vertically transmitted. In Aedes mosquitoes, NIRVS are widespread and produce PIWI-interacting RNAs (piRNAs). These are hypothesised to target incoming viral transcripts to modulate viral titre, perhaps rendering the organism a more efficient arbovirus vector. To explore the NIRVS landscape in alternative arbovirus vectors, we validated the NIRVS landscape in Aedes spp. and then identified novel NIRVS in six medically relevant arthropods and also in Drosophila melanogaster. We identified novel NIRVS in Phlebotomus papatasi, Culicoides sonorensis, Rhipicephalus microplus, Anopheles gambiae, Culex quinquefasciatus, and Ixodes scapularis. Due to their unexpected abundance, we further characterised NIRVS in the blacklegged tick I. scapularis (n = 143). Interestingly, NIRVS are not enriched in R. microplus, another hard tick, suggesting this is an Ixodes-specific adaptation. I. scapularis NIRVS are enriched in bunya- and orthomyxo-like sequences, reflecting that ticks are a dominant host for these virus groups. Unlike in mosquitoes, I. scapularis NIRVS are more commonly derived from the non-structural region (replicase) of negative-sense viruses, as opposed to structural regions (e.g. glycoprotein). Like other arthropods, I. scapularis NIRVS preferentially integrate into genomic piRNA clusters, and serve as a template for primary piRNA production in the commonly used embryonic I. scapularis ISE6 cell line. Interestingly, we identified a two-fold enrichment of non-long terminal repeat (non-LTR) retrotransposons, in genomic proximity to NIRVS, contrasting with studeis in Ae. aegypti, where LTR retrotransposons are instead associated with NIRVS formation. We characterised NIRVS phylogeny and integration patterns in the important vector, I. scapularis, revealing they are distinct from those in Aedes spp. Future studies will explore the possible antiviral mechanism conferred by NIRVS to I. scapularis,which may help the transmission of pathogenic arboviruses. Finally, this study explored NIRVS as an untapped wealth of viral diversity in arthropods.

Mitovirus and Mitochondrial Coding Sequences from Basal Fungus Entomophthora muscae

Fungi constituting the Entomophthora muscae species complex (members of subphylum Entomophthoromycotina, phylum Zoopagamycota) commonly kill their insect hosts and manipulate host behaviors in the process. In this study, we made use of public transcriptome data to identify and characterize eight new species of mitoviruses associated with several different E. muscae isolates. Mitoviruses are simple RNA viruses that replicate in host mitochondria and are frequently found in more phylogenetically apical fungi (members of subphylum Glomeromyoctina, phylum Mucoromycota, phylum Basidiomycota and phylum Ascomycota) as well as in plants. E. muscae is the first fungus from phylum Zoopagomycota, and thereby the most phylogenetically basal fungus, found to harbor mitoviruses to date. Multiple UGA (Trp) codons are found not only in each of the new mitovirus sequences from E. muscae but also in mitochondrial core-gene coding sequences newly assembled from E. muscae transcriptome data, suggesting that UGA (Trp) is not a rarely used codon in the mitochondria of this fungus. The presence of mitoviruses in these basal fungi has possible implications for the evolution of these viruses.

The role of co-infection and swarm dynamics in arbovirus transmission

Arthropod-borne viruses (arboviruses) are transmitted by hematophagous insects, primarily mosquitoes. The geographic range and prevalence of mosquito-borne viruses and their vectors has dramatically increased over the last 50 years. As a result, the most medically important arboviurses now co-exist in many regions, resulting in an increased frequency of co-infections in hosts and vectors. In addition to concurrent infections with human pathogens, mosquito-only viruses and/or enzootic viruses not associated with human disease are ubiquitous in mosquito populations. Moreover, mosquito-borne viruses are largely RNA viruses that exist within individual hosts as a diverse and dynamic swarm of closely related genotypes. Interactions among co-infecting viruses and genotypes can have profound effects on virulence, fitness and evolution. Here, we review our understanding of how these complex interactions influence transmission of mosquito-borne viruses, focusing on the often-neglected virus interactions in the mosquito vector, and identify gaps in our knowledge that should guide future studies.