A novel rhabdovirus, related to Merida virus, in field-collected mosquitoes from Anatolia and Thrace

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.

Metagenomics analysis reveals presence of the Merida-like virus in Georgia

Arbovirus surveillance is fundamental for the discovery of novel viruses and prevention of febrile vector-borne illnesses. Vector-borne pathogens can rapidly expand and adapt in new geographic and environmental conditions. In this study, metagenomic surveillance was conducted to identify novel viruses in the Country of Georgia. A total of 521 mosquitoes were captured near a military training facility and pooled from species Culex pipiens (Linnaeus) (87%) and Aedes albopictus (Skuse) (13%). We decided to further analyze the Culex pipiens mosquitoes, due to the more extensive number of samples collected. Our approach was to utilize an unbiased total RNA-seq for pathogen discovery in order to explore the mosquito virome. The viral reads from this analysis were mostly aligned to Insect-specific viruses from two main families, the Iflaviridae; a positive-stranded RNA virus and the Rhabdoviridae; a negative- and single-stranded RNA virus. Our pathogen discovery analysis revealed viral reads aligning to the Merida-like virus Turkey (MERDLVT) strain among the Rhabdoviridae. To further validate this result, we conducted a BLAST sequence comparison analysis of our samples with the MERDLVT strain. Our positive samples aligned to the MERDLVT strain with 96-100% sequence identity and 99.7-100% sequence coverage. A bootstrapped maximum-likelihood phylogenetic tree was used to evaluate the evolutionary relationships among these positive pooled specimens with the (MERDLVT) strain. The Georgia samples clustered most closely with two strains from Turkey, the Merida-like virus KE-2017a isolate 139-1-21 and the Merida-like virus Turkey isolate P431. Collectively, these results show the presence of the MERDLVT strain in Georgia.

Metavirome Analysis and Identification of Midge-Borne Viruses from Yunnan Province, China, in 2021

Midges are widely distributed globally and can transmit various human and animal diseases through blood-sucking. As part of this study, 259,300 midges were collected from four districts in Yunnan province, China, to detect the viral richness and diversity using metavirome analysis techniques. As many as 26 virus families were detected, and the partial sequences of bluetongue virus (BTV), dengue virus (DENV), and Getah virus (GETV) were identified by phylogenetic analysis and PCR amplification. Two BTV gene fragments, 866 bps for the VP2 gene of BTV type 16 and 655 bps for the VP5 gene of BTV type 21, were amplified. The nucleotide sequence identities of the two amplified BTV fragments were 94.46% and 98.81%, respectively, with two classical BTV-16 (GenBank: JN671907) and BTV-21 strains (GenBank: MK250961) isolated in Yunnan province. Furthermore, the BTV-16 DH2021 strain was successfully isolated in C6/36 cells, and the peak value of the copy number reached 3.13 × 107 copies/μL after five consecutive BHK-21 cell passages. Moreover, two 2054 bps fragments including the E gene of DENV genotype Asia II were amplified and shared the highest identity with the DENV strain isolated in New Guinea in 1944. A length of 656 bps GETV gene sequence encoded the partial capsid protein, and it shared the highest identity of 99.68% with the GETV isolated from Shandong province, China, in 2017. Overall, this study emphasizes the importance of implementing prevention and control strategies for viral diseases transmitted by midges in China.

Viruses of Freshwater Mussels during Mass Mortality Events in Oregon and Washington, USA

Freshwater mussels (Unionida) are globally imperiled, in part due to largely unexplained mass mortality events (MMEs). While recent studies have begun to investigate the possibility that mussel MMEs in the Eastern USA may be caused by infectious diseases, mussels in the Western USA have received relatively little attention in this regard. We conducted a two-year epidemiologic investigation of the role of viruses in ongoing MMEs of the Western pearlshell (Margaritifera falcata) and the Western ridged mussel (Gonidea angulata) in the Chehalis River and Columbia River watersheds in the Western USA. We characterized viromes of mussel hemolymph from 5 locations in 2018 and 2020 using metagenomic methods and identified 557 viruses based on assembled contiguous sequences, most of which are novel. We also characterized the distribution and diversity of a previously identified mussel Gammarhabdovirus related to pathogenic finfish viruses. Overall, we found few consistent associations between viruses and mussel health status. Variation in mussel viromes was most strongly driven by location, with little influence from date, species, or health status, though these variables together only explained ~1/3 of variation in virome composition. Our results demonstrate that Western freshwater mussels host remarkably diverse viromes, but no single virus or combination of viruses appears to be associated with morbidity or mortality during MMEs. Our findings have implications for the conservation of imperiled freshwater mussels, including efforts to enhance natural populations through captive propagation.

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.

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.

2021 Taxonomic update of phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales

In March 2021, following the annual International Committee on Taxonomy of Viruses (ICTV) ratification vote on newly proposed taxa, the phylum Negarnaviricota was amended and emended. The phylum was expanded by four families (Aliusviridae, Crepuscuviridae, Myriaviridae, and Natareviridae), three subfamilies (Alpharhabdovirinae, Betarhabdovirinae, and Gammarhabdovirinae), 42 genera, and 200 species. Thirty-nine species were renamed and/or moved and seven species were abolished. This article presents the updated taxonomy of Negarnaviricota as now accepted by the ICTV.

Metatranscriptomic analysis identifies different viral-like sequences in two neotropical Mansoniini mosquito species

Mosquitoes interact with a wide range of viruses including both arboviruses and insect-specific viruses. This study aimed to characterize the RNA viruses that are interacting with Mansonia wilsoni and Coquillettidia hermanoi mosquito species. The total RNA extracted from mosquito pools were sequenced on a Ion torrent platform. Viral contigs were identified against viral databases and their evolutionary relationship were reconstructed. We identified a total of 107 viral sequences, 11 of which were assigned as endogenous viral elements, and at least six known viral families were identified. Phylogenetic reconstructions were performed for 4 viral families. All Mansoniini viruses investigated through phylogenetic analysis are closely related to insect-specific viruses found in other mosquito species although with considerable divergence at the amino acid level, suggesting that we have detected new viral lineages. This study enhanced our understanding about the virome of two sylvatic Mansoniini mosquitoes.

Virome Diversity among Mosquito Populations in a Sub-Urban Region of Marseille, France

Some mosquito species have significant public health importance given their ability to transmit major diseases to humans and animals, making them the deadliest animals in the world. Among these, the Aedes (Ae.) genus is a vector of several viruses such as Dengue, Chikungunya, and Zika viruses that can cause serious pathologies in humans. Since 2004, Ae. albopictus has been encountered in the South of France, and autochthonous cases of Dengue, Chikungunya, and Zika diseases have recently been reported, further highlighting the need for a comprehensive survey of the mosquitoes and their associated viruses in this area. Using high throughput sequencing (HTS) techniques, we report an analysis of the DNA and RNA viral communities of three mosquito species Ae. albopictus, Culex (Cx.) pipiens, and Culiseta (Cs.) longiareolata vectors of human infectious diseases in a small sub-urban city in the South of France. Results revealed the presence of a significant diversity of viruses known to infect bacteria, plants, insects, and mammals. Several novel viruses were detected, including novel members of the Rhabdoviridae, Totiviridae, Iflaviviridae, Circoviridae, and Sobemoviridae families. No sequence related to major zoonotic viruses transmitted by mosquitoes was detected. The use of HTS on arthropod vector populations is a promising strategy for monitoring the emergence and circulation of zoonoses and epizooties. This study is a contribution to the knowledge of the mosquito microbiome.

Vector-borne viruses in Turkey: A systematic review and bibliography

Turkey serves as a natural hub for the dissemination of vector-borne viruses and provides many suitable habitats with diverse ecologies for introduction and establishment of new pathogens. This manuscript provides an updated systematic review and meta-analysis of the vector-borne viruses documented in Turkey. Following web-based identification, screening and eligibility evaluation, 291 published reports were reviewed. The publications were categorized and listed as a supplementary bibliography accompanying the manuscript. In brief, Crimean-Congo hemorrhagic fever virus (CCHFV) and West Nile virus (WNV) are currently documented as prominent tick and mosquito-borne viral pathogens in Turkey. CCHFV produces a significant number of infections annually, with severe outcome or death in a portion of cases. WNV gained attention following the clustering of cases in 2010. Exposure and infections with sandfly-borne phleboviruses, such as Toscana virus, are indigenous and widespread. Epidemiology, risk factors, symptomatic infections in susceptible hosts, vectors and reservoirs for these pathogens have been explored in detail. Detection of novel viruses in mosquitoes, sandflies and ticks from several regions is of particular interest, despite scarce information on their epidemiology and pathogenicity in vertebrates. Introduction and emergence of viruses transmitted by invasive Aedes mosquitoes constitute a threat, albeit only imported infections have so far been documented. Detection of Rift valley fever virus exposure is also of concern, due to its detrimental effects on livestock and spillover infections in humans. Vigilance to identify and diagnose probable cases as well as vector surveillance for established and potential pathogens is therefore, imperative.

Depicting the RNA Virome of Hematophagous Arthropods from Belgrade, Serbia

Hematophagous arthropods are important vectors for zoonotic pathogens. To date, a huge number of viruses have been identified in these arthropods, with a considerable proportion of them being human pathogens. However, the viromes of hematophagous arthropods are still largely unresearched. In this study, a number of arthropods were collected from Belgrade, Serbia including mosquitoes, ticks and bedbugs. The viromes of these arthropods were identified and characterized using Illumina MiSeq sequencing. In total, 21 viruses belonging to 11 families were characterized, with 11 of them representing novel species. These results may contribute to our knowledge of RNA viruses in arthropods and the discovery of novel human pathogens.

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.

Arboviral screening of invasive Aedes species in northeastern Turkey: West Nile virus circulation and detection of insect-only viruses

BACKGROUND

The recent reports of Aedes aegypti and Ae. albopictus populations in Turkey, in parallel with the territorial expansion identified in several surrounding countries, have raised concerns about the establishment and re-establishment of these invasive Aedes mosquitoes in Turkey. This cross-sectional study was performed to detect Aedes aegypti and Ae. albopictus in regions of recent incursions, and screen for viral pathogens known to be transmitted elsewhere by these species.

METHODOLOGY

Mosquitoes were collected at several locations in Artvin, Rize and Trabzon provinces of the Black Sea region during 2016-2017, identified morphologically, pooled and analyzed via generic or specific nucleic acid amplification assays. Viruses in positive pools were identified by product sequencing, cell culture inoculation and next generation sequencing (NGS) in selected specimens.

PRINCIPAL FINDINGS

The study group comprised 791 specimens. Aedes albopictus was the most abundant species in all locations (89.6%), followed by Ae. aegypti (7.8%) and Culex pipiens (2.5%). Mosquitoes were screened for viruses in 65 pools where fifteen (23.1%) were reactive. The infecting strains was identified as West Nile virus (WNV) in 5 pools (7.7%) with Ae. albopictus or Cx. pipiens mosquitoes. The obtained WNV sequences phylogenetically grouped with local and global lineage 1 clade 1a viruses. In 4 (6.2%) and 6 (9.2%) pools, respectively, cell fusing agent virus (CFAV) and Aedes flavivirus (AEFV) sequences were characterized. NGS provided a near-complete AEFV genome in a pool of Ae. albopictus. The strain is provisionally called "AEFV-Turkey", and functional analysis of the genome revealed several conserved motifs and regions associated with virus replication. Merida-like virus Turkey (MERDLVT), a recently-described novel rhabdovirus, was also co-detected in a Cx. pipiens pool also positive for WNV.

CONCLUSIONS/SIGNIFICANCE

Invasive Aedes mosquitoes are established in certain locations of northeastern Turkey. Herein we conclusively show the role of these species in WNV circulation in the region. Biosurveillance is imperative to monitor the spread of these species further into Asia Minor and to detect possible introduction of pathogens.

Genomic characterization, phylogenetic position and in situ localization of a novel putative mononegavirus in Lepeophtheirus salmonis

The complete genome sequence of a novel mononegavirus, Lepeophtheirus salmonis negative-stranded RNA virus 1 (LsNSRV-1), obtained from a salmonid ectoparasite, Lepeophtheirus salmonis was determined. The viral genome contains five open reading frames encoding three unknown proteins (ORF I, II and III), a putative glycoprotein (G), and a large (L) protein. Phylogenetic analysis placed LsNSRV-1 in the recently established mononegaviral family Artoviridae. LsNSRV-1 showed a prevalence of around 97% and was detected in all L. salmonis developmental stages. Viral genomic and antigenomic RNA was localized to nerve tissue, connective tissue, epithelial cells of the gut, subepidermal tissue, exocrine and cement glands, as well as the testis, vas deferens and spermatophore sac of male L. salmonis and the ovaries and oocytes of females. Viral RNA was detected in both the cytoplasm and the nucleoli of infected cells, and putative nuclear export and localization signals were found within the ORF I, III and L proteins, suggesting nuclear replication of LsNSRV-1. RNA interference (RNAi) was induced twice during development by the introduction of a double-stranded RNA fragment of ORF I, resulting in a transient knockdown of viral RNA. A large variation in the knockdown level was seen in adult males and off springs of knockdown animals, whereas the RNA level was more stable in adult females. Together with the localization of viral RNA within the male spermatophore and female oocytes and the amplification of viral RNA in developing embryos, this suggests that LsNSRV-1 is transmitted both maternally and paternally. Small amounts of viral RNA were detected at the site where chalimi were attached to the skin of Atlantic salmon (Salmo salar). However, as the RNAi-mediated treatment did not result in LsNSRV-1-negative offspring and the virus failed to replicate in the tested fish cell cultures, it is difficult to investigate the influence of secreted LsNSRV-1 on the salmon immune response.

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.

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.

West Nile virus, Anopheles flavivirus, a novel flavivirus as well as Merida-like rhabdovirus Turkey in field-collected mosquitoes from Thrace and Anatolia

Mosquitoes are involved in the transmission and maintenance of several viral diseases with significant health impact. Biosurveillance efforts have also revealed insect-specific viruses, observed to cocirculate with pathogenic strains. This report describes the findings of flavivirus and rhabdovirus screening, performed in eastern Thrace and Aegean region of Anatolia during 2016, including and expanding on locations with previously-documented virus activity. A mosquito cohort of 1545 individuals comprising 14 species were collected and screened in 108 pools via generic and specific amplification and direct metagenomics by next generation sequencing. Seven mosquito pools (6.4%) were positive in the flavivirus screening. West Nile virus lineage 1 clade 1a sequences were characterized in a pool Culex pipiens sensu lato specimens, providing the initial virus detection in Aegean region following 2010 outbreak. In an Anopheles maculipennis sensu lato pool, sequences closely-related to Anopheles flaviviruses were obtained, with similarities to several African and Australian strains of this new insect-specific flavivirus clade. In pools comprising Uranotaenia unguiculata (n=3), Cx. pipiens s.l. (n=1) and Aedes caspius (n=1) mosquitoes, sequences of a novel flavivirus, distantly-related to Flavivirus AV2011, identified previously in Spain and Turkey, were characterized. Moreover, DNA forms of the novel flavivirus were detected in two Ur. unguiculata pools. These sequences were highly-similar to the sequences amplified from viral RNA, with undisrupted reading frames, suggest the occurrence of viral DNA forms in natural conditions within mosquito hosts. Rhabdovirus screening revealed sequences of a recently-described novel virus, named the Merida-like virus Turkey (MERDLVT) in 5 Cx. pipiens s.l. pools (4.6%). Partial L and N gene sequences of MERDLVT were well-conserved among strains, with evidence for geographical clustering in phylogenetic analyses. Metagenomics provided the near-full genomic sequence in a specimen, revealing an identical genome organization and limited divergence from the prototype MERDLVT isolate.