Identification and Characterization of Two Novel RNA Viruses from Anopheles gambiae Species Complex Mosquitoes

Analysis of somatic piRNAs in the malaria mosquito Anopheles coluzzii reveals atypical classes of genic small RNAs

Piwi-interacting small RNAs (piRNA) play a key role in controlling the activity of transposable elements (TEs) in the animal germline. In diverse arthropod species, including the pathogen vectors mosquitoes, the piRNA pathway is also active in nongonadal somatic tissues, where its targets and functions are less clear. Here, we studied the features of small RNA production in head and thorax tissues of an uninfected laboratory strain of Anopheles coluzzii focusing on the 24-32-nt-long RNAs. Small RNAs derived from repetitive elements constitute a minor fraction while most small RNAs process from long noncoding RNAs (lncRNAs) and protein-coding gene mRNAs. The majority of small RNAs derived from repetitive elements and lncRNAs exhibited typical piRNAs features. By contrast, majority of protein-coding gene-derived 24-32 nt small RNAs lack the hallmarks of piRNAs and have signatures of nontemplated 3' end tailing. Most of the atypical small RNAs exhibit female-biased expression and originate from mitochondrial and nuclear genes involved in energy metabolism. We also identified atypical genic small RNAs in Anopheles gambiae somatic tissues, which further validates the noncanonical mechanism of their production. We discuss a novel mechanism of small RNA production in mosquito somatic tissues and the possible functional significance of genic small RNAs.

Metagenome analysis of viruses associated with Anopheles mosquitoes from Ramu Upazila, Cox's Bazar District, Bangladesh

Bangladesh has a warm climate and landscapes favourable for the proliferation of mosquitoes. Mosquito-borne pathogens including malaria and arthropod-borne viruses (arboviruses) remain a serious threat to the public health requiring constant vector control and disease surveillance. From November 2018 to April 2019, Anopheles mosquitoes were collected in three unions in the Ramu Upazila (sub-district) of Cox's Bazar District, Bangladesh. The mosquito specimens were combined into pools based on date of collection, household ID, and sex. Metagenome next-generation sequencing was conducted to elucidate diversity of virus sequences in each pool. Homology-based taxonomic classification and phylogenetic analyses identified a broad diversity of putative viruses from 12 known families, with additional unclassified viruses also likely present. Analysis of male mosquitoes showed some of these viruses are likely capable of being vertically transmitted. Moreover, many of the assembled virus sequences share homology and phylogenetic affinity with segments in sequenced Anopheles genomes, and may represent endogenous viral elements derived from a past evolutionary relationship between these putative viruses and their mosquito hosts.

Viral metagenomics of hematophagous insects collected in the Carajas mining complex, Pará State, Brazil

Hematophagous insects are vectors of viruses that cause diseases in humans and animals worldwide. Mosquitoes (Culicidae), biting midges (Ceratopogonidae), and sandflies (Psychodidae) were collected in three municipalities (Marabá, Canaã dos Carajás, and Curionópolis) in the state of Pará, Brazil, in 2019. Morphological keys were used for the taxonomic identification of insect species. High-throughput sequencing and metagenomic analysis were employed to characterize the viromes of the hematophagous insects. We characterized the virome of 839 insects grouped into 14 pools. A total of 729 million paired reads were generated, with 12 million viral sequences (3 % of the reads). The families Reoviridae, Myoviridae, Retroviridae, and Poxviridae were found in all samples of this study. Phylogenies of RNA-dependent RNA polymerase (RdRp) from viruses of the families Chuviridae, Dicistroviridae, Flaviviridae, Iflaviridae, Mesoniviridae, Phenuiviridae, and Rhabdoviridae were performed. In this study, the first isolation of the Guaico Culex Virus (GCXV) in the northern region of Brazil was obtained from a pool of Culex (Melanoconion) spp. mosquitoes collected in Curionópolis. The data obtained in this study demonstrate that the Carajás region has an ecosystem rich in viruses. Additional studies are needed to understand the dynamics of viruses in vectors, vertebrates, and the human population in the region.

Translating mosquito viromes into vector management strategies

Mosquitoes are best known for transmitting human and animal viruses. However, they also harbour mosquito-specific viruses (MSVs) as part of their microbiota. These are a group of viruses whose diversity and prevalence overshadow their medically relevant counterparts. Although metagenomics sequencing has remarkably accelerated the discovery of these viruses, what we know about them is often limited to sequence information, leaving much of their fundamental biology to be explored. Understanding the biology and ecology of MSVs can enlighten our knowledge of virus-virus interactions and lead to new innovations in the management of mosquito-borne viral diseases. We retrace the history of their discovery and discuss research milestones that would line the path from mosquito virome knowledge to vector management strategies.

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.

Determinants of Chikungunya and O'nyong-Nyong Virus Specificity for Infection of Aedes and Anopheles Mosquito Vectors

Mosquito-borne diseases caused by viruses and parasites are responsible for more than 700 million infections each year. Anopheles and Aedes are the two major vectors for, respectively, malaria and arboviruses. Anopheles mosquitoes are the primary vector of just one known arbovirus, the alphavirus o'nyong-nyong virus (ONNV), which is closely related to the chikungunya virus (CHIKV), vectored by Aedes mosquitoes. However, Anopheles harbor a complex natural virome of RNA viruses, and a number of pathogenic arboviruses have been isolated from Anopheles mosquitoes in nature. CHIKV and ONNV are in the same antigenic group, the Semliki Forest virus complex, are difficult to distinguish via immunodiagnostic assay, and symptomatically cause essentially the same human disease. The major difference between the arboviruses appears to be their differential use of mosquito vectors. The mechanisms governing this vector specificity are poorly understood. Here, we summarize intrinsic and extrinsic factors that could be associated with vector specificity by these viruses. We highlight the complexity and multifactorial aspect of vectorial specificity of the two alphaviruses, and evaluate the level of risk of vector shift by ONNV or CHIKV.

Expansion of viral genomes with viral protein genome linked copies

Virus protein-linked genome (VPg) proteins are required for replication. VPgs are duplicated in a subset of RNA viruses however their roles are not fully understood and the extent of viral genomes containing VPg copies has not been investigated in detail. Here, we generated a novel bioinformatics approach to identify VPg sequences in viral genomes using hidden Markov models (HMM) based on alignments of dicistrovirus VPg sequences. From metagenomic datasets of dicistrovirus genomes, we identified 717 dicistrovirus genomes containing VPgs ranging from a single copy to 8 tandem copies. The VPgs are classified into nine distinct types based on their sequence and length. The VPg types but not VPg numbers per viral genome followed specific virus clades, thus suggesting VPgs co-evolved with viral genomes. We also identified VPg duplications in aquamavirus and mosavirus genomes. This study greatly expands the number of viral genomes that contain VPg copies and indicates that duplicated viral sequences are more widespread than anticipated.

Novel Dicistroviruses in an Unexpected Wide Range of Invertebrates

Dicistroviruses are members of a rapidly growing family of small RNA viruses. Related sequences have been discovered in many environmental samples, indicating that our knowledge about dicistrovirus diversity and host range is still limited. In this study, we performed a systematic search against the publicly available transcriptome database, and identified large numbers of dicistrovirus-like sequences in a wide variety of eukaryotic species. The origins of these sequences were 108 invertebrates (including 77 insect species belonging to 18 orders) and 11 plants, revealing new associations between dicistroviruses and hosts. Finally, 83 transcripts corresponding to nearly-complete viral genomes were retrieved from the RNA-seq data, of which most sequences showed limited similarity to known dicistroviruses and might present previously unreported virus species. Phylogenetic analysis suggested that horizontal virus transfer has occurred between diverse hosts and has important implications for dicistrovirus evolution. The results will provide new insight into the hidden diversity of the Dicistroviridae, and help us to better understand the viral evolution, host range and the possible way of transmission.

A library preparation optimized for metagenomics of RNA viruses

Our understanding of the viral communities associated to animals has not yet reached the level attained on the bacteriome. This situation is due to, among others, technical challenges in adapting metagenomics using high-throughput sequencing to the study of RNA viromes in animals. Although important developments have been achieved in most steps of viral metagenomics, there is yet a key step that has received little attention: the library preparation. This situation differs from bacteriome studies in which developments in library preparation have largely contributed to the democratisation of metagenomics. Here, we present a library preparation optimized for metagenomics of RNA viruses from insect vectors of viral diseases. The library design allows a simple PCR-based preparation, such as those routinely used in bacterial metabarcoding, that is adapted to shotgun sequencing as required in viral metagenomics. We first optimized our library preparation using mock viral communities and then validated a full metagenomic approach incorporating our preparation in two pilot studies with field-caught insect vectors; one including a comparison with a published metagenomic protocol. Our approach provided a fold increase in virus-like sequences compared to other studies, and nearly-full genomes from new virus species. Moreover, our results suggested conserved trends in virome composition within a population of a mosquito species. Finally, the sensitivity of our approach was compared to a commercial diagnostic PCR for the detection of an arbovirus in field-caught insect vectors. Our approach could facilitate studies on viral communities from animals and the democratization of metagenomics in community ecology of 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.

Next-Generation Sequencing on Insectivorous Bat Guano: An Accurate Tool to Identify Arthropod Viruses of Potential Agricultural Concern

Viruses belonging to the Dicistroviridae family have attracted a great deal of attention from scientists owing to their negative impact on agricultural economics, as well as their recent identification as potential aetiological agents of febrile illness in human patients. On the other hand, some Dicistroviruses are also studied for their potential biopesticide properties. To date, Dicistrovirus characterized in African mainland remain scarce. By using High-Throughput Sequencing technology on insectivorous bat faeces (Hipposideros Caffer) sampled in a cave used by humans to collect bat guano (bat manure) as fertilizer in Zimbabwe, we characterized the full-length sequences of three Dicistrovirus belonging to the Cripavirus and Aparavirus genus: Big Sioux River Virus-Like (BSRV-Like), Acute Bee Paralysis Virus (ABPV), and Aphid Lethal Paralysis Virus (ALPV). Phylogenetic analyses of ORF-1 and ORF-2 genes showed a complex evolutionary history between BSRV and close viruses, as well as for the Aparavirus genus. Herewith, we provide the first evidence of the presence of Dicistrovirus in Zimbabwe and highlight the need to further document the impact of such viruses on crops, as well as in beekeeping activities in Zimbabwe which represent a crucial source of income for Zimbabwean people.

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.

Small RNA responses of Culex mosquitoes and cell lines during acute and persistent virus infection

RNA interference is a crucial antiviral mechanism in arthropods, including in mosquito vectors of arthropod-borne viruses (arboviruses). Although the exogenous small interfering RNA (siRNA) pathway constitutes an efficient antiviral response in mosquitoes, virus-derived P-element induced wimpy testis (PIWI)-interacting RNAs (piRNAs) have been implicated in the response to alpha-, bunya- and flaviviruses in Aedes spp. mosquitoes. Culex mosquitoes transmit several medically important viruses including West Nile virus (WNV), but are considerably less well studied than Aedes mosquitoes and little is known about antiviral RNA interference in Culex mosquitoes. Therefore, we sequenced small RNA (sRNA) libraries from different Culex cell lines and tissues infected with WNV. The clear majority of virus-derived sRNA reads were 21 nt siRNAs in all cell lines and tissues tested, with no evidence for a role of WNV-derived piRNAs. Additionally, we aligned sRNA reads from Culex quinquefasciatus Hsu cells to the insect-specific rhabdovirus, Merida virus, which persistently replicates in these cells. We found that a significant proportion of the sRNA response to Merida virus consisted of piRNAs. Since viral DNA forms have been implicated in siRNA and piRNA responses of Aedes spp. mosquitoes, we also tested for viral DNA forms in WNV infected Culex cells. We detected viral DNA in Culex tarsalis cells infected with WNV and, to a lesser amount, WNV and Merida virus-derived DNA in Culex quinquefasciatus Hsu cells. In conclusion, Hsu cells generated Merida virus-derived piRNAs, but our data suggests that the major sRNA response of Culex cells and mosquitoes to WNV infection is the exogenous siRNA response. It is also evident that sRNA responses differ significantly between specific virus-mosquito combinations. Future work using additional Culex-borne viruses may further elucidate how virus-derived piRNAs are generated in Culex cells and what role they may play in controlling replication of different viruses.

Exploiting the Legacy of the Arbovirus Hunters

In recent years, it has become evident that a generational gap has developed in the community of arbovirus research. This apparent gap is due to the dis-investment of training for the next generation of arbovirologists, which threatens to derail the rich history of virus discovery, field epidemiology, and understanding of the richness of diversity that surrounds us. On the other hand, new technologies have resulted in an explosion of virus discovery that is constantly redefining the virosphere and the evolutionary relationships between viruses. This paradox presents new challenges that may have immediate and disastrous consequences for public health when yet to be discovered arboviruses emerge. In this review we endeavor to bridge this gap by providing a historical context for the work being conducted today and provide continuity between the generations. To this end, we will provide a narrative of the thrill of scientific discovery and excitement and the challenges lying ahead.

Interaction of RNA viruses of the natural virome with the African malaria vector, Anopheles coluzzii

Mosquitoes are colonized by a little-studied natural virome. Like the bacterial microbiome, the virome also probably influences the biology and immunity of mosquito vector populations, but tractable experimental models are lacking. We recently discovered two novel viruses in the virome of wild Anopheles and in colonies of the malaria vector Anopheles coluzzii: Anopheles C virus and Anopheles cypovirus. Here, we describe biological interactions between these two viruses and An. coluzzii mosquitoes. Viral abundance varies reproducibly during mosquito development. DNA forms of these viruses were not detected, and thus viral persistence is likely based on vertical transmission of RNA genomes. At least Anopheles C virus is vertically transmitted by an intraembryonic route. Relative abundance of the two viruses is inversely correlated in individual mosquitoes. One possible mechanism for this could be interactions with host immunity, and functional genomic analysis indicated differential influence of at least the Toll and JAK/STAT immune signaling pathways upon the viruses. The nonrandom distributions and interactions with host immunity suggest that these and other members of the natural virome may constitute a source of unrecognized heterogeneity in mosquito vector populations.

Mosquito Small RNA Responses to West Nile and Insect-Specific Virus Infections in Aedes and Culex Mosquito Cells

Small RNA mediated responses are essential for antiviral defence in mosquitoes, however, they appear to differ per virus-vector combination. To further investigate the diversity of small RNA responses against viruses in mosquitoes, we applied a small RNA deep sequencing approach on five mosquito cell lines: Culex tarsalis CT cells, Aedes albopictus U4.4 and C6/36 cells, Ae. aegypti Aag2 cells (cleared from cell fusing agent virus and Culex Y virus (CYV) by repetitive dsRNA transfections) and Ae. pseudoscutellaris AP-61 cells. De novo assembly of small RNAs revealed the presence of Phasi Charoen-like virus (PCLV), Calbertado virus, Flock House virus and a novel narnavirus in CT cells, CYV in U4.4 cells, and PCLV in Aag2 cells, whereas no insect-specific viruses (ISVs) were detected in C6/36 and AP-61 cells. Next, we investigated the small RNA responses to the identified ISVs and to acute infection with the arthropod-borne West Nile virus (WNV). We demonstrate that AP-61 and C6/36 cells do not produce siRNAs to WNV infection, suggesting that AP-61, like C6/36, are Dicer-2 deficient. CT cells produced a strong siRNA response to the persistent ISVs and acute WNV infection. Interestingly, CT cells also produced viral PIWI-interacting (pi)RNAs to PCLV, but not to WNV or any of the other ISVs. In contrast, in U4.4 and Aag2 cells, WNV siRNAs, and pi-like RNAs without typical ping-pong piRNA signature were observed, while this signature was present in PCLV piRNAs in Aag2 cells. Together, our results demonstrate that mosquito small RNA responses are strongly dependent on both the mosquito cell type and/or the mosquito species and family of the infecting virus.

Metagenomic Analysis of Flaviviridae in Mosquito Viromes Isolated From Yunnan Province in China Reveals Genes From Dengue and Zika Viruses

More than 6,000 mosquitoes of six species from six sites were collected and tested for their virome using metagenomics sequencing and bioinformatic analysis. The identified viral sequences belonged to more than 50 viral families. The results were verified by PCR of selected viruses in all mosquitoes, followed by phylogenetic analysis. In the present study, we identified the partial dengue virus (DENV), Zika virus (ZIKV), and Japanese encephalitis virus (JEV) sequences in mosquitoes. Metagenomic analysis and the PCR amplification revealed three DENV sequences, one of which encodes a partial envelope protein. Two ZIKV sequences both encoding partial nonstructural protein 3 and one JEV sequence encoding the complete envelope protein were identified. There was variability in the viral titers of the newly isolated virus JEV-China/YN2016-1 of different passage viruses. The newly identified Zika virus gene from ZIKV-China/YN2016-1 was an Asian genotype and shared the highest nucleotide sequence identity (97.1%) with a ZIKV sequence from Thailand isolated in 2004. Phylogenetic analysis of ZIKV-China/YN2016-1 and ZIKV-China/YN2016-2 with known Flavivirus genes indicated that ZIKV has propagated in Yunnan province, China.

Identification of a piscine reovirus-related pathogen in proliferative darkening syndrome (PDS) infected brown trout (Salmo trutta fario) using a next-generation technology detection pipeline

The proliferative darkening syndrome (PDS) is an annually recurring disease that causes species-specific die-off of brown trout (Salmo trutta fario) with a mortality rate of near 100% in pre-alpine rivers of central Europe. So far the etiology and causation of this disease is still unclear. The objective of this study was to identify the cause of PDS using a next-generation technology detection pipeline. Following the hypothesis that PDS is caused by an infectious agent, brown trout specimens were exposed to water from a heavily affected pre-alpine river with annual occurrence of the disease. Specimens were sampled over the entire time period from potential infection through death. Transcriptomic analysis (microarray) and RT-qPCR of brown trout liver tissue evidenced strong gene expression response of immune-associated genes. Messenger RNA of specimens with synchronous immune expression profiles were ultra-deep sequenced using next-generation sequencing technology (NGS). Bioinformatic processing of generated reads and gap-filling Sanger re-sequencing of the identified pathogen genome revealed strong evidence that a piscine-related reovirus is the causative organism of PDS. The identified pathogen is phylogenetically closely related to the family of piscine reoviruses (PRV) which are considered as the causation of different fish diseases in Atlantic and Pacific salmonid species such as Salmo salar and Onchorhynchus kisutch. This study also highlights that the approach of first screening immune responses along a timeline in order to identify synchronously affected stages in different specimens which subsequently were ultra-deep sequenced is an effective approach in pathogen detection. In particular, the identification of specimens with synchronous molecular immune response patterns combined with NGS sequencing and gap-filling re-sequencing resulted in the successful pathogen detection of PDS.

Metagenomic Sequencing From Mosquitoes in China Reveals a Variety of Insect and Human Viruses

We collected 8,700 mosquitoes in three sites in China, which belonged to seven species. Their viromes were tested using metagenomic sequencing and bioinformatic analysis. The abundant viral sequences were detected and annotated belonging to more than 50 viral taxonomic families. The results were verified by PCR, followed by phylogenetic analysis. In the present study, we identified partial viral genes of dengue virus (DENV), a novel circovirus (CCV), densovirus (DNV), Japanese encephalitis virus (JEV), and Wuhan mosquito virus (WMV) in mosquitoes. Metagenomic analysis and PCR amplification revealed three DENV sequences, which were as homologous to the NS3 gene of DENV from Singapore isolated in 2005, with at least 91% nucleotide (nt) identity. Seven fragments of JEV encoding structural proteins were identified belonging to genotype I. They all shared high homology with structural protein genes of JEV isolated from Laos in 2009. The production of infectious virus particles of the newly isolated virus YunnanJEV2017-4 increased after passage from the BHK-21 cell line to the Vero cell line. Novel circovirus-related genes were identified and as being related to an unnamed gene of a mosquito circovirus (MCCV) sequence from the USA isolated in 2011, with at least 41% nt identity: this distant relationship suggests that the parent virus might belong to a novel circovirus genus. Additionally, numerous known viruses and some unknown viruses were also detected in mosquitoes from Yunnan province, China, which will be tested for propagation.

Highly focused transcriptional response of Anopheles coluzzii to O'nyong nyong arbovirus during the primary midgut infection

Background

Anopheles mosquitoes are efficient vectors of human malaria, but it is unknown why they do not transmit viruses as well as Aedes and Culex mosquitoes. The only arbovirus known to be consistently transmitted by Anopheles mosquitoes is O’nyong nyong virus (ONNV, genus Alphavirus, family Togaviridae). The interaction of Anopheles mosquitoes with RNA viruses has been relatively unexamined.

Results

We transcriptionally profiled the African malaria vector, Anopheles coluzzii, infected with ONNV. Mosquitoes were fed on an infectious bloodmeal and were analyzed by Illumina RNAseq at 3 days post-bloodmeal during the primary virus infection of the midgut epithelium, before systemic dissemination. Virus infection triggers transcriptional regulation of just 30 host candidate genes. Most of the regulated candidate genes are novel, without known function. Of the known genes, a significant cluster includes candidates with predicted involvement in carbohydrate metabolism. Two candidate genes encoding leucine-rich repeat immune (LRIM) factors point to possible involvement of immune protein complexes in the mosquito antiviral response. The primary ONNV infection by bloodmeal shares little transcriptional response in common with ONNV infection by intrathoracic injection, nor with midgut infection by the malaria parasites, Plasmodium falciparum or P. berghei. Profiling of A. coluzzii microRNA (miRNA) identified 118 known miRNAs and 182 potential novel miRNA candidates, with just one miRNA regulated by ONNV infection. This miRNA was not regulated by other previously reported treatments, and may be virus specific. Coexpression analysis of miRNA abundance and messenger RNA expression revealed discrete clusters of genes regulated by Imd and JAK/STAT, immune signaling pathways that are protective against ONNV in the primary infection.

Conclusions

ONNV infection of the A. coluzzii midgut triggers a remarkably limited gene regulation program of mostly novel candidate genes, which likely includes host genes deployed for antiviral defense, as well as genes manipulated by the virus to facilitate infection. Functional dissection of the ONNV-response candidate genes is expected to generate novel insight into the mechanisms of virus-vector interaction.

Electronic supplementary material

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

A Systematic Review of the Natural Virome of Anopheles Mosquitoes

Anopheles mosquitoes are vectors of human malaria, but they also harbor viruses, collectively termed the virome. The Anopheles virome is relatively poorly studied, and the number and function of viruses are unknown. Only the o’nyong-nyong arbovirus (ONNV) is known to be consistently transmitted to vertebrates by Anopheles mosquitoes. A systematic literature review searched four databases: PubMed, Web of Science, Scopus, and Lissa. In addition, online and print resources were searched manually. The searches yielded 259 records. After screening for eligibility criteria, we found at least 51 viruses reported in Anopheles, including viruses with potential to cause febrile disease if transmitted to humans or other vertebrates. Studies to date have not provided evidence that Anopheles consistently transmit and maintain arboviruses other than ONNV. However, anthropophilic Anopheles vectors of malaria are constantly exposed to arboviruses in human bloodmeals. It is possible that in malaria-endemic zones, febrile symptoms may be commonly misdiagnosed. It is also possible that anophelines may be inherently less competent arbovirus vectors than culicines, but if true, the biological basis would warrant further study. This systematic review contributes a context to characterize the biology, knowledge gaps, and potential public health risk of Anopheles viruses.

Genetic characterization of a novel picorna-like virus in Culex spp. mosquitoes from Mozambique

BACKGROUND

Mosquitoes are the potential vectors for a variety of viruses that can cause diseases in the human and animal populations. Viruses in the order Picornavirales infect a broad range of hosts, including mosquitoes. In this study, we aimed to characterize a novel picorna-like virus from the Culex spp. of mosquitoes from the Zambezi Valley of Mozambique.

METHODS

The extracted RNA from mosquito pools was pre-amplified with the sequence independent single primer amplification (SISPA) method and subjected to high-throughput sequencing using the Ion Torrent platform. Reads that are classified as Iflaviridae, Picornaviridae and Dicistroviridae were assembled by CodonCode Aligner and SPAdes. Gaps between the viral contigs were sequenced by PCR. The genomic ends were analyzed by 5' and 3' RACE PCRs. The ORF was predicted with the NCBI ORF finder. The conserved domains were identified with ClustalW multiple sequence alignment, and a phylogenetic tree was built with MEGA. The presence of the virus in individual mosquito pools was detected by RT-PCR assay.

RESULTS

A near full-length viral genome (9740 nt) was obtained in Culex mosquitoes that encoded a complete ORF (3112 aa), named Culex picorna-like virus (CuPV-1). The predicted ORF had 38% similarity to the Hubei picorna-like virus 35. The sequence of the conserved domains, Helicase-Protease-RNA-dependent RNA polymerase, were identified by multiple sequence alignment and found to be at the 3' end, similar to iflaviruses. Phylogenetic analysis of the putative RdRP amino acid sequences indicated that the virus clustered with members of the Iflaviridae family. CuPV-1 was detected in both Culex and Mansonia individual pools with low infection rates.

CONCLUSIONS

The study reported a highly divergent, near full-length picorna-like virus genome from Culex spp. mosquitoes from Mozambique. The discovery and characterization of novel viruses in mosquitoes is an initial step, which will provide insights into mosquito-virus interaction mechanisms, genetic diversity and evolution.

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.

Genomes of viral isolates derived from different mosquitos species

Eleven viral isolates derived mostly in albopictus C6/36 cells from mosquito pools collected in Southeast Asia and the Americas between 1966 and 2014 contained particles with electron microscopy morphology typical of reoviruses. Metagenomics analysis yielded the near complete genomes of three novel reoviruses, Big Cypress orbivirus, Ninarumi virus, and High Island virus and a new tetravirus, Sarawak virus. Strains of previously characterized Sathuvarachi, Yunnan, Banna and Parry's Lagoon viruses (Reoviridae), Bontang virus (Mesoniviridae), and Culex theileri flavivirus (Flaviviridae) were also characterized. The availability of these mosquito virus genomes will facilitate their detection by metagenomics or PCR to better determine their geographic range, extent of host tropism, and possible association with arthropod or vertebrate disease.

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.

Advances in the application of high-throughput sequencing in invertebrate virology

Over the last decade, advances in high-throughput sequencing technologies have revolutionised biological research, making it possible for DNA/RNA sequencing of any organism of interest to be undertaken. Sequencing approaches are now routinely used in the detection and characterisation of (novel) viruses, investigation of host-pathogen interactions, and effective development of disease treatment strategies. For the sequencing and identification of viruses of interest, metagenomics approaches using infected host tissue are frequently used, as it is not always possible to culture and isolate these pathogens. High-throughput sequencing can also be used to investigate host-pathogen interactions by investigating (temporal) transcriptomic responses of both the host and virus, potentially leading to the discovery of novel opportunities for treatment and drug targets. In addition, viruses in environmental samples (e.g. water or soil samples) can be identified using eDNA/metagenomics approaches. The promise that recent developments in sequencing brings to the field of invertebrate virology are not devoid of technical challenges, including the need for better laboratory and bioinformatics strategies to sequence and assemble virus genomes within complex tissue or environmental samples, and the difficulties associated with the annotation of the large number of novel viruses being discovered.

Metavisitor, a Suite of Galaxy Tools for Simple and Rapid Detection and Discovery of Viruses in Deep Sequence Data

Metavisitor is a software package that allows biologists and clinicians without specialized bioinformatics expertise to detect and assemble viral genomes from deep sequence datasets. The package is composed of a set of modular bioinformatic tools and workflows that are implemented in the Galaxy framework. Using the graphical Galaxy workflow editor, users with minimal computational skills can use existing Metavisitor workflows or adapt them to suit specific needs by adding or modifying analysis modules. Metavisitor works with DNA, RNA or small RNA sequencing data over a range of read lengths and can use a combination of de novo and guided approaches to assemble genomes from sequencing reads. We show that the software has the potential for quick diagnosis as well as discovery of viruses from a vast array of organisms. Importantly, we provide here executable Metavisitor use cases, which increase the accessibility and transparency of the software, ultimately enabling biologists or clinicians to focus on biological or medical questions.

Mayaro virus: a forest virus primed for a trip to the city?

Mayaro virus (MAYV) is an emerging arthropod-borne virus (arbovirus). Infection by MAYV can produce Mayaro virus disease (MAYVD) which is usually a clinically diagnosed, acute, febrile illness associated with prolonged and painful joint inflammation and swelling. MAYVD may be clinically indistinguishable from dengue, chikungunya fever, malaria, rabies, measles or other arboviral diseases. The full spectrum of disease, sequelae, routes of infection, virus shedding and any rarer means of transmission remain undefined. MAYVD cases in humans have so far been localised to Central and South America, particularly regions in and around the Amazon basin. MAYV usually circulates in a sylvan cycle of forest mosquitoes and vertebrates, however it has also been found in more urban locations alongside anthropophilic (preferring humans) insect vectors. If transmission via anthropophilic mosquitoes becomes more efficient following viral change, or existing vectors change their habitat and biting habits, the risk of urban establishment and further spread into non-forested areas will grow. Surveillance, testing and vector control remain key to monitoring and preventing global spread and establishment. The possibility of MAYV becoming further urbanized is worthy of note, consideration and action to ensure MAYV does not spread beyond the forests and establish in the world's cities.

West African Anopheles gambiae mosquitoes harbor a taxonomically diverse virome including new insect-specific flaviviruses, mononegaviruses, and totiviruses

Anopheles gambiae are a major vector of malaria in sub-Saharan Africa. Viruses that naturally infect these mosquitoes may impact their physiology and ability to transmit pathogens. We therefore used metagenomics sequencing to search for viruses in adult Anopheles mosquitoes collected from Liberia, Senegal, and Burkina Faso. We identified a number of virus and virus-like sequences from mosquito midgut contents, including 14 coding-complete genome segments and 26 partial sequences. The coding-complete sequences define new viruses in the order Mononegavirales, and the families Flaviviridae, and Totiviridae. The identification of a flavivirus infecting Anopheles mosquitoes broadens our understanding of the evolution and host range of this virus family. This study increases our understanding of virus diversity in general, begins to define the virome of a medically important vector in its natural setting, and lays groundwork for future studies examining the potential impact of these viruses on anopheles biology and disease transmission.