Next Article in Journal
Genotyping Porcine Circovirus 3 (PCV-3) Nowadays: Does It Make Sense?
Previous Article in Journal
Molecular Insights into Host and Vector Manipulation by Plant Viruses
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 3
10.3390/v12030264
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deciphering the Virome of Culex vishnui Subgroup Mosquitoes, the Major Vectors of Japanese Encephalitis, in Japan

by
Astri Nur Faizah
1,2,
Daisuke Kobayashi
2,3,
Haruhiko Isawa
2,*,
Michael Amoa-Bosompem
2,4,
Katsunori Murota
2,5,
Yukiko Higa
2,
Kyoko Futami
6,
Satoshi Shimada
7,
Kyeong Soon Kim
8,
Kentaro Itokawa
9,
Mamoru Watanabe
2,
Yoshio Tsuda
2,
Noboru Minakawa
6,
Kozue Miura
1,
Kazuhiro Hirayama
1,* and
Kyoko Sawabe
2
1
Laboratory of Veterinary Public Health, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 1-1-1 Yayoi, Bunkyo-ku, Tokyo 113-8657, Japan
2
Department of Medical Entomology, National Institute of Infectious Diseases, 1-23-1 Toyama, Shinjuku-ku, Tokyo 162-8640, Japan
3
Department of Research Promotion, Japan Agency for Medical Research and Development, 20F Yomiuri Shimbun Bldg. 1-7-1 Otemachi, Chiyoda-ku, Tokyo 100-0004, Japan
4
Department of Environmental Parasitology, Tokyo Medical and Dental University, 1-5-45 Yushima, Bunkyo-ku, Tokyo 113-8510, Japan
5
Kyushu Research Station, National Institute of Animal Health, NARO, 2702 Chuzan, Kagoshima 891-0105, Japan
6
Department of Vector Ecology and Environment, Institute of Tropical Medicine, Nagasaki University, 1-12-4 Sakamoto, Nagasaki 852-8523, Japan
7
Department of Virology, Institute of Tropical Medicine, Nagasaki University, 1-12-4 Sakamoto, Nagasaki 852-8523, Japan
8
Joint Department of Veterinary Medicine, Faculty of Agriculture, Tottori University, 4-101 Koyama-cho Minami, Tottori 680-8550, Japan
9
Pathogen Genomics Center, National Institute of Infectious Diseases, 1-23-1 Toyama, Shinjuku-ku, Tokyo 162-8640, Japan
*
Authors to whom correspondence should be addressed.
Viruses 2020, 12(3), 264; https://0-doi-org.brum.beds.ac.uk/10.3390/v12030264
Submission received: 17 January 2020 / Revised: 18 February 2020 / Accepted: 26 February 2020 / Published: 28 February 2020
(This article belongs to the Section Insect Viruses)
Browse Figures
Versions Notes

Abstract

:
Japanese encephalitis (JE) remains a public health concern in several countries, and the Culex mosquito plays a central role in its transmission cycle. Culex mosquitoes harbor a wide range of viruses, including insect-specific viruses (ISVs), and can transmit a variety of arthropod-borne viruses (arboviruses) that cause human and animal diseases. The current trend of studies displays enhanced efforts to characterize the mosquito virome through bulk RNA sequencing due to possible arbovirus–ISV interactions; however, the extent of viral diversity in the mosquito taxon is still poorly understood, particularly in some disease vectors. In this study, arboviral screening and RNA virome analysis of Culex tritaeniorhynchus and C. pseudovishnui, which are part of the Culex vishnui subgroup mosquitoes, were performed. Results from these two mosquito species, known as the major vectors of JE virus (JEV) in Asia, collected in three prefectures in Japan were also compared with the sympatric species C. inatomii. A total of 27 viruses, including JEV, were detected from these Culex mosquitoes. Molecular and phylogenetic analyses of the detected viruses classified 15 of the 27 viruses as novel species, notably belonging to the Flaviviridae, Rhabdoviridae, Totiviridae, and Iflaviridae families. The successful isolation of JEV genotype I confirmed its continuous presence in Japan, suggesting the need for periodic surveillance. Aside from JEV, this study has also reported the diversity of the RNA virome of disease vectors and broadened the knowledge on mosquito virome profiles containing both arbovirus and ISV. Mosquito taxon seemed to contribute largely to the virome structure (e.g., virome composition, diversity, and abundance) as opposed to the geographical location of the mosquito species. This study therefore offers notable insights into the ecology and evolution of each identified virus and viral family. To the authors’ knowledge, this is the first study to characterize the viromes of the major JE vectors in Japan.

1. Introduction

Approximately 67,900 Japanese encephalitis (JE) cases typically occur annually in 24 JE-endemic countries, with an overall incidence of 1.8 per 100,000 [1]. As the name implies, it has historically been associated with Japan, as the first clinical case and isolation of JE were both recorded in Japan in 1871 and 1935, respectively [2,3]. The disease is caused by infection with JE virus (JEV), a mosquito-borne flavivirus, which is mainly transmitted by mosquitoes of the Culex vishnui subgroup [4,5]. The role of Culex tritaeniorhynchus, a member of the C. vishnui subgroup, as the most important JE vector was first reported in 1938 [6]. Culex pseudovishnui, another species in the C. vishnui subgroup, was later found to be involved in the transmission [4] and is currently regarded as a major vector of JE.
The application of next-generation sequencing (NGS) technology has enabled broad surveys of viral diversity [7] and thus has significantly sped up mosquito virome analysis in addition to disease surveillance. Although most of the detected viruses are considered insect specific and have no direct impact on public health [8], some of them may have the potential to serve as “natural competitors” to pathogens within mosquitoes. It has been hypothesized that the effects of these insect-specific viruses (ISVs), which can be used as biocontrols, include superinfection exclusion, antiviral immune response enhancement, and vertical transmission in mosquitoes [9]. Furthermore, a number of studies have reported the potential ability of ISV to modulate the transmission of viruses that are pathogenic to vertebrates by reducing the competence of the vector in transmitting arboviruses or modifying the host immune response to subsequent arbovirus infection [10,11,12]. As more evidence showed that the microbiome of the mosquito can alter the susceptibility to certain arboviruses [13], mosquito virome studies have significantly risen.
Emerging efforts for mosquito virome determination and characterization have been recorded in some parts of the world [14,15,16], and the virome mapping in Culex mosquitoes is no exception. To date, there are at least four studies examining the virome of C. tritaeniorhynchus mosquitoes [7,15,17,18]. All four studies have provinces in China as collection sites, whereas, in Japan, where JEV circulation has been going on for roughly 1.5 centuries, no virome study of the mosquito vectors has ever been done.
With a bulk RNA sequencing approach, such as NGS, the combination of virome profiling and disease surveillance became feasible; together, they may improve the understanding of the virus community in mosquitoes. This study comprehensively describes the virome profiles of two JE major vectors, C. tritaeniorhynchus and C. pseudovishnui, in comparison to that of another Culex species, Culex inatomii, collected in several prefectures in Japan. Using NGS, 27 viral genomes were uncovered, 15 of which were newly identified. These results revealed mosquito viral diversity, including the disease surveillance target, the JE.

2. Materials and Methods

2.1. Mosquito Collection and Species Identification

Mosquito collections were conducted in three prefectures of Japan: Ishikawa, Tottori, and Nagasaki, using CDC miniature light traps (The John W. Hock Co., Gainesville, FL, USA) with CO2, sweep nets, and/or aspirators. The locations selected for mosquito collection were based on area descriptions with a presumed link to JEV circulation, such as near pigpens (Isahaya City in Nagasaki Prefecture (32°49′ N, 130°03′ E)), rice field areas (several cities in Ishikawa Prefecture (37°29′ N, 136°74′ E)), and a bird migration hotspot (Yonago waterbirds sanctuary in Tottori Prefecture (35°26′ N, 133°17′ E)), as shown in Figure 1. The identification of Culex spp. was according to established identification keys [19]. After identification, mosquitoes were sorted by species and stored at −80 °C until further experiment. Only C. tritaeniorhynchus and C. pseudovishnui mosquitoes, which belong to the C. vishnui subgroup, and C. inatomii as the outgroup were processed in this study. The biological classification (taxonomic rank) of Culex mosquitoes used in this study is shown in Figure 2.

2.2. Metagenomic Analysis of RNA Virome in Culex Mosquitoes

A total of 652 mosquitoes, sorted into 13 pools of 21 to 25 mosquitoes, were processed for viral metagenomic analysis (Table 1) based on a previously described method [20]. Briefly, mosquito pools were homogenized, filtered, and then nuclease treated before RNA extraction using ISOGEN II (Nippon Gene, Toyama, Japan). Subsequently, the sequencing libraries were constructed and the amount of cDNAs was quantified by quantitative real-time polymerase chain reaction (qRT-PCR). All libraries were sequenced on an NGS (Illumina MiniSeq, San Diego, CA, USA). CLC Genomics Workbench version 11 (Qiagen, Aarhus C, Denmark) was used for processing nucleotide reading datasets. Sequences with >500 bp in length were chosen for the de novo assembly approach and the resulting datasets were compared with other nucleotide and amino acid sequences by BLASTn and BLASTx on the NCBI database. Standard database (nucleotide collection (nr/nt)) and highly similar sequences (megablast) optimization were selected for the BLASTn program, while standard (1) genetic code and non-redundant protein sequences (nr) as he target database were used for BLASTx program. Sequences with an 0.0 e-value and high query cover with other viral species were considered as known, while ‘unknown’ or potential viral sequences were checked for open reading frames (ORFs) in GENETYX ver. 13 software (Genetyx Corp., Tokyo, Japan) and/or NCBI ORFfinder, with parameters of minimal ORF 75 nt in length, standard genetic code, and ‘any sense codon’ for ORF start. HHpred interactive server (https://toolkit.tuebingen.mpg.de/tools/hhpred) for protein homology detection and NCBI Conserved Domain Database version 3.17 with an expected value threshold of 1 × 10−2 were also used for viral sequence detection, particularly for unknown short contigs with ORF(s). Subsequently, the short viral contigs with unassembled overlaps were merged, and translated proteins were aligned with closely related viruses as a reference. Any gaps found were filled by RT-PCR and Sanger sequencing using the ABI PRISM 3130 Genetic Analyzer (Thermofisher, Waltham, MA, USA). To verify the assembly results, reads were mapped back to NGS datasets on CLC Genomic Workbench and reference sequences.

2.3. Virus Isolation

A total of 2323 mosquitoes were employed for virus isolation (Table 1). The mosquito C6/36 cell line derived from Aedes albopictus Skuse (European Collection of Authenticated Cell Cultures, Darmstadt, Germany) was used for virus isolation as described previously [21]. Pools of mosquitoes (consisted of 21–25 individual females per pool) were homogenized, filter sterilized, and inoculated onto cells followed by incubation for 2 h to allow virus adsorption. After the addition of fresh medium, cell cultures were incubated at 28 °C with the same conditions for about 7 days. The final harvest of supernatants was conducted after at least two blind passages and stored at −80 °C until use. To identify the isolated virus, supernatants from the final harvest were once again subjected to NGS using the MiniSeq system as described above and subsequently confirmed by RT-PCR. Total reads were obtained and then assembled de novo in CLC Genomics Workbench (Supplementary Table S1) Time course assay in C6/36 mosquito cells was also performed to ensure virus isolation.

2.4. Pools Confirmation Using RT-PCR and Sequence Analysis of Viral RNA

The confirmation of virus-positive pools was done by RT-PCR. Viral RNA was extracted from cell culture supernatants using ISOGEN II. RT-PCR was conducted using TaKaRa PrimeScript One-Step RT-PCR Kit version 2 (Takara Bio, Kusatsu, Shiga, Japan) according to the manufacturer’s instructions using gene-specific primers designed from detected contigs (Supplementary Table S2). PCR amplifications were carried out in a total volume of 10 μL consisting of 0.5 μL of template (cDNA), 0.5 μM of forward and reverse primers, 0.5 μL of PrimeScript One Step Enzyme Mix, 5 μL 2× One Step Buffer, and Milli-Q ddH2O to supplement the system. Thermal cycling was initiated with incubation at 50 °C for 30 s and denaturation step at 94 °C for 5 min; followed by 35 cycles of 94 °C for 30 s, 53 °C for 30 s, and 72 °C for 30 s. The amplified products were visualized by agarose gel electrophoresis followed by capillary sequencing using ABI PRISM 3130 Genetic Analyzer.

2.5. Viral Genome Characterization and Phylogenetic Analyses

A number of contigs and assembled putative viral sequences from NGS were finally obtained (see Section 2.2.). ORFs and amino acid translations were determined by GENETYX ver. 13 software (Genetyx Corp.) and NCBI Orffinder (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/orffinder/) with a 300-nt minimum length. The genome structure was annotated using NCBI Conserved Domain Database version 3.17 (expected value threshold of 1 × 10−2) and the NCBI viral genome database as references.
To further characterize each viral isolate, the sequences were aligned by MAFFT version 7 (https://mafft.cbrc.jp/alignment/server/) employing G-INS-1 (slow) for a more accurate results and then trimmed by Gblocks version 0.91b (http://molevol.cmima.csic.es/castresana/Gblocks_server.html) [22] with a less stringent selection option (smaller final blocks, gap positions within the final blocks, and less strict flanking positions were allowed). Subsequently phylogenetic trees were constructed with MEGA version 6.06 using the maximum likelihood algorithm with 1000 bootstrap resampling. The substitution model and rates setting for the maximum likelihood tree were determined in advance by the best DNA/protein model searching tool. Viral genomes with less than or equal to 83% nucleotide identity to previously described viruses were identified as novel. The nucleotide sequence accession numbers are provided in Table 2 and Supplementary Table S3.

3. Results

3.1. RNA Virome Analysis by Metagenomic Sequencing

Viral metagenomic sequencing was performed on 13 pools of Culex mosquitoes, including both the well-known JE vector group (C. vishnui subgroup) and a unknown JE vector group (C. inatomii). The viromes of C. tritaeniorhynchus mosquitoes collected from three locations were almost identical, with seven out of nine viral families present across locations (Figure 3A). Regardless of the total number of samples used, it is interesting how virus families, such as Totiviridae, Partitiviridae, Flaviviridae, Iflaviridae, and Xinmoviridae, appeared concurrently, suggesting a robust relationship of the virome with the host, C. tritaeniorhynchus.
In comparison to the C. tritaeniorhynchus virome, the result of the C. pseudovishnui and C. inatomii viromes seemed to clarify the impact of the host taxon in the virome structure, as the three mosquito species have a different set of virus families (Figure 3B). Despite the differences, C. tritaeniorhynchus shared some virus families with the closely related species C. pseudovishnui, three families in total, including an unclassified virus family.
From a total of 27 identified viruses, 18, 11, and 3 viruses were identified from C. tritaeniorhynchus, C. pseudovishnui, and C. inatomii, respectively, suggesting viral abundance variation in each mosquito species (Figure 4, Table 3). Eight viruses were detected from all pools of C. tritaeniorhynchus mosquitoes collected from three prefectures, belonging to Partitiviridae, Flaviviridae, Totiviridae, and unclassified virus families, respectively (Figure 4, in the center of overlapping areas). The result showed corresponding patterns with the previous result of the virus family ratio, predisposing mosquito taxon over geographical location as the virome-shaping cause.
Culex pseudovishnui collected from one collection site in Nagasaki also displayed an interesting pattern in terms of shared viruses, as it shared three viruses with C. tritaeniorhynchus (Table 3). Being the outgroup, only three viruses were detected in C. inatomii, those belonging to Totiviridae and unclassified families.
Taken together, the viruses identified in this study comprised double-stranded RNA (dsRNA) viruses (Totiviridae, Partitiviridae, and Chrysoviridae families), positive-sense single-stranded RNA [(+)ssRNA] viruses (Flaviviridae and Iflaviridae families, and unclassified groups), and negative-sense ssRNA [(−)ssRNA] viruses (Xinmoviridae and Rhabdoviridae families, and unclassified groups). Further details of each virus are described below.

3.2. Dataset Breakdown of Mosquito RNA Viruses

The metagenome sequencing method uncovered a wide range of viruses from a total of 652 mosquitoes grouped into 13 pools (Table 1). In total, BLAST analyses of NGS viral reads identified 27 RNA viral genomes, 15 of which were newly described (Table 2, Supplementary Table S3). All characterized viruses in this study were phylogenetically within an order and related or belonged to viral families listed in the International Committee on Taxonomy of Viruses.

3.2.1. dsRNA Viruses

Totiviridae

Totiviruses are known as a group of dsRNA viruses that infect fungi, protozoa, or invertebrates [23]. In this study, most mosquito pools have a totivirus sequence present. Culex vishnui subgroup mosquitoes have the exact same totivirus beyond collection sites (designated as Culex vishnui subgroup totivirus, CvsTV), whereas C. inatomii mosquitoes have a different totivirus (designated as Culex inatomii totivirus, CiTV) (Figure 5A). CvsTV was highly distributed in the pools tested, and it was present in Culex vishnui subgroup pools (Table 3). CvsTV presented a close relationship with dsRNA virus environmental sample (Figure 5B), whereas CiTV formed a cluster with Australian Anopheles totivirus.

Chrysoviridae

Chrysoviruses are known to infect fungi [24], but recent studies have revealed an association with insects. Hubei chryso-like virus 1 (HbCl-1) was described previously in mosquitoes [25]. In this study, a new strain of this viral sequence was found in most pools of C. tritaeniorhynchus originating from Ishikawa (Figure 5).

Partitiviridae

Partitiviruses are widely known to be associated with plants and fungi; however, recent studies have suggested that their known host range might be expanded, and some members in this genus were discovered in arthropods [26,27]. In this study, two partitivirus sequences were found in C. tritaeniorhynchus: New strains of Hubei partiti-like virus 22 and a novel partitivirus designated as Culex tritaeniorhynchus partitivirus (CtPV). Culex pseudovishnui has a different novel partitivirus named Culex pseudovishnui partitivirus (CpPV). CtPV formed a cluster with Hubei partiti-like virus 19, Hubei partiti-like virus 18, and Galbut virus, whereas CpPV displayed a correlation with Hubei partiti-like virus 56 and Changjiang partiti-like virus 1 (Figure 6).

3.2.2. (+)ssRNA viruses

The majority of viruses identified in this study were (+)ssRNA viruses. As many as 15 (+)ssRNA viral genomes were identified. Some were previously described as ISVs, some were novels, and one was a mosquito-borne arbovirus.

Flaviviridae

This virus family has humans and mammals asits natural host, which makes virus discovery of this family noteworthy. Flaviviridae members found by NGS are mosquito flavivirus, Culex tritaeniorhynchus flavi-like virus (CtFLV), and JEV. A strain of mosquito flavivirus, described as Yamadai flavivirus, was found in two C. pseudovishnui pools. CtFLV belonged to an unclassified flavivirus clade together with viruses discovered from arthropods: Gamboa mosquito virus, Macrosiphum euphorbiae virus 1, Shayang fly virus 4, and Shuangao lacewing virus 2 (Figure 7), based on the NS5 region.
JEV strain 17CxIT-I4-D31 was successfully characterized as genotype I based on envelope region, and it was isolated from a C. tritaeniorhynchus sample from Ishikawa. Based on the evolutionary tree (Figure 8), JEV identified in this study was closely related to isolates from relatively distant areas, such as JEVs detected in Miyazaki Prefecture in 2009, Yamaguchi Prefecture in 2013, and Zhejiang, China in 2007. As widely known, there are five genotypes of JEV in the world; currently, genotype I is the predominant genotype [28,29].

Negevirus-Related Viruses

Negevirus was recently described as a diverse group of ISVs isolated from mosquitoes and phlebotomine sandflies [30]. In this study, two novel negevirus-like sequences were found in C. vishnui subgroup mosquitoes. A novel virus was identified from all pools of C. tritaeniorhynchus and designated as Culex tritaeniorhynchus negev-like virus (CtNLV). Another novel virus was identified from C. pseudovishnui and designated as Culex pseudovishnui negev-like virus (CpNLV). Phylogenetic analysis revealed that CtNLV formed a cluster with Mill Lade virus and Aedes camptorhynchus negev-like, whereas CpNLV formed a cluster with Yongsan negev-like virus (Figure 9).

Tymovirales

Tymovirales is an order of non-segmented (+)ssRNA viruses, which currently includes five families. Most members of the order infect plants, but a few species are from plant pathogenic fungi [31]. Generally, members of the family Tymoviridae have plants as their natural hosts with insects as a mechanical vector. Studies have recently discovered tymo-related viruses, such as Tarnsjo virus and Guadeloupe Culex-tymo-like virus in mosquitoes, thus implying the relationship in nature. A tymo-like virus found in this study was detected in two pools of C. pseudovishnui and designated as Culex pseudovishnui tymo-like virus (CpTLV). CpTLV formed a cluster with two recently identified viruses: Tarnsjo virus and Guadeloupe Culex-tymo-like virus (Figure 9).

Luteoviridae-Related Viruses

Luteoviruses are a family of viruses commonly known to include many economically important plant pathogens [32]. Recently, luteo-like viruses (related but phylogenetically divergent compared to recognized members) have been discovered in diverse invertebrates, such as dragonflies, spiders, and mosquitoes, indicating that this group of viruses has wider range than previously known [27,33]. In this study, three known luteo-like viruses were identified in C. tritaeniorhynchus and/or C. pseudovishnui: Wenzhou tombus-like virus 11 (WTLV11), Hubei mosquito virus 2 (HMV2), and Hubei mosquito virus 4 (HMV4) (Figure 9). A novel luteo-like virus designated as C. inatomii luteo-like virus (CiLLV) was identified in C. inatomii. Based on the RdRp region, CiLLV was closely related to HMV2 (Figure 9); likewise, the nucleotide sequence of segment 1 and segment 2 displayed similar results (Supplementary Figure S1).

Sobemovirus-Related Viruses

For sobemoviruses, plants are their sole natural host. However, bulk RNA sequencing has increased the scope of their host range, with some viruses that belong to this group found in vertebrates and arthropods [33]. In this study, two known sobemo-like viruses, Wenzhou sobemo-like virus 3 (WSLV3) and bat sobemovirus, were found in all pools of C. vishnui subgroup mosquitoes but not in the C. inatomii library.
WSLV3 was reported to have 2841 bp nucleotides [33], shorter than the average sobemo-like viruses that have 4000 to 5000 bp nucleotides. In contrast, only the coat protein of bat sobemovirus was submitted by Dacheux et al. [34] without any other segments nor the polymerase region. Furthermore, multiple appearances of these two viruses in the same mosquito pools have inferred the possibility of the two being the same virus (as illustrated in Figure 10), although more studies are required to confirm this finding.

Iflaviridae

Iflaviruses are widely known as a group of picorna-like viruses that infect arthropods [35]. In this study, three iflaviruses were found: A novel iflavirus designated as Isahaya Culex iflavirus (ICIFV), which was detected in most pools of C. tritaeniorhynchus and C. pseudovishnui collected in Nagasaki; another strain of Yongsan iflavirus; and a novel iflavirus designated Yonago Culex iflavirus (YCIFV), which was detected from C. tritaeniorhynchus pools from Tottori. The phylogenetic tree displayed ICIFV as being in the same cluster with Bee iflavirus 1 and Wuhan insect virus 13, whereas YCIFV formed a cluster together with Kinkell virus and Wuhan fly virus 4 China (Figure 11).

3.2.3. (−)ssRNA Viruses

Five putative (−)ssRNA viruses, including major taxonomic categories, were discovered. Among these, one virus is a known mosquito virus, whereas the remaining four viruses are newly described in this study.

Xinmoviridae

The recently updated classification of Mononegavirales taxon has a new family named Xinmoviridae, established for the floating genus Anphevirus [36,37]. In this study, a new member of the virus family, designated as C. tritaeniorhynchus anphevirus (CtAV), was found in most pools of C. tritaeniorhynchus. Phylogenetic analysis showed a close relationship of CtAV with Guadeloupe mosquito mononega-like virus and Culex mononega-like virus 2 (Figure 12).

Rhabdoviridae

Rhabdoviruses are a diverse set of (−)ssRNA viruses known to infect both animals and plants [38], with frequent host switching during their evolution history [39]. Two rhabdoviruses were detected in this study: Culex tritaeniorhynchus rhabdovirus was found in C. tritaeniorhynchus only from Ishikawa and Tottori, whereas a novel virus designated as Culex pseudovishnui rhabdo-like virus (CpRLV) was solely identified in C. pseudovishnui. CpRLV displayed a close relationship with Tongilchon virus (Figure 12).

Bunyavirales

Bunyavirales is an order of segmented (−)ssRNA viruses, which currently includes more than 10 families [40,41]. In this study, sequences related to bunyavirus L and M genes were detected from two C. pseudovishnui pools. The new bunya-like virus was tentatively named C. pseudovishnui bunya-like virus (CpBLV). The phylogenetic tree based on the RdRp sequences encoded by the L segment revealed CpBLV as being in the same clade as RdRp of the Bunyaviridae environmental sample and Culex bunyavirus 2, which possibly belong to the Phenuiviridae family (Figure 12). Phylogenetic analysis of M and S segments of CpBLV were also provided (Supplementary Figure S2)

Aspiviridae-Related Viruses

In this study, a virus related to the Aspiviridae (formerly called Ophioviridae) family, was identified, which was designated Culex tritaeniorhynchus aspi-like virus (CtALV). CtALV was detected in a C. tritaeniorhynchus pool collected in Tottori. The phylogenetic tree showed that CtALV is closely related to Wilkie ophio-like virus 1 detected from A. camptorhynchus mosquitoes in Australia [33] and ophiovirus-related viruses of plant pathogenic fungi [42] (Figure 12).

3.3. Attempt to Isolate Viruses Using a Mosquito Cell Line

C6/36 cells derived from A. albopictus were used in this study to isolate the identified viruses (Table 1) and to disclose any latent viruses undetected by NGS. After three blind passages to propagate viruses, supernatants were harvested before being subsequently examined with time course assay for validation of the virus isolation. Through this method, one medically important virus and five insect-related viruses were successfully isolated by C6/36 cells: JEV, mosquito flavivirus, Culex flavivirus (CxFV), and a novel iflavirus called YCIFV.
Although the cell line used was mosquitoes, the isolated viruses were not limited to ISVs. Moreover, the fact that not every ISV was isolatable illustrated the cell line had variable sensitivity to virus infection. Lineages of the C6/36 cell line have been widely used to study the relationship between arboviruses and mosquito vectors [43] and they have been shown to be susceptible to a wide range of arboviruses, partially due to the lack of a functional RNA interference response [44].
Three blind passages of virus isolation in mosquito cell lines amplified two viruses not detected by NGS: CxFV and mosquito flavivirus in the Nagasaki group’s C. tritaeniorhynchus. The fact that those undetected viruses belong to the Flavivirus genus, which contains many medically important viruses, highlighted the sense that, in arbovirus surveillance, the method of virus amplification through mosquito cell lines should not be overlooked.

3.4. Pools Confirmation and Additional Individual Screening

Using the NGS result as the lead, primers were designed and RNA viruses in pools were confirmed by RT-PCR; thus, the virome of each mosquito species and collection sites can be analyzed and putative seasonal changes of virome composition can also be portrayed (Supplementary Figure S3). There were two patterns observed in seasonal virome composition: Continuous presence (detected in every month during mosquito emergence of early summer to late autumn) and occasional presence (irregular detection). To establish the proportion of highly prevalent viruses in mosquitoes, additional data of individual infection rates were provided. Results showed percentages of infectivity that ranged from 4% to 58% (Supplementary Table S4), suggesting the active circulation of viruses infecting most of the mosquitoes in the population. Nevertheless, further studies over a longer period of time (3–5 years) would be required to confirm the seasonal changes.

4. Discussion

It is estimated that approximately 67,900 JE cases typically occur annually in 24 JE-endemic countries [1], including Japan. A survey on the natural infection of JEV in humans and horses indicated that JE remains prevalent in Japan [45,46] and C. tritaeniorhynchus belonging to the vishnui subgroup is the main vector in Japan [47]. This study successfully combined mosquito-borne arbovirus surveillance with virome profiling in major JE vectors collected in three prefectures in Japan.
The virome of three Culex mosquito species, C. tritaeniorhynchus and C. pseudovishnui renowned as important vectors for JEV in Asia and C. inatomii as the non-vector species, was analyzed to better understand arbovirus–ISV interaction. The abundance and diversity of possible ISVs constituted a significant component of Culex mosquito viromes; for example, the virome of Culex spp. in Zambezi, Mozambique was dominated (88.5%) by them [48], whereas a study of the C. tritaeniorhynchus virome in Hubei, China also showed a conformable predominance (88%) of insect viruses [18]. In line with the reckonings, the result displayed a practically similar constitution, as much as 84%, 82%, and 100% of viruses were insect viruses, or had a high correlation with the mosquito virome previously identified, in C. tritaeniorhynchus, C. pseudovishnui, and C. inatomii, respectively (Supplementary Figure S4), implying a consistent virome composition in Culex mosquitoes, with slight variations possibly associated with the number of samples used.
Metagenomic sequencing results suggested a distinct virus distribution depending on the mosquito taxonomy. Here, it was shown that two mosquitoes belonging to the same subgroup, C. tritaeniorhynchus and C. pseudovishnui, shared a number of virus families (Totiviridae, Partitiviridae, Flaviviridae, Iflaviridae, and Rhabdoviridae), whereas another Culex belonging to a different subgenus, C. inatomii, shared only one virus with mosquitoes collected in the same area. This finding corresponds to a study comparing viral diversity in some species of mosquitoes, which showed correlation between mosquito species and their virome structures [49].
Recurring detection of certain predominant viruses from one mosquito species, C. tritaeniorhynchus, originated from various collection sites, and time has corroborated the taxon impact on virome shaping even more. The absence of these predominant viruses in other species could imply disparity in host susceptibility, although this would need further clarification.
Not only does this study display a pattern of host taxon as one of the virome-shaping factors, but it also reveals another point—that some viruses could be detected in discrete Culex mosquito species inhabiting the same area. This phenomenon was vividly demonstrated by detection of HMV4, which was detected from the two mosquito species (C. tritaeniorhynchus and C. inatomii) collected from the same sites in Tottori. In some way, particular viruses present in this study were also discovered in geographically distant localities, such as HMV2 and WSLV3, which were already detected in C. tritaeniorhynchus in China [17]. In the beginning, virus predominance in Culex sp. mosquitoes has long been speculated, before nucleotide identities (99% in both viruses) remarked an untypical limited genetic transformation despite the geographical distance. Again, this correlated with a previous finding that some viruses infected hosts over a wide geographical area and were possibly introduced by windblown, cyclones, or humans rather than the spread of mosquitoes across continents over time [33,50,51,52]. Clearly, however, more extensive studies need to be undertaken to determine if these virome patterns are consistent beyond mosquito species and geographic regions.
Aside from the patterns observed, this study also discusses the viral association with hosts and possible interrelation with other coexisting viruses. The arbovirus transmitted by mosquitoes detected in the study was JEV, which was identified in C. tritaeniorhynchus from Ishikawa Prefecture. No other virus was considered as arbovirus based on their closest related viruses or phylogenetic tree. An interesting aspect exhibited from virus isolation might be that another flavivirus, CxFV, was not detected in the mosquito homogenate before virus isolation. Nagasaki Prefecture, where CxFV isolates were collected, is a reputable JE-endemic area in Japan, as JEV detection is frequent in annual disease surveillance of amplifiers and vectors. JEV was not isolated nor identified in 2017 Nagasaki samples, thus making it difficult to decipher the interrelation between the flaviviruses. Although an in vitro study by Kuwata et al. [53] implied that JEV superinfection might be harmful to CxFV-infected C. tritaeniorhynchus mosquitoes, this study was unable to demonstrate the correlation between CxFV occurrence and the detection or otherwise of JEV.
The only collection site where JEV was identified is Monzenmachi-toge, Ishikawa Prefecture, where the landscape is a vast rice field surrounded by forest with the river and sea in close proximity. The result was supported by annual disease surveillance conducted in September 2017 by the Japan National Institute of Infectious Disease (2017) [54], which has indicated that 50% to 60% of swine sera had JE antibodies, although, in the previous 2 years (2015 and 2016), the range was only between 0% and 10%. An increase of JEV antibodies in swine sera contributed greatly to JEV detection in mosquito, as swine plays a major role as a JEV amplifier [55].
Phylogenetic analysis of the complete genome showed that the isolated JEV strain 17CxIT-I4-D31 was closely related to JEV isolates detected in Yamaguchi in 2013 rather than those from Ishikawa previously or other surrounding areas, suggesting that JEV may be transported from west Japan into central Japan. The phenomenon of this virus spread has also occurred in Getah virus (genus Alphavirus, family Togaviridae). A study reported the virus being transported in the same direction from west Japan into central Japan, with possible routes of transportation, such as bird migration and windblown mosquitoes [56].
Detection of JEV in Ishikawa, where other ISVs are circulating simultaneously (Supplementary Figure S3), could signify the natural occurrence of coexistence of the virome. However, this must be interpreted carefully, as this study used mosquito pools instead of individual mosquitoes.
Although this work used mosquitos as the sample, surprisingly, viral hosts other than mosquitos appeared from the result. As mentioned previously, where identified viruses were either insect specific or human arbovirus, another virus host, which is a disease vector, emerged. Bat sobemovirus (named Bordeaux sobemovirus) was first identified from pooled organs (lungs, liver, and brain) of Pipistrellus pipistrellus, an insectivorous bat, in France in 2009 [34]. Yet, it found its way into the C. vishnui subgroup mosquitoes in Japan. This finding is not novel, as other studies similarly found the linkage between the bat and arthropod through the virome [57,58]. However, it should be noted that the bat might have a key role in the sylvatic cycle of flaviviruses [59]; additionally, there has been evidence of JEV presence in the bat [60,61]. Therefore, substantial appearances of bat sobemovirus in all libraries of C. vishnui subgroup mosquitoes may strengthen the possible relationship between bat and JE mosquito vectors.
Nevertheless, the bat sobemovirus nucleotide sequence identified by Dacheux et al. (2014) [34] was a viral coat protein with no polymerase or other segments. In this study, bat sobemovirus sequences were found in pools where WSLV3 was also detected. More studies are required to clarify the connection between the two viral sequences; however, the trace this study gave has indicated the possibility of the two being the same virus.
Importantly, this mosquito virome study was carried out using metagenomic sequencing in combination with mosquito cell line virus isolation. Efforts have been put into metagenomic analysis using NGS for disease surveillance, as the method allows the simultaneous identification of viruses from a single mosquito in a single reaction [62]. However, NGS methods have disadvantages compared to other molecular methods of virus detection, as they are less sensitive than qRT-PCR for the detection of samples with low virus titers, particularly when using mosquitoes as samples.
The combination of high-throughput sequencing and virus isolation in this study was meant to retrieve the complete mosquito virome, including underlying viruses that exist at very low copy numbers, whereupon the results can be interpreted. C6/36 cells were selected, as they are easy to maintain and highly permissive to numerous arboviruses [63], although the limitation is that they may fail to accurately model mosquito–arbovirus interactions [44], and the output given for this virome identification is of importance, as they successfully uncovered viruses undetected by preceding NGS. The detection of two insect-specific flaviviruses (ISF), CxFV and MFV, in only the third supernatants should not be overlooked. This is because they belong to a medically important virus genus and thus may imply potentially harmful viruses persisting at low copy numbers, undetected by NGS.
Finally, complete analysis was performed on the virome profile of major JE vectors in Japan, including an arbovirus, ISVs, and other viruses. This is the first report on the Culex mosquito virome study in Japan. In addition to JEV and several known ISVs, other novel viral sequences were identified from these mosquitoes. Notwithstanding, further studies will be required to confirm the real host of the identified viruses as the metatranscriptomic studies may also include non-mosquito RNA. By studying one of the medically important disease vectors, this study offers important insights into the ecology and evolution of mosquito viromes consisting of both arbovirus and insect viruses.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4915/12/3/264/s1, Figure S1: (A) Phylogenetic analysis of Luteo-sobemo cluster using nucleotide sequences of segment 1. (B) Phylogenetic analysis of Luteo-sobemo cluster using nucleotide sequences of segment 2. Figure S2: (A) Phylogenetic analysis of Bunyavirus-like cluster using nucleotide sequences of M segment. (B) Phylogenetic analysis of Bunyavirus-like cluster using nucleotide sequences of S segment. Figure S3: (A) Seasonal changes in virome composition of C. tritaeniorhynchus mosquitoes in Tottori in 2017. (B) Seasonal changes in virome composition of C. tritaeniorhynchus mosquitoes in Ishikawa in 2017. Figure S4: Virome composition based on viruses’ presumptive host. Table S1: Percentage of viral reads from NGS. Table S2: Primers used in the study. Table S3: Virome identification result (additional viral segments). Table S4: Individual infection rate of several viruses.

Author Contributions

Conceptualization, A.N.F., D.K., H.I. and K.S.; Methodology, A.N.F., D.K., H.I., Y.H., K.I.; Validation, A.N.F., D.K., H.I., K.S., K.M. (Kozue Miura), K.M. (Katsunori Murota) and K.H.; Formal analysis, A.N.F., D.K. and H.I.; Investigation, A.N.F., D.K., M.A.-B., K.M. (Kozue Miura), K.M. (Katsunori Murota) Y.H., K.F., S.S., K.I., M.W., Y.T. and K.S.K.; Resource, M.W., K.S.K., H.I., N.M. and K.S.; Data curation, A.N.F., D.K. and H.I.; Writing—original draft, A.N.F.; Writing—review and editing: A.N.F., D.K., M.A.-B., H.I. and K.S.; Visualization, A.N.F., D.K., H.I., K.M. (Kozue Miura) and K.H.; Supervision, A.N.F., H.I., K.S., K.M. (Katsunori Murota), K.M. (Kozue Miura) and K.H.; Project Administration, A.N.F, D.K., H.I., N.M., K.M. (Kozue Miura), K.H. and K.S.; Funding acquisition, D.K., H.I., K.S. and K.H. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Regulatory Science Research from Ministry of Health, Labor and Welfare, the Research Program on Emerging and Re-emerging Infectious Diseases from Japan Agency for Medical Research and development (AMED), and JSPS KAKENHI Grant Numbers JP18K19220 and JP18H02856.

Acknowledgments

The authors would like to thank Enago (www.enago.jp) for the English language review.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Campbell, G.; Hills, S.; Fischer, M.; Jacobson, J.; Hoke, C.; Hombach, J.; Marfin, A.; Solomon, T.; Tsai, T.; Tsui, V.; et al. Estimated global incidence of Japanese encephalitis: A systematic review. Bull. World Health Organ. 2011, 89, 766–774. [Google Scholar] [CrossRef] [PubMed]
  2. Erlanger, T.E.; Weiss, S.; Keiser, J.; Utzinger, J.; Wiedenmayer, K. Past, present, and future of Japanese encephalitis. Emerg. Infect. Dis. 2009, 15, 1–7. [Google Scholar] [CrossRef] [PubMed]
  3. Solomon, T.; Dungc, N.M.; Kneend, R.; Gainsboroughc, M.; Vaughnf, D.W.; Khanhc, V.T. Japanese encephalitis. J. Neurol. Neurosurg. Psychiatry 2000, 68, 405–415. [Google Scholar] [CrossRef] [PubMed]
  4. Mourya, D.T.; Mishra, A.C.; Soman, R.S. Transmission of Japanese encephalitis virus in Culex pseudovishnui & C. tritaeniorhynchus mosquitoes. Indian J. Med. Res. 1991, 93, 250–252. [Google Scholar]
  5. Van den Hurk, A.F.; Ritchie, S.A.; Mackenzie, J.S. Ecology and geographical expansion of Japanese encephalitis virus. Annu. Rev. Entomol. 2009, 54, 17–35. [Google Scholar] [CrossRef] [Green Version]
  6. Buescher, E.L.; Scherer, W.F.; Rosenberg, M.Z.; Gresser, I.; Hardy, J.L.; Bullock, H.R. Ecologic studies of Japanese encephalitis virus in Japan. II. Mosquito infection. Am. J. Trop. Med. Hyg. 1959, 8, 651–664. [Google Scholar] [CrossRef]
  7. Atoni, E.; Wang, Y.; Karungu, S.; Waruhiu, C.; Zohaib, A.; Obanda, V.; Agwanda, B.; Mutua, M.; Xia, H.; Yuan, Z. Metagenomic virome analysis of Culex mosquitoes from Kenya and China. Viruses 2018, 10, 30. [Google Scholar] [CrossRef] [Green Version]
  8. Hall, R.A.; Bielefeldt-Ohmann, H.; McLean, B.J.; O’Brien, C.A.; Colmant, A.M.; Piyasena, T.B.; Harrison, J.J.; Newton, N.D.; Barnard, R.T.; Prow, N.A.; et al. Commensal viruses of mosquitoes: Host restriction, transmission, and interaction with arboviral pathogens. Evol. Bioinform. Online 2017, 12, 35–44. [Google Scholar] [CrossRef]
  9. Öhlund, P.; Lundén, H.; Blomström, A.L. Insect-specific virus evolution and potential effects on vector competence. Virus Genes 2019, 55, 127–137. [Google Scholar] [CrossRef] [Green Version]
  10. Hall-Mendelin, S.; McLean, B.J.; Bielefeldt-Ohmann, H.; Hobson-Peters, J.; Hall, R.A.; van den Hurk, A.F. The insect-specific Palm Creek virus modulates West Nile virus infection in and transmission by Australian mosquitoes. Parasit. Vectors 2016, 9, 414. [Google Scholar] [CrossRef] [Green Version]
  11. Vasilakis, N.; Tesh, R.B. Insect-specific viruses and their potential impact on arbovirus transmission. Curr. Opin. Virol. 2015, 15, 69–74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Bolling, B.G.; Weaver, S.C.; Tesh, R.B.; Vasilakis, N. Insect-specific virus discovery: Significance for the arbovirus community. Viruses 2015, 7, 4911–4928. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Jupatanakul, N.; Sim, S.; Dimopoulos, G. The insect microbiome modulates vector competence for arboviruses. Viruses 2014, 6, 4294–4313. [Google Scholar] [CrossRef] [PubMed]
  14. Zakrzewski, M.; Rašić, G.; Darbro, J.; Krause, L.; Poo, Y.S.; Filipović, I.; Parry, R.; Asgari, S.; Devine, G.; Suhrbier, A. Mapping the virome in wild-caught Aedes aegypti from Cairns and Bangkok. Sci. Rep. 2018, 8, 4690. [Google Scholar] [CrossRef] [Green Version]
  15. Xiao, P.; Li, C.; Zhang, Y.; Han, J.; Guo, X.; Xie, L.; Tian, M.; Li, Y.; Wang, M.; Liu, H.; et al. Metagenomic Sequencing from mosquitoes in China reveals a variety of insect and human viruses. Front. Cell. Infect. Microbiol. 2018, 8, 364. [Google Scholar] [CrossRef]
  16. Cholleti, H.; Hayer, J.; Fafetine, J.; Berg, M.; Blomström, A.L. Genetic characterization of a novel picorna-like virus in Culex spp. mosquitoes from Mozambique. Virol J. 2018, 15, 71. [Google Scholar] [CrossRef] [Green Version]
  17. Xia, H.; Wang, Y.; Shi, C.; Atoni, E.; Zhao, L.; Yuan, Z. Comparative metagenomic profiling of viromes associated with four common mosquito species in China. Virol. Sin. 2018, 33, 59–66. [Google Scholar] [CrossRef]
  18. Shi, C.; Liu, Y.; Hu, X.; Xiong, J.; Zhang, B.; Yuan, Z. A metagenomic survey of viral abundance and diversity in mosquitoes from Hubei province. PLoS ONE 2015, 10, e0129845. [Google Scholar] [CrossRef]
  19. Tanaka, K.; Mizusawa, K.; Saugstad, E.S. A revision of the adult and larval mosquitoes of Japan (including the Ryukyu Archipelago and the Ogasawara Islands) and Korea (Diptera: Culicidae). Contrib. Am. Entomol. Inst. (Ann Arbor) 1979, 16, 1–987. [Google Scholar]
  20. Kobayashi, D.; Murota, K.; Itokawa, K.; Ejiri, H.; Amoa-Bosompem, M.; Faizah, A.N.; Watanabe, M.; Maekawa, Y.; Hayashi, T.; Noda, S.; et al. RNA virome analysis of questing ticks from Hokuriku District, Japan, and the evolutionary dynamics of tick-borne phleboviruses. Ticks Tick Borne Dis. 2020, 11. [Google Scholar] [CrossRef]
  21. Hoshino, K.; Isawa, H.; Tsuda, Y.; Sawabe, K.; Kobayashi, M. Isolation and characterization of a new insect flavivirus from Aedes albopictus and Aedes flavopictus mosquitoes in Japan. Virology 2009, 391, 119–129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Castresana, J. Selection of Conserved Blocks from Multiple Alignments for Their Use in Phylogenetic Analysis. Mol. Biol. Evol. 2000, 17, 540–552. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Ghabrial, S.A.; Nibert, M.L. Victorivirus, a new genus of fungal viruses in the family Totiviridae. Arch. Virol. 2009, 154, 373–379. [Google Scholar] [CrossRef] [PubMed]
  24. Ghabrial, S.A.; Castón, J.R.; Coutts, R.H.A.; Hillman, B.I.; Jiang, D.; Kim, D.H.; Moriyama, H.; Ictv report consortium. ICTV Virus Taxonomy Profile: Chrysoviridae. J. Gen. Virol. 2018, 99, 19–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Shi, M.; Lin, X.D.; Tian, J.H.; Chen, L.J.; Chen, X.; Li, C.X.; Qin, X.C.; Li, J.; Cao, J.P.; Eden, J.S.; et al. Redefining the invertebrate RNA virosphere. Nature 2016, 540, 539–543. [Google Scholar] [CrossRef]
  26. Nibert, M.L.; Ghabrial, S.A.; Maiss, E.; Lesker, T.; Vainio, E.J.; Jiang, D.; Suzuki, N. Taxonomic reorganization of family Partitiviridae and other recent progress in partitivirus research. Virus Res. 2014, 188, 128–141. [Google Scholar] [CrossRef]
  27. Pettersson, J.H.; Shi, M.; Bohlin, J.; Eldholm, V.; Brynildsrud, O.B.; Paulsen, K.M.; Andreassen, Å.; Holmes, E.C. Characterizing the virome of Ixodes ricinus ticks from northern Europe. Sci. Rep. 2017, 7, 10870. [Google Scholar] [CrossRef] [Green Version]
  28. Solomon, T.; Ni, H.; Beasley, D.W.C.; Ekkelenkamp, M.; Cardosa, M.J.; Barrett, A.D.T. Origin and evolution of Japanese encephalitis virus in southeast asia. J. Virol. 2003, 77, 3091–3098. [Google Scholar] [CrossRef] [Green Version]
  29. Ma, S.P.; Yoshida, Y.; Makino, Y.; Tadano, M.; Ono, T.; Ogawa, M. Short report: A major genotype of Japanese encephalitis virus currently circulating in Japan. Am. J. Trop. Med. Hyg. 2003, 69, 151–154. [Google Scholar] [CrossRef]
  30. Nunes, M.R.T.; Contreras-Gutierrez, M.A.; Guzman, H.; Martins, L.C.; Barbirato, M.F.; Savit, C.; Balta, V.; Uribe, S.; Vivero, R.; Suaza, J.D.; et al. Genetic characterization, molecular epidemiology, and phylogenetic relationships of insect-specific viruses in the taxon Negevirus. Virology 2017, 504, 152–167. [Google Scholar] [CrossRef]
  31. Adams, M.J.; Kreuze, J.F.; Martelli, G.P. Order Tymovirales. In Virus taxonomy, Ninth report of the International Committee on Taxonomy of Viruses; Andrew, M.Q.K., Adams, M.J., Carstens, E.B., Lefkowitz, E.J., Eds.; Elsevier Academic Press: London, UK, 2012; pp. 901–903. [Google Scholar] [CrossRef]
  32. Ali, M.; Hameed, S.; Tahir, M. Luteovirus: Insights into pathogenicity. Arch. Virol. 2014, 159, 2853–2860. [Google Scholar] [CrossRef] [PubMed]
  33. Shi, M.; Neville, P.; Nicholson, J.; Eden, J.S.; Imrie, A.; Holmes, E.C. High-resolution metatranscriptomics reveals the ecological dynamics of mosquito-associated RNA viruses in Western Australia. J. Virol. 2017, 91, e00680-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Dacheux, L.; Cervantes-Gonzalez, M.; Guigon, G.; Thiberge, J.M.; Vandenbogaert, M.; Maufrais, C.; Caro, V.; Bourhy, H. A preliminary study of viral metagenomics of French bat species in contact with humans: Identification of new mammalian viruses. PLoS ONE 2014, 9, e87194. [Google Scholar] [CrossRef] [PubMed]
  35. Chen, Y.P.; Nakashima, N.; Christian, P.D.; Bakonyi, T.; Bonning, B.C.; Valles, S.M.; Lightner, D. Family Iflaviridae. In Virus Taxonomy, Ninth Report of the International Committee on Taxonomy of Viruses; Andrew, M.Q.K., Adams, M.J., Carstens, E.B., Lefkowitz, E.J., Eds.; Elsevier Academic Press: London, UK, 2012; pp. 846–849. [Google Scholar] [CrossRef]
  36. Maes, P.; Amarasinghe, G.K.; Ayllón, M.A.; Basler, C.F.; Bavari, S.; Blasdell, K.R.; Briese, T.; Brown, P.A.; Bukreyev, A.; Balkema-Buschmann, A.; et al. Taxonomy of the order Mononegavirales: Second update 2018. Arch. Virol. 2019, 164, 1233–1244. [Google Scholar] [CrossRef] [Green Version]
  37. Amarasinghe, G.K.; Aréchiga Ceballos, N.G.; Banyard, A.C.; Basler, C.F.; Bavari, S.; Bennett, A.J.; Blasdell, K.R.; Briese, T.; Bukreyev, A.; Caì, Y.; et al. Taxonomy of the order Mononegavirales: Update 2018. Arch. Virol. 2018, 163, 2283–2294. [Google Scholar] [CrossRef] [Green Version]
  38. Bourhy, H.; Cowley, J.A.; Larrous, F.; Holmes, E.C.; Walker, P.J. Phylogenetic relationships among rhabdoviruses inferred using the L polymerase gene. J. Gen. Virol. 2005, 86, 2849–2858. [Google Scholar] [CrossRef]
  39. Longdon, B.; Murray, G.G.; Palmer, W.J.; Day, J.P.; Parker, D.J.; Welch, J.J.; Obbard, D.J.; Jiggins, F.M. The evolution, diversity, and host associations of rhabdoviruses. Virus Evol. 2015, 1, vev014. [Google Scholar] [CrossRef] [Green Version]
  40. Walter, C.T.; Barr, J.N. Recent advances in the molecular and cellular biology of bunyaviruses. J. Gen. Virol. 2011, 92, 2467–2484. [Google Scholar] [CrossRef]
  41. Abudurexiti, A.; Adkins, S.; Alioto, D.; Alkhovsky, S.V.; Avšič-Županc, T.; Ballinger, M.J.; Bente, D.A.; Beer, M.; Bergeron, É.; Blair, C.D.; et al. Taxonomy of the order Bunyavirales: Update 2019. Arch. Virol. 2019, 164, 1949–1965. [Google Scholar] [CrossRef] [Green Version]
  42. Osaki, H.; Sasaki, A.; Nomiyama, K.; Tomioka, K. Multiple virus infection in a single strain of Fusarium poae shown by deep sequencing. Virus Genes 2016, 52, 835–847. [Google Scholar] [CrossRef]
  43. Walker, T.; Jeffries, C.L.; Mansfield, K.L.; Johnson, N. Mosquito cell lines: History, isolation, availability and application to assess the threat of arboviral transmission in the United Kingdom. Parasit. Vectors 2014, 7, 382. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Brackney, D.E.; Scott, J.C.; Sagawa, F.; Woodward, J.E.; Miller, N.A.; Schilkey, F.D.; Mudge, J.; Wilusz, J.; Olson, K.E.; Blair, C.D.; et al. C6/36 Aedes albopictus cells have a dysfunctional antiviral RNA interference response. PLoS Negl. Trop. Dis. 2010, 4, e856. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Infectious Disease Surveillance Center, National Institute of Infectious Diseases (NIID) Japan: Japanese encephalitis, Japan, 2003–2008. Infect. Agents Surveill. Rep. (IASR) 2009, 30, 147–148. Available online: https://idsc.niid.go.jp/iasr/30/352/tpc352.html (accessed on 14 January 2020).
  46. Infectious Disease Surveillance Center, NIID Japan: Japanese encephalitis, Japan, 2007–2016. Infectious Infect. Agents Surveill. Rep. (IASR) 2017, 38, 151–152. Available online: https://www.niid.go.jp/niid/en/iasr-vol38-e/865-iasr/7483-450te.html (accessed on 14 January 2020).
  47. Tadano, M.; Kanemura, K.; Hasegawa, H.; Makino, Y.; Fukunaga, T. Epidemiological and ecological studies of Japanese encephalitis in Okinawa, subtropical area in Japan. I. Investigations on antibody levels to Japanese encephalitis virus in swine sera and vector mosquito in Okinawa, Miyako and Ishigaki islands. Microbiol. Immunol. 1994, 38, 117–122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Cholleti, H.; Hayer, J.; Abilio, A.P.; Mulandane, F.C.; Verner-Carlsson, J.; Falk, K.I.; Fafetine, J.M.; Berg, M.; Blomström, A.L. Discovery of novel viruses in mosquitoes from the Zambezi valley of Mozambique. PLoS ONE 2016, 11, e0162751. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Shi, C.; Beller, L.; Deboutte, W.; Yinda, K.C.; Delang, L.; Vega-Rúa, A.; Failloux, A.B.; Matthijnssens, J. Stable distinct core eukaryotic viromes in different mosquito species from Guadeloupe, using single mosquito viral metagenomics. Microbiome 2019, 7, 121. [Google Scholar] [CrossRef]
  50. Ritchie, S.A.; Rochester, W. Wind-blown mosquitoes and introduction of Japanese encephalitis into Australia. Emerg. Infect. Dis. 2001, 7, 900–903. [Google Scholar] [CrossRef]
  51. Inglis, T.J.; O’Reilly, L.; Merritt, A.J.; Levy, A.; Heath, C.H. The aftermath of the Western Australian melioidosis outbreak. Am. J. Trop. Med. Hyg. 2011, 84, 851–857. [Google Scholar] [CrossRef] [Green Version]
  52. Johansen, C.A.; van den Hurk, A.F.; Pyke, A.T.; Zborowski, P.; Phillips, D.A.; Mackenzie, J.S.; Ritchie, S.A. Entomological investigations of an outbreak of Japanese encephalitis virus in the Torres Strait, Australia, in 1998. J. Med. Entomol. 2001, 38, 581–588. [Google Scholar] [CrossRef] [Green Version]
  53. Kuwata, R.; Isawa, H.; Hoshino, K.; Sasaki, T.; Kobayashi, M.; Maeda, K.; Sawabe, K. Analysis of mosquito-borne flavivirus superinfection in Culex tritaeniorhynchus (Diptera: Culicidae) cells persistently infected with Culex Flavivirus (Flaviviridae). J. Med. Entomol. 2015, 52, 222–229. [Google Scholar] [CrossRef] [PubMed]
  54. National Epidemiological Surveillance of Vaccine-Preventable Diseases: Seroprevalence of Japanese Encephalitis in Swine during Summer Season. Available online: https://www.niid.go.jp/niid/ja/y-graphs/1600-yosoku-index-e.html (accessed on 7 January 2020).
  55. Scherer, W.F.; Buescher, E.L. Ecologic studies of Japanese encephalitis virus in Japan: I. Introduction. Am. J. Trop. Med. Hyg. 1959, 8, 644–650. [Google Scholar] [CrossRef] [PubMed]
  56. Kobayashi, D.; Isawa, H.; Ejiri, H.; Sasaki, T.; Sunahara, T.; Futami, K.; Tsuda, Y.; Katayama, Y.; Mizutani, T.; Minakawa, N.; et al. Complete genome sequencing and phylogenetic analysis of a Getah virus strain (Genus Alphavirus, Family Togaviridae) isolated from Culex tritaeniorhynchus mosquitoes in Nagasaki, Japan in 2012. Vector Borne Zoonotic Dis. 2016, 16, 769–776. [Google Scholar] [CrossRef] [PubMed]
  57. Bennett, A.J.; Bushmaker, T.; Cameron, K.; Ondzie, A.; Niama, F.R.; Parra, H.J.; Mombouli, J.V.; Olson, S.H.; Munster, V.J.; Goldberg, T.L. Diverse RNA viruses of arthropod origin in the blood of fruit bats suggest a link between bat and arthropod viromes. Virology 2019, 528, 64–72. [Google Scholar] [CrossRef] [PubMed]
  58. Li, L.; Victoria, J.G.; Wang, C.; Jones, M.; Fellers, G.M.; Kunz, T.H.; Delwart, E. Bat guano virome: Predominance of dietary viruses from insects and plants plus novel mammalian viruses. J. Virol. 2010, 84, 6955–6965. [Google Scholar] [CrossRef] [Green Version]
  59. Kading, R.C.; Schountz, T. Flavivirus infections of bats: Potential role in Zika virus ecology. Am. J. Trop. Med. Hyg. 2016, 95, 993–996. [Google Scholar] [CrossRef] [Green Version]
  60. Saepulloh, M.; Dharmayanti, N.L.P.I.; Adjid, R.M.A.; Ratnawati, A.; Sendow, I. The presence of Japanese encephalitis virus infection in Pteropus sp. in West Kalimantan. Proc. Int. Sem. Livest. Prod. Vet. Technol. 2016, 549–553. [Google Scholar] [CrossRef]
  61. Wang, J.L.; Pan, X.L.; Zhang, H.L.; Fu, S.H.; Wang, H.Y.; Tang, Q.; Wang, L.F.; Liang, G.D. Japanese encephalitis viruses from bats in Yunnan, China. Emerg. Infect. Dis. 2009, 15, 939–942. [Google Scholar] [CrossRef]
  62. Hall-Mendelin, S.; Allcock, R.; Kresoje, N.; van den Hurk, A.F.; Warrilow, D. Detection of arboviruses and other micro-organisms in experimentally infected mosquitoes using massively parallel sequencing. PLoS ONE 2013, 8, e58026. [Google Scholar] [CrossRef] [Green Version]
  63. White, L.A. Susceptibility of Aedes albopictus C6/36 cells to viral infection. J. Clin. Microbiol. 1987, 25, 1221–1224. [Google Scholar] [CrossRef] [Green Version]
Viruses 12 00264 g001 550
Figure 1. Map showing the mosquito collection sites in three prefectures of Japan: Ishikawa, Tottori, and Nagasaki.
Figure 1. Map showing the mosquito collection sites in three prefectures of Japan: Ishikawa, Tottori, and Nagasaki.
Viruses 12 00264 g001
Viruses 12 00264 g002 550
Figure 2. Classification of mosquitoes used in the present study. Three JE major vectors (C. tritaeniorhynchus, C. pseudovishnui, and C. vishnui) belong to the C. vishnui subgroup. Note that the C. vishnui subgroup and the outlier (C. inatomii) differ in the subgenus level.
Figure 2. Classification of mosquitoes used in the present study. Three JE major vectors (C. tritaeniorhynchus, C. pseudovishnui, and C. vishnui) belong to the C. vishnui subgroup. Note that the C. vishnui subgroup and the outlier (C. inatomii) differ in the subgenus level.
Viruses 12 00264 g002
Viruses 12 00264 g003 550
Figure 3. (A) Viral read ratio of C. tritaeniorhynchus identified virus families based on collection sites (Ishikawa, Tottori, and Nagasaki prefectures) along with the natural hosts of each virus family for better illustration. In general, C. tritaeniorhynchus collected from different locations showed the same virus families with variations in NGS reads. (B) The viral read ratio of all identified virus families based on mosquito species. The three mosquito species (C. tritaeniorhynchus, C. pseudovishnui, and C. inatomii) displayed different sets of virus families. Accordingly, the virome diversity of C. tritaeniorhynchus and C. pseudovishnui is richer than C. inatomii, which belongs to a different subgenus.
Figure 3. (A) Viral read ratio of C. tritaeniorhynchus identified virus families based on collection sites (Ishikawa, Tottori, and Nagasaki prefectures) along with the natural hosts of each virus family for better illustration. In general, C. tritaeniorhynchus collected from different locations showed the same virus families with variations in NGS reads. (B) The viral read ratio of all identified virus families based on mosquito species. The three mosquito species (C. tritaeniorhynchus, C. pseudovishnui, and C. inatomii) displayed different sets of virus families. Accordingly, the virome diversity of C. tritaeniorhynchus and C. pseudovishnui is richer than C. inatomii, which belongs to a different subgenus.
Viruses 12 00264 g003
Viruses 12 00264 g004 550
Figure 4. Venn diagram showing viruses identified from C. tritaeniorhynchus collected from three prefectures. Note the eight viruses in the overlapping region.
Figure 4. Venn diagram showing viruses identified from C. tritaeniorhynchus collected from three prefectures. Note the eight viruses in the overlapping region.
Viruses 12 00264 g004
Viruses 12 00264 g005 550
Figure 5. (A) Genome organization of two novel totiviruses found in this study. Note the difference in the overlapping region. (B) Maximum likelihood phylogeny of Totiviridae and Chrysoviridae clade using ~123 amino acid sequences of RdRp conserved domains selected by Gblocks. Values on branches indicate bootstrap support based on 1000 bootstrap replicates. Bootstrap values <70 are not shown. ► indicates novel viruses and ➢ indicates a new strain identified in this study.
Figure 5. (A) Genome organization of two novel totiviruses found in this study. Note the difference in the overlapping region. (B) Maximum likelihood phylogeny of Totiviridae and Chrysoviridae clade using ~123 amino acid sequences of RdRp conserved domains selected by Gblocks. Values on branches indicate bootstrap support based on 1000 bootstrap replicates. Bootstrap values <70 are not shown. ► indicates novel viruses and ➢ indicates a new strain identified in this study.
Viruses 12 00264 g005
Viruses 12 00264 g006 550
Figure 6. (A) Genome organization of two novel partitiviruses found in this study. (B) Maximum likelihood phylogeny of the Partitiviridae using ~163 amino acid sequences of RdRp conserved domains selected by Gblocks. Values on branches indicate bootstrap support based on 1000 bootstrap replicates. Bootstrap values <70 are not shown. Viruses identified in this study are in bold, with ► indicating novel viruses and ➢ indicating a new strain identified in this study.
Figure 6. (A) Genome organization of two novel partitiviruses found in this study. (B) Maximum likelihood phylogeny of the Partitiviridae using ~163 amino acid sequences of RdRp conserved domains selected by Gblocks. Values on branches indicate bootstrap support based on 1000 bootstrap replicates. Bootstrap values <70 are not shown. Viruses identified in this study are in bold, with ► indicating novel viruses and ➢ indicating a new strain identified in this study.
Viruses 12 00264 g006
Viruses 12 00264 g007 550
Figure 7. (A) Genome organization of a novel flavi-like virus found in this study. (B) Maximum likelihood phylogeny of Flaviviridae using ~297 amino acid sequences of the NS5 region. Values on branches indicate bootstrap support based on 1000 bootstrap replicates. Bootstrap values <70 are not shown. Viruses identified in this study are in bold, with ► indicating novel viruses and ➢ indicating a new strain identified in this study.
Figure 7. (A) Genome organization of a novel flavi-like virus found in this study. (B) Maximum likelihood phylogeny of Flaviviridae using ~297 amino acid sequences of the NS5 region. Values on branches indicate bootstrap support based on 1000 bootstrap replicates. Bootstrap values <70 are not shown. Viruses identified in this study are in bold, with ► indicating novel viruses and ➢ indicating a new strain identified in this study.
Viruses 12 00264 g007
Viruses 12 00264 g008 550
Figure 8. Maximum likelihood phylogeny of JEV using ~10,995 nucleotide sequences. Values on branches indicate bootstrap support based on 1000 bootstrap replicates. Bootstrap values <70 are not shown. ➢ indicates a new strain identified in this study.
Figure 8. Maximum likelihood phylogeny of JEV using ~10,995 nucleotide sequences. Values on branches indicate bootstrap support based on 1000 bootstrap replicates. Bootstrap values <70 are not shown. ➢ indicates a new strain identified in this study.
Viruses 12 00264 g008
Viruses 12 00264 g009 550
Figure 9. (A) Genome organization of four novel (+)ssRNA viruses found in this study. (B) Phylogenetic analysis of several (+)ssRNA virus families was performed with the maximum likelihood method using ~4396 amino acid sequences of the RdRp region. Node bootstraps were calculated with 1000 replicates. Bootstrap values <70 are not shown. Viruses identified in this study are in bold, with ► indicating novel viruses and ➢ indicating a new strain identified in this study.
Figure 9. (A) Genome organization of four novel (+)ssRNA viruses found in this study. (B) Phylogenetic analysis of several (+)ssRNA virus families was performed with the maximum likelihood method using ~4396 amino acid sequences of the RdRp region. Node bootstraps were calculated with 1000 replicates. Bootstrap values <70 are not shown. Viruses identified in this study are in bold, with ► indicating novel viruses and ➢ indicating a new strain identified in this study.
Viruses 12 00264 g009
Viruses 12 00264 g010 550
Figure 10. Genome organization of two viruses found in this study. The appearance of two viruses in the same mosquito pools has inferred the possibility of the two sequences being the same virus. Motts Mill virus (accession number MH384280) is used as reference genome.
Figure 10. Genome organization of two viruses found in this study. The appearance of two viruses in the same mosquito pools has inferred the possibility of the two sequences being the same virus. Motts Mill virus (accession number MH384280) is used as reference genome.
Viruses 12 00264 g010
Viruses 12 00264 g011 550
Figure 11. (A) Genome organization of two novel iflaviruses found in this study. (B) Phylogenetic analysis of Iflaviridae was performed with the maximum likelihood method using ~971 amino acid sequences of RdRp conserved domains selected by Gblocks. Node bootstraps were calculated with 1000 replicates. Bootstrap values <70 are not shown. Viruses identified in this study are in bold, with ► indicating novel viruses and ➢ indicating a new strain identified in this study.
Figure 11. (A) Genome organization of two novel iflaviruses found in this study. (B) Phylogenetic analysis of Iflaviridae was performed with the maximum likelihood method using ~971 amino acid sequences of RdRp conserved domains selected by Gblocks. Node bootstraps were calculated with 1000 replicates. Bootstrap values <70 are not shown. Viruses identified in this study are in bold, with ► indicating novel viruses and ➢ indicating a new strain identified in this study.
Viruses 12 00264 g011
Viruses 12 00264 g012 550
Figure 12. (A) Genome organization of four novel (−)ssRNA viruses found in this study. (B) Phylogenetic analysis of several (−)ssRNA virus families was performed with the maximum likelihood method using ~4556 amino acid sequences of the RdRp region. Node bootstraps were calculated with 1000 replicates. Bootstrap values <70 are not shown. Viruses identified in this study are in bold, with ► indicating novel viruses and ➢ indicating a new strain identified in this study.
Figure 12. (A) Genome organization of four novel (−)ssRNA viruses found in this study. (B) Phylogenetic analysis of several (−)ssRNA virus families was performed with the maximum likelihood method using ~4556 amino acid sequences of the RdRp region. Node bootstraps were calculated with 1000 replicates. Bootstrap values <70 are not shown. Viruses identified in this study are in bold, with ► indicating novel viruses and ➢ indicating a new strain identified in this study.
Viruses 12 00264 g012
Table
Table 1. The profile of mosquito samples employed in metagenomic analysis and virus isolation.
Table 1. The profile of mosquito samples employed in metagenomic analysis and virus isolation.
Mosquito SpeciesMosquito Collection SiteGPS CoordinatesCollection MonthMethods
Metagenomic AnalysisVirus Isolation
No. of Pools TestedNo. of Female Mosquitoes TestedNo. of Pools TestedNo. of Female Mosquitoes TestedJEV-Positive Pools
C. tritaeniorhynchusMonzen-machi, Wajima city, Ishikawa37°29′ N, 136°74′ EJun–Oct5250123001
2017
C. tritaeniorhynchusYonago waterbirds sanctuary, Yonago city, Tottori35°26′ N, 133°17′ EJun–Sep315082000
2017
C. tritaeniorhynchusTeramine farm, Isahaya city, Nagasaki32°49′ N, 130°03′ EAug2100541,3500
2017
C. pseudovishnuiTeramine farm, Isahaya city, Nagasaki32°49′ N, 130°03′ EAug2100174210
2017
C. inatomiiYonago waterbirds sanctuary, Yonago city, Totttori35°26′ N, 133°17′ EJun1522520
2017
Total 136529323231
Table
Table 2. Virome identification result (representative genomes or RdRp segments only).
Table 2. Virome identification result (representative genomes or RdRp segments only).
Virus
Category		VIRUS TAXONVirus NameClosely Related Viruses
(Accession no.)		BlastnBlastxLocation †/Species Origin ‡SegmentAccession no.
[This Study]		Complete CDs
Identity
(%) *		Query
Cover (%)		e-ValueIdentity
(%) *		Query
Cover (%)		e-Value
dsRNATotiviridaeCulex vishnui subgroup totivirus (CvsTV)dsRNA virus environmental sample (AJT39583)---46390.0Japan/CtrCps-LC514295Y
TotiviridaeCulex inatomii totivirus (CiTV)Australian Anopheles Totivirus (ASU43981)---61370.0Tottori/Cnt-LC514398Y
PartitiviridaeCulex pseudovishnui partitivirus (CpPV)Hubei partiti-like virus 56 (APG78242)---56850.0Nagasaki/Cps1LC514399Y
PartitiviridaeHubei partiti-like virus 22 (HPLV22)(APG78283)99980.099980.0Tottori/Ctr
Nagasaki/Ctr		1LC514400
LC514401		Y
Y
PartitiviridaeCulex tritaeniorhynchus partitivirus (CtPV)Hubei partiti-like virus 19 (APG78260)---49841 × 10−170Japan/Ctr1LC514402Y
ChrysoviridaeHubei chryso-like virus 1 (HCLV1)(ASA47395)84990.0---Ishikawa/Ctr1LC514396Y
(+)ssRNAFlaviviridaeCulex tritaeniorhynchus flavi-like virus (CtFLV)Shayang fly virus 4 (YP_009179225)---41710.0Japan/Ctr-LC514290N
FlaviviridaeMosquito flavivirus (MFV)(BAR88121)98990.099930.0Nagasaki/Ctr
Nagasaki/Cps		-LC513840
LC513841		Y
Y
FlaviviridaeJapanese Encephalitis Virus (JEV)(AB981184)991000.099970.0Ishikawa/Ctr-LC513838Y
FlaviviridaeCulex Flavivirus (CxFV)(BAM74417)971000.099930.0Nagasaki/Ctr-LC513839Y
Negevirus-relatedCulex tritaeniorhynchus negev-like virus (CtNLV)Mill Lade virus (QAY29259)---47220.0Japan/Ctr-LC507097Y
Negevirus-relatedCulex pseudovishnui negev-like virus (CpNLV)Yongsan negev-like virus 1 (AXV43886)---87730.0Nagasaki/Cps-LC512731Y
(+)ssRNATymoviralesCulex pseudovishnui tymo-like virus (CpTLV)Tarnsjo virus (QGA70928)79760.074840.0Nagasaki/Cps-LC512732Y
Luteoviridae-relatedHubei mosquito virus 2 (HMV2)APG7562899990.0100340.0Ishikawa/Ctr
Tottori/Ctr		1
1		LC513829
LC513830		Y
Y
Luteoviridae-relatedCulex inatomii luteo-like virus (CiLLV)Hubei mosquito virus 2 (APG75628)83990.0---Tottori/Cnt1LC513833Y
Tombusviridae-relatedHubei mosquito virus 4 (HMV4)(APG76308)91820.0---Tottori/Ctr
Nagasaki/Cnt		-LC512733
LC512734		Y
N
90980.0---Tottori/Ctr
Nagasaki/Cnt		-LC512735
LC512736		N
N
Tombusviridae-relatedWenzhou tombus-like virus 11 (WTLV11)(YP_009342051)99990.0100280.0Nagasaki/Ctr
Nagasaki/Cps		-LC512737
LC512738		Y
Y
IflaviridaeIsahaya Culex Iflavirus (ICIFV)Wuhan fly virus 4 (YP 009342337)---43710.0Nagasaki/Cps-LC513835Y
IflaviridaeYongsan Iflavirus 1 (YIFV1)(AXV43887)961000.097960.0Tottori/Ctr-LC513837N
IflaviridaeYonago Culex Iflavirus (YCIFV)Wuhan insect virus 13 (YP 009342321)---38820.0Tottori/Ctr-LC513836Y
Sobemovirus-relatedWenzhou sobemo-like virus 3 (WSLV3)(APG75759)91920.096470.0Ishikawa/Ctr
Tottori/Ctr
Nagasaki/Ctr
Nagasaki/Cps		-LC512854
LC512855
LC512856
LC512857		N
N
N
Y
Sobemovirus-relatedBat sobemovirus (BSV)(AGN73380)90420.094424×10−133Ishikawa/Ctr
Tottori/Ctr
Nagasaki/Ctr
Nagasaki/Cps		-LC512858
LC512859
LC512860
LC512861		N
N
Y
Y
(−)ssRNAXinmoviridaeCulex tritaeniorhynchus anphevirus (CtAV)Guadeloupe mosquito mononega-like virus (QEM39171.1)---38470.0Ishikawa/Ctr
Tottori/Ctr		-LC514054
LC514055		Y
Y
RhabdoviridaeCulex tritaeniorhynchus rhabdovirus (CTRV)(BAU46576)991000.099560.0Japan/Ctr-LC514403Y
RhabdoviridaeCulex pseudovishnui rhabdo-like virus (CpRLV)Tongilchon virus 1 (YP_009182186)76550.088540.0Nagasaki/Cps
Nagasaki/Cps		-LC514056
LC514057		Y
Y
BunyaviralesCulex pseudovishnui bunya-like virus (CpBLV)Bunyaviridae environmental sample (AJT39594)80980.0---Nagasaki/Cps
Nagasaki/Cps		LLC514291
LC514293		N
Y
Aspiviridae-relatedCulex tritaeniorhynchus Aspi-like virus (CtALV)Fusarium poae negative-stranded virus (YP_009272911)---30800.0Tottori/Ctr-LC514058Y
Novel viruses are highlighted in blue. * Only the RNA-dependent RNA polymerase (RdRp) region’s identity is shown, except for bat sobemovirus. † Location is one of the following: Ishikawa, Tottori, Nagasaki, or Japan (when sequences obtained from different locations were identical). ‡ Species origin is abbreviations of Ctr (C. tritaeniorhynchus), Cps (C. pseudovishnui), and Cnt (C. inatomii). complete CDs: Yes (Y) and No (N).
Table
Table 3. Viral diversity based on mosquito species and collection sites.
Table 3. Viral diversity based on mosquito species and collection sites.
Virus NameMosquito Groups *
Ctr-IshikawaCtr-TottoriCtr-NagasakiCps-NagasakiCnt-Tottori
WSLV3++++−
BSV++++−
CvsTV++++−
HPLV22+++−−
CtPV+++−−
CtNLV+++−−
CtAV+++−−
CtFLV+++−−
HCLV1++−−−
CTRV++−−−
HMV2++−−−
HMV4−+−−+
ICIFV−−++−
WTLV11−−++−
MFV−−++−
CxFV−−+−−
JEV+−−−−
YCIFV−+−−−
CtALV−+−−−
YIFV−+−−−
CpPV−−−+−
CpNLV−−−+−
CpBLV−−−+−
CpRLV−−−+−
CpTLV−−−+−
CiTV−−−−+
CiLLV−−−−+
* mosquito group abbreviations are Ctr (C. tritaeniorhynchus), Cps (C. pseudovishnui), and Cnt (C. inatomii).

Share and Cite

MDPI and ACS Style

Faizah, A.N.; Kobayashi, D.; Isawa, H.; Amoa-Bosompem, M.; Murota, K.; Higa, Y.; Futami, K.; Shimada, S.; Kim, K.S.; Itokawa, K.; et al. Deciphering the Virome of Culex vishnui Subgroup Mosquitoes, the Major Vectors of Japanese Encephalitis, in Japan. Viruses 2020, 12, 264. https://0-doi-org.brum.beds.ac.uk/10.3390/v12030264

AMA Style

Faizah AN, Kobayashi D, Isawa H, Amoa-Bosompem M, Murota K, Higa Y, Futami K, Shimada S, Kim KS, Itokawa K, et al. Deciphering the Virome of Culex vishnui Subgroup Mosquitoes, the Major Vectors of Japanese Encephalitis, in Japan. Viruses. 2020; 12(3):264. https://0-doi-org.brum.beds.ac.uk/10.3390/v12030264

Chicago/Turabian Style

Faizah, Astri Nur, Daisuke Kobayashi, Haruhiko Isawa, Michael Amoa-Bosompem, Katsunori Murota, Yukiko Higa, Kyoko Futami, Satoshi Shimada, Kyeong Soon Kim, Kentaro Itokawa, and et al. 2020. "Deciphering the Virome of Culex vishnui Subgroup Mosquitoes, the Major Vectors of Japanese Encephalitis, in Japan" Viruses 12, no. 3: 264. https://0-doi-org.brum.beds.ac.uk/10.3390/v12030264

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop