Next Article in Journal
Biomarkers of Progression after HIV Acute/Early Infection: Nothing Compares to CD4+ T-cell Count?
Previous Article in Journal
Causation of Acute Flaccid Paralysis by Myelitis and Myositis in Enterovirus-D68 Infected Mice Deficient in Interferon αβ/γ Receptor Deficient Mice
Previous Article in Special Issue
Myelin Oligodendrocyte Glycoprotein-Independent Rubella Infection of Keratinocytes and Resistance of First-Trimester Trophoblast Cells to Rubella Virus In Vitro
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metagenomic Virome Analysis of Culex Mosquitoes from Kenya and China

1
Key Laboratory of Agricultural and Environmental Microbiology, Wuhan Institute of Virology, Chinese Academy of Sciences, Wuhan 430071, China
2
Key Laboratory of Special Pathogens, Wuhan Institute of Virology, Chinese Academy of Sciences, Wuhan 430071, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
4
Veterinary Services Department, Kenya Wildlife Service, P.O. Box 40241, Nairobi 00100, Kenya
5
Mammalogy Section, National Museum of Kenya, P.O. Box 40658, Nairobi 00100, Kenya
6
Entomology Section, National Museum of Kenya, P.O. Box 40658, Nairobi 00100, Kenya
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2018, 10(1), 30; https://doi.org/10.3390/v10010030
Submission received: 4 December 2017 / Revised: 3 January 2018 / Accepted: 9 January 2018 / Published: 12 January 2018
(This article belongs to the Collection Global Viral Threats)

Abstract

:
Many blood-feeding arthropods are known vectors of viruses that are a source of unprecedented global health concern. Mosquitoes are an integral part of these arthropod vectors. Advancements in next-generation sequencing and bioinformatics has expanded our knowledge on the richness of viruses harbored by arthropods. In the present study, we applied a metagenomic approach to determine the intercontinental virome diversity of Culex quinquefasciatus and Culex tritaeniorhynchus in Kwale, Kenya and provinces of Hubei and Yunnan in China. Our results showed that viromes from the three locations were strikingly diverse and comprised 30 virus families specific to vertebrates, invertebrates, plants, and protozoa as well as unclassified group of viruses. Though sampled at different times, both Kwale and Hubei mosquito viromes were dominated by vertebrate viruses, in contrast to the Yunnan mosquito virome, which was dominated by insect-specific viruses. However, each virome was unique in terms of virus proportions partly influenced by type of ingested meals (blood, nectar, plant sap, environment substrates). The dominant vertebrate virus family in the Kwale virome was Papillomaviridae (57%) while in Hubei it was Herpesviridae (30%) and the Yunnan virome was dominated by an unclassified viruses group (27%). Given that insect-specific viruses occur naturally in their hosts, they should be the basis for defining the viromes. Hence, the dominant insect-specific viruses in Kwale, Hubei, and Yunnan were Baculoviridae, Nimaviridae and Iflaviridae, respectively. Our study is preliminary but contributes to growing and much needed knowledge, as mosquito viromes could be manipulated to prevent and control pathogenic arboviruses.

Graphical Abstract

1. Introduction

Global health security is increasingly threatened by emerging and re-emerging pathogens whose consequences to human and animal health are dire and devastating to economies. Most of these pathogens have their origin in diverse hosts, mainly wildlife, and are vectored by multiple blood-feeding arthropods. Part of the global health challenge is that some of these arthropods are widely distributed and highly efficient vectors for many pathogens. Specifically, the mosquito is one blood-feeding arthropod, known to feed on myriad hosts, and is an efficient vector for a multitude of virulent pathogens such as Plasmodium spp., filarial nematodes, Zika virus, dengue virus, Rift Valley fever virus and yellow fever virus. The diversity and richness of mosquito flora is only beginning to be understood following the application of next-generation sequencing (NGS) and vector-enabled metagenomics (VEM). These advanced molecular platforms have revealed that in addition to known pathogens, mosquitoes harbor a large proportion of organisms whose pathogenicity in vertebrate hosts is not yet known. Previous metagenomic studies have revealed that mosquitoes harbor a highly diverse microbiome composed of bacteria, protozoa, and viruses, of which some are specific to arthropods, protozoans, and plants [1,2,3,4]. The microbiome community of the mosquito is shaped by three broad factors; the vector itself, in vivo microorganisms in the mosquito body, and lastly the mosquito meal sources.
There is a developing specific interest to understand mosquito viromes, particularly in terms of virus diversity and composition. This arises from the fact that most mosquito-borne viruses bear significant threat to public and animal health, whereas a large proportion (<70% amino acid identity) of recovered viruses are of unknown pathogenicity [5,6]. Surveillance of animal and human viruses to determine the composition and proportion of circulating viruses in specific mosquitoes within a region is a necessary challenge because the outcome is fundamental and critical for formulating strategies on prediction, preparedness, prevention, and eventual control of any disease. Thus, mosquito virome studies have been proposed as a promising robust tool that could be applied in virus surveillance [1]. This is based on the fact that both DNA and RNA viruses can be detected from mosquito blood meal for up to 24 h [7], coupled by the reality that the diversity includes viruses not biologically vectored by the specific mosquito.
Besides viruses that can replicate in vertebrates, mosquito viromes seem to harbor a large proportion of viruses that naturally infect insects and can replicate only in insect cell lines in vitro. Even though these insect-specific viruses seem not to replicate in vertebrates [8], some of them are phylogenetically related to known pathogenic arboviruses and it is speculated that they may evolve to become novel pathogenic viruses in vertebrates [9]. Further, due to their close evolutionary association with known arboviruses, they are being hypothesized for potential biocontrol applications (modulation of transmission of arboviruses) or as a basis for new vaccine formulations [9,10,11,12]. Nevertheless, due to their abundance and diversity (including the virus families Togaviridae, Flaviviridae, Rhabdoviridae, Bunyaviridae and Mesoniviridae), they constitute a significant component of mosquito viromes. In addition, it seems that each virome is dominated by particular viruses, which hypothetically may be transient, but significant in distinguishing or describing diverse viromes. For example, the virome of Culex spp. in Zambezi, Mozambique was dominated (88.5%) by insect-specific viruses, of the Iflaviridae family [13]. In China, a previous study of Culex tritaeniorhynchus virome in Hubei province also showed a similarly high dominance (88%) of insect-specific viruses [4].
Increased frequency and speed of human travel and commercial trade has been implicated in the expansion of mosquito vectors and their pathogens from indigenous habitats to new areas. The well-known Aedes albopictus first made its incursion to South Africa through trade in old car tires [14] and has since been established in Africa. We hypothesized that because of the increased interaction between China and Kenya, there is high likelihood for exchange of both mosquito and viruses. For instance, Ochieng et al. (2013) isolated a flavivirus from Rabai town, coastal Kenya that had 77% nucleotide similarity to Chaoyang virus from Shenyang, China. Therefore, our objective was to determine viromes of Culex mosquitoes in specific sites in China and Kenya. The Culex mosquito is an efficient vector for several arboviruses such as West Nile virus, Rift Valley fever virus, and Japanese encephalitis, as well as parasitic bancroftian filariasis. African mosquito virome studies are in their nascent stage, with few exploratory investigations, though outcomes from such studies show remarkable diversity and composition [2,13]. Our study will provide important preliminary information on virus status in populations of Culex spp. in both Kenya and China.

2. Materials and Methods

2.1. Study Sampling Sites and Related Important Ecological Features

The study areas in Kenya were located in the Kwale and Kilifi counties, along the Indian Ocean coastline. The study areas in China were located in the Hubei province in south-central china and the Yunnan provinces in south-western China.

2.1.1. Kwale County—Kenya

Kwale county is located on the south-eastern part of Kenya along the shores of the Indian ocean (Figure 1). It lies at 5 m (18 ft) above sea level. It has a tropical humid and temperate climate all year. The annual mean temperature is 26.3 °C (79.3 °F). In the months of January and February, the average temperature is 31 °C (88 °F) and an average precipitation of 22 mm (Source: Kenya Meteorological Department). The forest cover values for the three counties flagging the Indian ocean coastline are 5.44%, 7.67% and 33.9% for Kwale county, Kilifi country and Lamu county, respectively (Source: http://www.opendata.go.ke). The region has experienced several arbovirus outbreaks due to its high human population together with high inflow and outflow of visiting tourists. Also, the close vector/human/wildlife proximity might be a contributing factor to the emerging outbreaks in the region.

2.1.2. Kilifi County—Kenya

Kilifi county lies along the coastal belt, in the south-eastern part of Kenya. The region receives an average annual rainfall of 900–1100 mm. The annual temperatures ranges from 21 °C to 30 °C in the coastal belt but rise to maximum of 34 °C in the hinterland (Source: http://www.kilifi.go.ke). There is a precipitation difference of 246 mm between the driest and wettest months. The annual temperature variation is 3.8 °C. Due to its proximity to the Indian Ocean, tourism and fishing are the main economic activities hence experiencing a high human traffic.

2.1.3. Hubei Province—China

Hubei Province is located in south-central China (Figure 1). It has extreme weather conditions, (more so during the summer) and is humid all year round due to high number of water bodies in the region. It features a subtropical monsoon climate with an average annual temperature of about 15 °C (59 °F). The province has four distinct seasons with a burning hot summer (June, July and August) and a chilly winter (December, January and February). The hottest month, July, averages 27–30 °C and the coldest month, January, 1–5 °C. The average annual precipitation is between 800–1600 mm. In Hubei, summer is the wettest season with 300–700 mm in average precipitation, while winter is the driest season with average rainfall at 30–190 mm. The average temperature for the months of July to September is 30 °C (83 °F) (Source: http://en.hubei.gov.cn).

2.1.4. Yunnan Province—China

Yunnan Province is located in the south-western part of China (Figure 1) and it has a characteristic feature of “spring-like” weather with very good forest cover. The region sits at an altitude of 1000 m above sea level with abundant heat and light, a spring-like climate. The annual average temperature is 17 °C. The average temperature for the months of July to Sept is 22 °C and average precipitation is 21 mm. Yunnan province in China serves as the main business trading route between countries in Southeast Asia and China and it shares a lengthy border with three countries: Myanmar, Laos and Vietnam. This region is seen as the main epicenter of emerging or re-emerging viral diseases due to its multi-ethnic society and poorly developed public health structures, as well as abundant and diverse wild animal resources and their illegal trading [15,16,17,18].

2.2. Sample Collection and Taxonomic Identification

Mosquito trappings for adult stages were set at randomly selected sites in Kwale and Kilifi county, on the Kenyan coast. The mosquitoes were trapped alive using carbon dioxide mosquito traps (Mosquito Magnet® Woodstream Corp., Brampton, ON, Canada) from outdoor sites between dusk and dawn, in the months of January to February 2017. In China, adult mosquitoes were sampled from two provinces, namely the Hubei province (Coordinates: 30.8845° N, 112.5923° E) and Yunnan province (Coordinates: 24.97564° N, 101.4848° E) (Figure 1). In Hubei, sampling was done between the month of July and September of the year 2014 while in the Yunnan province it was done from July to September of the year 2015. In these two sites in china, electrical mosquito aspirators (Yalin, Zhejiang, China) were used in mosquito capture from homesteads and surrounding regions. In all the study sites, the mosquitoes were killed by dipping in liquid nitrogen. The killed mosquitoes were morphologically identified as described by Edwards (1941) and Hopkins (1952) [19,20]. The mosquitoes were grouped into pools of 50 and stored at −80 °C until further analysis.

2.3. Viral RNA Isolation, cDNA Library Preparation, Amplification and Sequencing

The mosquito pools were homogenized in 1 ml Roswell Park Memorial Institute medium (RPMI) and the homogenate filtered through a 0.22 μM polycarbonate filter membrane (Millipore, Billerica, MA, USA) to exclude bacteria and other debris. RNA was extracted from 140 μL of homogenate using the QIAamp viral RNA extraction kit according to manufacturer’s instructions (Qiagen, GmBH, Hilden, Germany). The RNA integrity was checked using the Agilent bioanalyzer 2100 (Agilent Tech, Santa Clara, CA, USA) and further purified using the RNA clean XP kit (Beckham Coulter Inc., Brea, CA, USA). Strand-specific libraries were prepared using the TruSeq® Stranded Total RNA Sample Preparation kit (Illumina, San Diego, CA, USA) following the manufacturer’s instructions. Briefly, ribosomal RNA was removed from total RNA using Ribo-Zero rRNA removal beads. Following purification, the mRNA was fragmented into small pieces using divalent cations under 94 °C for 8 min. The cleaved RNA fragments were copied into first strand cDNA using reverse transcriptase and random primers. This was followed by second strand cDNA synthesis using DNA Polymerase I and RNase H. These cDNA fragments then went through an end repair process, the addition of a single “A” base, and then ligation of the adapters. The products were then purified and enriched with PCR to create the final cDNA library. Purified libraries were quantified by Qubit® 2.0 Fluorometer (Life Technologies, Waltham, MA, USA) and validated by Agilent 2100 bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) to confirm the insert size and calculate the mole concentration. Cluster was generated by cBot with the library diluted to 10 pM and then were sequenced on the Illumina HiSeq 2500 (Illumina, San Diego, CA, USA). The library construction and sequencing was performed at Shanghai Biotechnology Corporation.

2.4. Bioinformatics Analysis

Bioinformatics analysis was performed as described by Guillaume et al. [21] but with slight modifications. Briefly, library adapters were trimmed from the sequences in the fastq file and the reads below 50 bp were discarded. The host genome was removed by mapping to the Culex genomes with Bowtie2 [22]. The sequence reads that did not match the host genome were retained and subsequently used for de novo assembly in the “Trinity” tool [23]. These de novo assembled contigs were then subjected to basic local alignment search tools (BLASTn and BLASTx). The databases used for BLASTn and BLASTx were built from the viral reference sequences as described in Metavisitor pipeline [21]. An E-value cut-off of 0.01 was used. Following BLAST analysis, accession numbers were parsed to NCBI taxonomy files in R using package “taxonomizr” [24].
Taxonomic classification of identified viruses was done by conducting BLASTn and BLASTx on National Center For Biotechnology Information (NCBI) virus database and the obtained Blast XML files were imported into MEGAN6 program [25] and default parameters used. For comparison of viral content from the three study sites, assembled contigs recognized as viral after BLASTn and BLASTx were further classified into their respective taxonomic families or into a general group of unclassified viruses if not assigned a taxonomic family by the International Committee on Taxonomy of Viruses (ICTV).

2.5. Confirmation of Viral Sequences

Analyzed viral sequences were used to design specific PCR primers to confirm the presence of some selected viruses from the original RNA. PCR amplifications were carried out in a total volume of 25 μL consisting of 2 μL of template (cDNA), 0.25 μM of forward and reverse primers, 12.5 μL of PrimerSTAR Max Premix 2× (TaKaRa, Kusatsu, Japan), and Milli-Q ddH2O to supplement the system. Thermal cycling was initiated with a denaturation step at 94 °C for 5 min; followed by 35 cycles of 94 °C for 30 s, 58–60 °C 30 s and 72 °C for 1 min; and a final extension at 72 °C for 10 min. The PCR products were extracted with a Gel extraction kit (Omega Bio-Tek, Norcross, GA, USA) and sequenced at Tsingke Biological Technology. The primers used are presented in Table S1.

3. Results

3.1. Proportion of Viruses in Host Groups

In Kenya, we did not trap Culex spp. mosquitoes in sites within Kilifi County except sites in Kwale County where we dominantly trapped Culex quinquefasciatus. In China, we dominantly trapped Culex tritaeniorhynchus in all the sites within Hubei and Yunnan provinces. We sequenced 47 mosquito pools from Kwale, 100 pools from Hubei, and 100 pools from Yunnan. Mosquito samples from Kwale, Kenya generated a total of 21,747,508 raw reads. Samples from Hubei and Yunnan provinces, China had 6,714,707 and 36,277,174 raw reads, respectively. The average sequence reads lengths were 150 bp, 125 bp, and 398 bp for Kwale, Hubei and Yunnan respectively (Table 1).
Our results showed that mosquitoes harbored viruses from diverse host groups: vertebrate viruses, invertebrate viruses, plant viruses, mycoviruses, phages, and protozoan viruses as well as groups that were classified as environmental and unidentified viruses (Figure 2). Across the three locations, abundant viruses (proportions > 5%) were those specific to host groups of vertebrates, invertebrates and plants, with exception of Kwale, which in addition had abundant protozoan viruses. Specifically, vertebrate host viruses were the most abundant in Kwale and Hubei. Invertebrate-specific viruses were the most abundant in mosquitoes of Yunnan Province (Figure 2).

3.2. Virus Families

Across the three locations, the mosquitoes harbored a total of 30 virus families and one general group of unclassified viruses (Figure 3). Vertebrate virus families include Adenoviridae, Anelloviridae, Circoviridae, Coronaviridae, Flaviviridae, Herpesviridae, Papillomaviridae, Reoviridae, Retroviridae and Togaviridae. Kwale had 14 virus families, with Papillomaviridae being the most abundant (Figure 3). Hubei had 18 virus families with Herpesviridae being the most abundant and Yunnan had 11 virus families, with the most abundant being unclassified viruses (Figure 3).
In Kwale, vertebrate, invertebrate, plant and protozoan viruses were the most abundant (Figure 2). The virome of Kwale mosquito was comprised of 14 virus families, with Papillomaviridae (57%) being most dominant, followed by Herspesviridae (8.17%). Among the invertebrate specific viruses, Baculoviridae was the notable virus family, whereas Partitiviridae (0.6%) was the leading plant virus. Kwale had a notable proportion of protozoan viruses, which were dominated by Pandoravirus salinus and Pandoravirus dulcis (unclassified viruses). These viruses were not identified in samples from China.
In the Hubei province, vertebrate, invertebrate, plant and unclassified viruses were the most abundant (Figure 2). Herpesviridae (30.16%) was the dominant vertebrate virus in the virome of Hubei followed by Adenoviridae (14%). Nimaviridae (3.6%) was the leading virus family in the invertebrate host group, while Genomoviridae (9.4%) was leading in part of plant viruses from Hubei.
The virome of mosquito in the Yunnan province was dominantly composed of invertebrate-host related viruses (58.7%) (Figure 2) of which the majority were insect-specific viruses that included families Iflaviridae, Dicistroviridae, and Mesoniviridae, and some belonged to the unclassified virus groups. Our results showed that the Wuhan mosquito virus (unclassified viruses) and Culex tritaeniorhynchus rhabdovirus (Rhabdoviridae) were the most common virus species. In addition, Yunnan had a higher proportion of unclassified viruses (27.17%) compared to Kwale and Hubei. Under vertebrate viruses, the family Reoviridae (24.6%) was the most conspicuous group and it was dominated by the Banna virus (33%) and Kadipiro virus (11.6%) (Figure 4).

3.3. Viromes across the Three Locations

In the three study sites, the three most abundant virus families were Flaviviridae (Kwale—1.6%, Hubei—2.7% and Yunnan—2.5%) and Herpesviridae (Kwale—8.2%, Hubei—30%, and Yunnan—0.5%) as well as the group of unclassified viruses (Kwale—9.7%, Hubei—2.3%, and Yunnan—27.14%).
Under the family Flaviviridae, both Hubei (2.7% reads) and Kwale (1.6% reads) had flaviviruses that infected vertebrates and also some insect flaviviruses. Specifically, in Kwale we identified the dengue 2 virus (77 reads) which was further confirmed by PCR. Culex flavivirus (CxFV) was the only insect-specific flavivirus identified in Kwale. Yunnan Province exclusively had insect-specific flaviviruses (2.5%), mainly dominated by the Quang Binh virus (QBV). In Hubei province, flavivirus was similarly dominated by QBV.
Additionally, the family Mimiviridae (protozoan-host) was present in all the three study sites, while Totiviridae (protozoan-host) was present in Yunnan and Kwale mosquitoes. Mycoviruses, bacteriophages, and environmental viruses contributed a significantly small percentage in all the three study sites.

4. Discussion

Emerging studies on mosquito viromes depict a highly diverse community of viruses whose existence, role, and transmission are less understood. Our results are thus consistent with previous studies as they revealed a rich diversity of both known pathogenic groups of viruses to humans and animals as well as viruses specific to other organisms (plants, protozoans, and invertebrates). Given that mosquito trappings in the three locations in Kenya and China were carried out at different times, our findings should be interpreted with that in mind. From our findings, Kwale, Kenya and Hubei, China (Figure 1) were comparable in virome profiles in that both were dominated by vertebrate viruses while virome profiles in Yunnan were dominated by invertebrate viruses, specifically insect-specific viruses (Figure 2). It thus appears that Kwale and Hubei had a high presence of vertebrates in the environment coupled by the probability that the sampled mosquitoes were mostly gravid and questing for blood meal. Interestingly, this was the first mosquito virome study in Kenya and provides a preliminary status.
Across the three locations, the mosquito viromes harbored a total of 30 virus families and one general group of unclassified viruses (Figure 3). The virus families included those specific to vertebrates, invertebrates, plants, protozoa, and bacteriophages. Since the three study sites in Kenya and China (Figure 1) were ecologically rich in mammalian and plant diversity (and were in close proximity to human settlements), the diversity was likely shaped by factors such as sex, whether female mosquitoes were gravid, the indiscriminate mosquito diet range (that includes plants and blood), and the micro-ecology supporting diversity of microorganisms, as well as contaminants in the environmental substrates they used for resting [26,27,28,29,30]. In the present study, we identified vertebrate viruses known to naturally infect mammals and replicate in mammalian cell lines, including Adenoviridae, Anelloviridae, Circoviridae, Coronaviridae, Flaviviridae, Herpesviridae, Papillomaviridae, Reoviridae, Retroviridae, Togaviridae and Rhabdoviridae. These virus families are of medical and veterinary importance and some have been identified as part of mosquito viromes elsewhere [3,31]. Across the three study sites, Hubei had the highest number of vertebrate-specific viruses, followed by Kwale, and the least were found in Yunnan (Figure 3), a variation that could be attributed to variation in the number of sampled blood-fed female mosquitoes. We did not quantify the sampled mosquitoes by sex, which we consider as one caveat of this study.
Among the vertebrate virus group in Yunnan (based on the top ten most abundant virus species), Banna virus and Kadipiro virus (family Reoviridae) had the most abundant viral reads in this set (Figure 4). These two are arboviruses of public health interest because they can infect vertebrates, including humans. Banna virus (BAV) and Kadipiro virus both contain double-stranded RNA with 12 segments and they were both classified under the genus Seadonavirus by the ICTV in 2016 [32]. Other known seadonaviruses are Liao ning virus [33], Balaton virus [34] and Mangshi virus [35]. In previous studies, BAV has been isolated within China from cattle, pigs, ticks, mosquitoes [36,37,38], and from the sera of febrile patients [39]. Although BAV has not been successfully propagated in mammalian cell lines, its seropositivity implies it can replicate in humans where it is speculated to cause encephalitis. BAV distribution has been shown to extend from tropical zones to the north temperate climate regions of Asia. Outside China, Banna virus has been isolated in Vietnam [40], Indonesia [41], and Inner Mongolia [42]. Wang and his team have also isolated BAV in Yunnan province, at China’s border with Laos and Myanmar and they suggest that BAV together with other lesser-known arboviruses might be silently circulating and infecting humans in the region [17]. Kadipiro virus was first isolated in 1981 from Culex fuscocephalus mosquitoes captured in Java, Indonesia. In China, the first isolation of Kadipiro virus was from Culex tritaeniorhynchus, Anopheles sinensis, and Armigeres subalbatus mosquitoes, all from the Yunnan province [43]. Similarly, as with BAV, propagation of Kadipiro virus in mammalian cell lines has not been reported yet.
Our results showed that viruses in the family Herpesviridae represent a significant component of the Culex mosquito virome, given that it was present in all the locations in China and Kenya. In addition, the virus family was most abundant in Hubei, coming second in Kwale, and was least prevalent in Yunnan (Figure 3). Herpesviruses are endogenous to a wide variety of vertebrate and invertebrate organisms and they infect a wide range of organisms, although through their ability to modulate host immunity, they are able to cause lifelong latent infections [44,45]. Our study identified several herpesviruses, with human betaherpesvirus 7 being dominant in Kwale while cercopithecine betaherpesvirus 5 was dominant in Hubei (Figure 4). Some of the human herpesviruses cause pathology in humans [46]. Out of the numerous herpesviruses latent in non-human primates only cercopithecine herpesvirus 1 is known to cause fatal illness in humans [47,48]. The presence and abundance of these herpesviruses suggest human and cercopithecine as easily available blood sources for mosquitoes in Hubei and Kwale. Further, mosquitoes are suggested to play a key role in the transmission of herpesviruses, such as the human herpesvirus-8 in humans [49,50,51]. Previous metagenome studies on mosquitoes [4,31], bats [52], and the most recent deep transcriptome sequencing of more than 220 invertebrate species [53] have all reported presence of herpesviruses from the three taxonomic subfamilies Alphaherpesvirinae, Betaherpesvirinae and Gammaherpesvirinae [45,54].
Further, our results showed that Papillomaviridae was a specifically dominant virus family in Kwale, while it was absent in Yunnan (Figure 3). Papillomaviruses are known to infect virtually all amniote vertebrates but they do not trigger a strong immune response [55,56]. In their natural hosts and their related species, papillomaviruses cause papillomas (benign tumors) in specific body sites such as the skin and mucosal epithelia. In both Kwale and Hubei, we identified alphapapillomavirus 9 (Figure 4), a species group associated with human cervical cancer which is responsible for 75% of all cases of cervical cancer worldwide [57]. However, overall, 80% of healthy human skin latently harbors diverse papillomaviruses [58]. It is therefore plausible that during blood meal from humans, the mosquito ingested these viruses.
Although not among the top ten virus species (Figure 4), but rather due to its public health importance, we detected dengue virus serotype 2 (flavivirus; Flaviviridae) in the Kwale mosquito virome. The virus causes variable symptoms ranging from mild febrile illness to the more severe dengue hemorrhagic fever, which can progress to a fatal form, dengue shock syndrome (DSS). Aedes aegyptii and Aedes albopictus are the two competent vectors of the virus worldwide [59,60,61]; hence, Culex mosquitoes ingested the virus from viremic hosts. The Kenyan coastal towns, which include Kwale, have experienced intermittent outbreaks due to dengue 2 in the past and presence of the virus suggests continuous low-level circulation.
It is worth noting that Culex tritaeniorhynchus is the main vector of the endemic Japanese encephalitis virus (JEV) and other zoonotic viruses like Banna virus and Kadipiro virus in Yunnan, Hubei and other provinces of China [17,62,63,64]. Although West Nile virus (WNV) is not endemic in China, it has been severally reported in the Xinjiang province [65,66,67]. In Kenya, there has been sporadic isolations of WNV in Culex univittatus from the Rift Valley [68], wild birds in the Tana River and Garissa Counties [69] and in ticks from north-eastern Kenya [70], We did not detect these important public health viruses from our mosquito viromes, which suggests that the sampled mosquitoes sourced blood meals from non-infected hosts. Perhaps a longer period of mosquito sampling is necessary to pick the virus especially if prevalence in the vertebrate hosts is low. Nevertheless, from our top ten vertebrate viruses, we inferred that the mosquitoes fed from a variety of vertebrate hosts. (Figure 4). Discovery of insect-specific viruses is on the rise, yet their role and impact on arbovirus transmission or potential pathogenicity is equivocal. However, they constitute a significant proportion in most of the mosquito viromes. Our study identified the virus families Iflaviridae, Dicistroviridae and Mesoniviridae as the most common, which is consistent with virome diversity and proportions of Culex spp. in Zambezi, Mozambique and Hubei, China [4,13]. Interestingly, previous studies showed that both Culex mosquito viromes in Zambezi and Hubei were dominated (up to 88%) by the insect-specific viruses, Iflaviridae [4,13]. In addition, we found that the Yunnan and Hubei mosquito viromes also had the Quang Binh virus (QBV), which is an invertebrate flavivirus. Still, Hubei was defined by the abundant Nimaviridae, that is usually composed of a single genus, Whispovirus, known to be of global occurrence and fatally infects crustaceans [71]. In contrast, the abundant insect-specific virus in Kwale virome was the Baculoviridae, a virus family that is known to infect diverse invertebrates but has also been reported to cause epizootics in field populations of Culex spp. [72]. Species from this family are being tested for use as mosquito control [73]. In addition, we identified Culex flavivirus (CxFV), a virus that has previously been isolated across the world.
In the present study, we identified a notable proportion of plant-specific viruses across viromes in the three locations (Figure 2). We identified the family Tymoviridae—a group that constitutes globally distributed viruses—that causes significant losses in crops due to plant diseases [74]. Mosquitoes are not known to be vectors of plant-specific viruses even though due to feeding on plant sap and nectar, they ingest these viruses. Recently, a plant-specific virus, Culex originated Tymoviridae-like virus (CuTLV) (Tymoviridae), isolated from Culex spp. in China, showed cytopathic effect on mosquito cell lines, implying a potential for being vectored by mosquitoes [75].
In summary, mosquito virome studies are very few, implying that there is an enormous diversity of viruses yet to be uncovered. Our study was preliminary but it revealed that at a larger habitat scale, geolocation as well as the sampling time period are important factors that could enhance recovery of a greater diversity of mosquito viruses. We determined that each virome was unique in terms of virus proportions partly influenced by type of ingested meals (blood, nectar, plant sap, environmental substrates). Although vertebrate and other viruses ingested from other hosts are significant components of the virome, we suggest that since insect-specific viruses occur naturally in their hosts, they should be the basis for defining the mosquito virome. For instance, in this study, the dominant insect-specific viruses in Kwale, Hubei, and Yunnan were Baculoviridae, Nimaviridae and Iflaviridae, respectively, and hence such dominant viruses might have a role in the component community of the host. Our study is preliminary but due to emerging theories that suggest that manipulation of specific components of mosquito viromes could yield potential agents for control and prevention of pathogenic arboviruses and also the global emerging trends of mosquito-borne diseases, more studies are still necessary to uncover the global patterns of composition and diversity of mosquito viromes.

Supplementary Materials

The following are available online at www.mdpi.com/1999-4915/10/1/30/s1, Table S1: PCR primers used in validation of sequenced metagenomic data.

Acknowledgments

We are grateful to the entomology section, National Museums of Kenya and Hubei Provincial Center for Disease Control and Prevention (Hubei CDC) for assistance in sampling. This work was supported by the Ministry of Science and Technology of China (Science and Technology Basic Work Program 2013FY113500) and Sino-Africa Joint Research Center, Chinese Academy of Sciences (SAJC201605).

Author Contributions

Zhiming Yuan, Han Xia and Evans Atoni conceived and designed the experiments; Evans Atoni and Yujuan Wang performed the experiments; Evans Atoni, Yujuan Wang, Han Xia, Samuel Karungu, Ali Zohaib and Cecilia Waruhiu analyzed the data; Han Xia, Vincent Obanda, Bernard Agwanda, Morris Mutua and Zhiming Yuan contributed reagents/materials/analysis tools; Evans Atoni, Yujuan Wang, Vincent Obanda, Han Xia and Zhiming Yuan wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brinkmann, A.; Nitsche, A.; Kohl, C. Viral Metagenomics on Blood-Feeding Arthropods as a Tool for Human Disease Surveillance. Int. J. Mol. Sci. 2016, 17, 1743. [Google Scholar] [CrossRef] [PubMed]
  2. Fauver, J.R.; Grubaugh, N.D.; Krajacich, B.J.; Weger-Lucarelli, J.; Lakin, S.M.; Fakoli, L.S.; Bolay, F.K.; Diclaro, J.W.; Dabiré, K.R.; Foy, B.D.; et al. West African Anopheles Gambiae Mosquitoes Harbor a Taxonomically Diverse Virome Including New Insect-Specific Flaviviruses, Mononegaviruses, and Totiviruses. Virology 2016, 498, 288–299. [Google Scholar] [CrossRef] [PubMed]
  3. Frey, K.G.; Biser, T.; Hamilton, T.; Santos, C.J.; Pimentel, G.; Mokashi, V.P.; Bishop-Lilly, K.A. Bioinformatic Characterization of Mosquito Viromes within the Eastern United States and Puerto Rico: Discovery of Novel Viruses. Evol. Bioinform. 2016, 12, 1–12. [Google Scholar] [CrossRef] [PubMed]
  4. 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] [PubMed]
  5. Delwart, E.L. Viral Metagenomics. Rev. Med. Virol. 2007, 17, 115–131. [Google Scholar] [CrossRef] [PubMed]
  6. Moormann, A.M.; Snider, C.J.; Chelimo, K. The Company Malaria Keeps: How Co-Infection with Epstein-Barr Virus Leads to Endemic Burkitt Lymphoma. Curr. Opin. Infect. Dis. 2011, 24, 435–441. [Google Scholar] [CrossRef] [PubMed]
  7. Grubaugh, N.D.; Sharma, S.; Krajacich, B.J.; Fakoli, L.S., III; Bolay, F.K.; Diclaro, J.W., II; Johnson, W.E.; Ebel, G.D.; Foy, B.D.; Brackney, D.E. Xenosurveillance: A Novel Mosquito-Based Approach for Examining the Human-Pathogen Landscape. PLoS Negl. Trop. Dis. 2015, 9, e0003628. [Google Scholar] [CrossRef] [PubMed]
  8. Stollar, V.; Thomas, V.L. An Agent in the Aedes Aegypti Cell Line (Peleg) Which Causes Fusion of Aedes albopictus Cells. Virology 1975, 64, 367–377. [Google Scholar] [CrossRef]
  9. 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]
  10. Bolling, B.G.; Olea-Popelka, F.J.; Eisen, L.; Moore, C.G.; Blair, C.D. Transmission Dynamics of an Insect-Specific Flavivirus in a Naturally Infected Culex pipiens Laboratory Colony and Effects of Co-Infection on Vector Competence for West Nile Virus. Virology 2012, 427, 90–97. [Google Scholar] [CrossRef] [PubMed]
  11. Burivong, P.; Pattanakitsakul, S.-N.; Thongrungkiat, S.; Malasit, P.; Flegel, T.W. Markedly Reduced Severity of Dengue Virus Infection in Mosquito Cell Cultures Persistently Infected with Aedes albopictus Densovirus (AalDNV). Virology 2004, 329, 261–269. [Google Scholar] [CrossRef] [PubMed]
  12. Karpf, A.R.; Lenches, E.; Strauss, E.G.; Strauss, J.H.; Brown, D.T. Superinfection Exclusion of Alphaviruses in Three Mosquito Cell Lines Persistently Infected with Sindbis Virus. J. Virol. 1997, 71, 7119–7123. [Google Scholar] [PubMed]
  13. Cholleti, H.; Hayer, J.; Abilio, A.P.; Mulandane, C.; Verner-carlsson, J.; Falk, K.I.; Fafetine, J.M.; Berg, M.; Blomström, A. Discovery of Novel Viruses in Mosquitoes from the Zambezi Valley of Mozambique. PLoS ONE 2016, 11, e0162751. [Google Scholar] [CrossRef] [PubMed]
  14. Cornel, A.J.; Hunt, R.H. Aedes albopictus in Africa? First Records of Live Specimens in Imported Tires in Cape Town. J. Am. Mosq. Control Assoc. 1991, 7, 107–108. [Google Scholar] [PubMed]
  15. Ge, X.; Li, Y.; Yang, X.; Zhang, H.; Zhou, P.; Zhang, Y.; Shi, Z. Metagenomic Analysis of Viruses from Bat Fecal Samples Reveals Many Novel Viruses in Insectivorous Bats in China. J. Virol. 2012, 86, 4620–4630. [Google Scholar] [CrossRef] [PubMed]
  16. He, B.; Li, Z.; Yang, F.; Zheng, J.; Feng, Y.; Guo, H.; Li, Y.; Wang, Y.; Su, N.; Zhang, F.; et al. Virome Profiling of Bats from Myanmar by Metagenomic Analysis of Tissue Samples Reveals More Novel Mammalian Viruses. PLoS ONE 2013, 8, e61950. [Google Scholar] [CrossRef]
  17. Wang, J.; Zhang, H.; Sun, X.; Fu, S.; Wang, H.; Feng, Y.; Wang, H.; Tang, Q.; Liang, G.D. Distribution of Mosquitoes and Mosquito-Borne Arboviruses in Yunnan Province Near the China-Myanmar-Laos Border. Am. J. Trop. Med. Hyg. 2011, 84, 738–746. [Google Scholar] [CrossRef] [PubMed]
  18. 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] [PubMed]
  19. Edwards, F. Mosquitoes of the Ethiopian Region III—Culicine Adults and Pupae | Mosquito Taxonomic Inventory; British Museum (Natural History): London, UK, 1941. [Google Scholar]
  20. Hopkins, G. Mosquitoes of the Ethiopian Region I. Larval Bionomics of Mosquitoes and Taxonomy of Culicine Larvae; British Museum (Natural History): London, UK, 1952. [Google Scholar]
  21. Carissimo, G.; van den Beek, M.; Vernick, K.D.; Antoniewski, C.; Schweer, T.; Yarza, P.; Schneider, D. Metavisitor, a Suite of Galaxy Tools for Simple and Rapid Detection and Discovery of Viruses in Deep Sequence Data. PLoS ONE 2017, 12, e0168397. [Google Scholar] [CrossRef] [PubMed]
  22. Langmead, B.; Salzberg, S.L. Fast Gapped-Read Alignment with Bowtie 2. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed]
  23. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-Length Transcriptome Assembly from RNA-Seq Data without a Reference Genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  24. GitHub. Taxonomizr. Available online: https://github.com/sherrillmix/taxonomizr (accessed on 23 August 2017).
  25. Huson, D.H.; Auch, A.F.; Qi, J.; Schuster, S.C. MEGAN Analysis of Metagenomic Data. Genome Res. 2007, 17, 377–386. [Google Scholar] [CrossRef] [PubMed]
  26. Molaei, G.; Andreadis, T.G.; Armstrong, P.M.; Dennett, J.A.; Real, S.V.; Sargent, C.; Bala, A.; Randle, Y.; Guzman, H.; Wuithiranyagool, T.; et al. Host Feeding Pattern of Culex quinquefasciatus (Diptera: Culicidae) and Its Role in Transmission of West Nile Virus in Harris Country, Texas. Annu. J. Trop. Med. Hyg. 2007, 77, 73–81. [Google Scholar]
  27. Burkett-Cadena, N.D.; Bingham, A.M.; Porterfield, C.; Unnasch, T.R. Innate Preference or Opportunism: Mosquitoes Feeding on Birds of Prey at the Southeastern Raptor Center. J. Vector Ecol. 2014, 39, 21–31. [Google Scholar] [CrossRef] [PubMed]
  28. Börstler, J.; Jöst, H.; Garms, R.; Krüger, A.; Tannich, E.; Becker, N.; Schmidt-Chanasit, J.; Lühken, R. Host-Feeding Patterns of Mosquito Species in Germany. Parasit. Vectors 2016, 9, 318. [Google Scholar] [CrossRef] [PubMed]
  29. Sawabe, K.; Isawa, H.; Hoshino, K.; Sasaki, T.; Roychoudhury, S.; Higa, Y.; Kasai, S.; Tsuda, Y.; Nishiumi, I.; Hisai, N.; et al. Host-Feeding Habits of Culex pipiens and Aedes albopictus (Diptera: Culicidae) Collected at the Urban and Suburban Residential Areas of Japan. J. Med. Entomol. 2010, 47, 442–450. [Google Scholar] [CrossRef] [PubMed]
  30. Takken, W.; Verhulst, N.O. Host Preferences of Blood-Feeding Mosquitoes. Annu. Rev. Entomol. 2013, 58, 433–453. [Google Scholar] [CrossRef] [PubMed]
  31. Ng, T.F.F.; Willner, D.L.; Lim, Y.W.; Schmieder, R.; Chau, B.; Nilsson, C.; Anthony, S.; Ruan, Y.; Rohwer, F.; Breitbart, M. Broad Surveys of DNA Viral Diversity Obtained through Viral Metagenomics of Mosquitoes. PLoS ONE 2011, 6, e20579. [Google Scholar] [CrossRef] [PubMed]
  32. International Committee on Taxonomy of Viruses (ICTV). Virus Taxonomy: The Classification and Nomenclature of Viruses the Online (10th) Report of the ICTV; ICTV: Budapest, Hungary, 2016. [Google Scholar]
  33. Attoui, H.; Mohd Jaafar, F.; Belhouchet, M.; Tao, S.; Chen, B.; Liang, G.; Tesh, R.B.; de Micco, P.; de Lamballerie, X. Liao Ning Virus, a New Chinese Seadornavirus That Replicates in Transformed and Embryonic Mammalian Cells. J. Gen. Virol. 2006, 87, 199–208. [Google Scholar] [CrossRef] [PubMed]
  34. Reuter, G.; Boros, Á.; Delwart, E.; Pankovics, P. Novel Seadornavirus (Family Reoviridae) Related to Banna Virus in Europe. Arch. Virol. 2013, 158, 2163–2167. [Google Scholar] [CrossRef] [PubMed]
  35. Wang, J.; Li, H.; He, Y.; Zhou, Y.; Meng, J.; Zhu, W.; Chen, H.; Liao, D.; Man, Y. Isolation and Genetic Characterization of Mangshi Virus: A Newly Discovered Seadornavirus of the Reoviridae Family Found in Yunnan Province, China. PLoS ONE 2015, 10, e0143601. [Google Scholar] [CrossRef] [PubMed]
  36. Li, Q.P. First Isolation of New Orbivirus (BANNA) from Ticks and Infected Cattle Sera in Xinjiang. Endem. Dis. Bull. 1992, 7, 64–69. [Google Scholar]
  37. Xu, P.; Wang, Y.; Zuo, J.; Che, Y.; Peng, H.; Huang, Z.; Tang, G.; Xu, P.; Institute of Virology; et al. Recovery of the Same Type of Virus as Human New Orbivirus from Sera of Cattles and Pigs Collected in Yunnan Province. Chin. J. Virol. 1990, 6, 327–331. [Google Scholar]
  38. Xu, P.; Wang, Y.; Zuo, J.; Lin, J.; Xu, P. New Orbiviruses Isolated from Patients with Unknown Fever and Encephalitis in Yunnan Province. Chin. J. Virol. 1990, 6, 27–33. [Google Scholar]
  39. Tao, S.; Chen, B. Studies of Coltivirus in China. Chin. Med. J. (Engl.) 2005, 118, 581–586. [Google Scholar] [PubMed]
  40. Nabeshima, T.; Nga, P.T.; Guillermo, P.; Parquet, M.; Yu, F.; Thuy, N.T.; Trang, B.M.; Hien, N.T.; Nam, V.S.; Inoue, S.; et al. Isolation and Molecular Characterization of Banna Virus from Mosquitoes, Vietnam. Emerg. Infect. Dis. 2008, 14, 1276–1279. [Google Scholar] [CrossRef] [PubMed]
  41. Brown, S.E.; Gorman, B.M.; Tesh, R.B.; Knudson, D.L. Coltiviruses Isolated from Mosquitoes Collected in Indonesia. Virology 1993, 196, 363–367. [Google Scholar] [CrossRef] [PubMed]
  42. Cao, Y.; Fu, S.; Tian, Z.; He, Y.; Wang, H.; Wang, H.; Yang, H.; Tao, B.; Liang, G. Isolation and Identification of Banna Virus from Mosquito for the First Time in Inner Mongolia. Zhonghua Shi Yan He Lin Chuang Bing Du Xue Za Zhi 2009, 23, 106–108. [Google Scholar] [PubMed]
  43. Sun, X.; Meng, W.; Fu, S.; Feng, Y.; Zhai, Y.; Wang, J.; Wang, H.; Lv, X.; Liang, G. The First Report of Kadipiro Virus Isolation in China. Chin. J. Virol. 2009, 25, 173–177. [Google Scholar]
  44. Davison, A.J. Overview of Classification; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  45. Roizman, B.; Carmichael, L.E.; Deinhardt, F.; de-The, G.; Nahmias, A.J.; Plowright, W.; Rapp, F.; Sheldrick, P.; Takahashi, M.; Wolf, K. Herpesviridae. Definition, Provisional Nomenclature, and Taxonomy. The Herpesvirus Study Group, the International Committee on Taxonomy of Viruses. Intervirology 1981, 16, 201–217. [Google Scholar] [CrossRef] [PubMed]
  46. Campadelli-Fiume, G.; Mirandola, P.; Menotti, L. Human Herpesvirus 6: An Emerging Pathogen. Emerg. Infect. Dis. 1999, 5, 353–366. [Google Scholar] [CrossRef] [PubMed]
  47. Elmore, D.; Eberle, R. Monkey B Virus (Cercopithecine Herpesvirus 1). Comp. Med. 2008, 58, 11–21. [Google Scholar] [PubMed]
  48. Takashima, Y.; Otsuka, H. Pathogenesis of Animal Herpesviruses to Human. Nihon Rinsho 2000, 58, 957–961. [Google Scholar] [PubMed]
  49. McHardy, J.; Williams, E.H.; Geser, A.; de-Thé, G.; Beth, E.; Giraldo, G. Endemic Kaposi’s Sarcoma: Incidence and Risk Factors in the West Nile District of Uganda. Int. J. Cancer 1984, 33, 203–212. [Google Scholar] [CrossRef] [PubMed]
  50. Coluzzi, M.; Manno, D.; Guzzinati, S.; Tognazzo, S.; Zambon, P.; Arcà, B.; Costantini, C.; Ascoli, V. The Bloodsucking Arthropod Bite as Possible Cofactor in the Transmission of Human Herpesvirus-8 Infection and in the Expression of Kaposi’s Sarcoma Disease. Parassitologia 2002, 44, 123–129. [Google Scholar] [PubMed]
  51. Ascoli, V.; Manno, D.; Coluzzi, M. Geographic Variation in Human Herpesvirus 8 Seroprevalence and Possible Association with Exposure to Bites from Blood-Sucking Arthropods. J. Infect. Dis. 2006, 194, 401–402. [Google Scholar] [CrossRef] [PubMed]
  52. Wodak, E.; Richter, S.; Bagó, Z.; Revilla-Fernández, S.; Weissenböck, H.; Nowotny, N.; Winter, P. Detection and Molecular Analysis of West Nile Virus Infections in Birds of Prey in the Eastern Part of Austria in 2008 and 2009. Vet. Microbiol. 2011, 149, 358–366. [Google Scholar] [CrossRef] [PubMed]
  53. Shi, M.; Lin, X.; Tian, J.; Chen, L.; Chen, X.; Li, C.; Qin, X.; Li, J.; Cao, J.; Eden, J.; et al. Article Redefining the Invertebrate RNA Virosphere. Nature 2016, 540, 539–543. [Google Scholar] [CrossRef] [PubMed]
  54. Davison, A.J. Herpesvirus Systematics. Vet. Microbiol. 2010, 143, 52–69. [Google Scholar] [CrossRef] [PubMed]
  55. Bravo, I.G.; Félez-Sánchez, M. Papillomaviruses: Viral Evolution, Cancer and Evolutionary Medicine. Evol. Med. Public Health 2015, 2015, 32–51. [Google Scholar] [CrossRef] [PubMed]
  56. De Villiers, E.-M.; Fauquet, C.; Broker, T.R.; Bernard, H.-U.; zur Hausen, H. Classification of Papillomaviruses. Virology 2004, 324, 17–27. [Google Scholar] [CrossRef] [PubMed]
  57. Muñoz, N.; Bosch, F.X.; de Sanjosé, S.; Herrero, R.; Castellsagué, X.; Shah, K.V.; Snijders, P.J.F.; Meijer, C.J.; International Agency for Research on Cancer Multicenter Cervical Cancer Study Group. Epidemiologic Classification of Human Papillomavirus Types Associated with Cervical Cancer. N. Engl. J. Med. 2003, 348, 518–527. [Google Scholar] [CrossRef] [PubMed]
  58. Antonsson, A.; Forslund, O.; Ekberg, H.; Sterner, G.; Hansson, B.G. The Ubiquity and Impressive Genomic Diversity of Human Skin Papillomaviruses Suggest a Commensalic Nature of These Viruses. J. Virol. 2000, 74, 11636–11641. [Google Scholar] [CrossRef] [PubMed]
  59. Lutomiah, J.; Barrera, R.; Makio, A.; Mutisya, J.; Koka, H.; Owaka, S.; Koskei, E.; Nyunja, A.; Eyase, F.; Coldren, R.; et al. Dengue Outbreak in Mombasa City, Kenya, 2013–2014: Entomologic Investigations. PLoS Negl. Trop. Dis. 2016, 10, e0004981. [Google Scholar] [CrossRef] [PubMed]
  60. Ngoi, C.N.; Price, M.A.; Fields, B.; Bonventure, J.; Ochieng, C.; Mwashigadi, G.; Hassan, A.S.; Thiong, A.N.; Micheni, M.; Mugo, P.; et al. Dengue and Chikungunya Virus Infections among Young Febrile Adults Evaluated for Acute HIV-1 Infection in Coastal Kenya. PLoS ONE 2016, 11, e0167508. [Google Scholar] [CrossRef] [PubMed]
  61. Sang, R.C. Dengue in Africa. J. Med. Entomol. 2005, 2, 3–6. [Google Scholar]
  62. Sun, X.; Fu, S.; Gong, Z.; Ge, J.; Meng, W.; Feng, Y.; Wang, J.; Zhai, Y.; Wang, H.; Nasci, R.; et al. Distribution of Arboviruses and Mosquitoes in Northwestern Yunnan Province, China. Vector-Borne Zoonotic Dis. 2009, 9, 623–630. [Google Scholar] [CrossRef] [PubMed]
  63. Liu, H.; Lu, H.-J.; Liu, Z.-J.; Jing, J.; Ren, J.-Q.; Liu, Y.-Y.; Lu, F.; Jin, N.-Y. Japanese Encephalitis Virus in Mosquitoes and Swine in Yunnan Province, China 2009–2010. Vector-Borne Zoonotic Dis. 2013, 13, 41–49. [Google Scholar] [CrossRef] [PubMed]
  64. Feng, Y.; Zhang, H.L.; Yang, W.H.; Zhang, Y.Z.; Huang, L.J.; Deng, S.Z.; Sun, Y.J.; Yang, D.J.; Zhou, J.H. Molecular Epidemiology of Japanese Encephalitis Viruses Isolated in Yunnan Province, 1977–2010. Zhonghua Liu Xing Bing Xue Za Zhi 2016, 37, 1519–1525. [Google Scholar] [PubMed]
  65. Lu, Z.; Fu, S.-H.; Cao, L.; Tang, C.-J.; Zhang, S.; Li, Z.-X.; Tusong, M.; Yao, X.-H.; Zhang, H.-L.; Wang, P.-Y.; et al. Human Infection with West Nile Virus, Xinjiang, China, 2011. Emerg. Infect. Dis. 2014, 20, 1421–1423. [Google Scholar] [CrossRef] [PubMed]
  66. Li, X.-L.; Fu, S.-H.; Liu, W.-B.; Wang, H.-Y.; Lu, Z.; Tong, S.-X.; Li, Z.-X.; Nasci, R.S.; Kosoy, O.; Cui, Y.; et al. West Nile Virus Infection in Xinjiang, China. Vector-Borne Zoonotic Dis. 2013, 13, 131–133. [Google Scholar] [CrossRef] [PubMed]
  67. Cao, L.; Fu, S.; Lv, Z.; Tang, C.; Cui, S.; Li, X.; Gao, X.; Li, M.; Cao, Y.; Lei, W.; et al. West Nile Virus Infection in Suspected Febrile Typhoid Cases in Xinjiang, China. Emerg. Microbes Infect. 2017, 6, e41. [Google Scholar] [CrossRef] [PubMed]
  68. Miller, B.R.; Nasci, R.S.; Godsey, M.S.; Savage, H.M.; Lutwama, J.J.; Lanciotti, R.S.; Peters, C.J. First Field Evidence for Natural Vertical Transmission of West Nile Virus in Culex Univittatus Complex Mosquitoes from Rift Valley Province, Kenya. Am. J. Trop. Med. Hyg. 2000, 62, 240–246. [Google Scholar] [CrossRef] [PubMed]
  69. Nyamwaya, D.; Wang’ondu, V.; Amimo, J.; Michuki, G.; Ogugo, M.; Ontiri, E.; Sang, R.; Lindahl, J.; Grace, D.; Bett, B. Detection of West Nile Virus in Wild Birds in Tana River and Garissa Counties, Kenya. BMC Infect. Dis. 2016, 16, 696. [Google Scholar] [CrossRef] [PubMed]
  70. Lwande, O.W.; Venter, M.; Lutomiah, J.; Michuki, G.; Rumberia, C.; Gakuya, F.; Obanda, V.; Tigoi, C.; Odhiambo, C.; Nindo, F.; et al. Whole Genome Phylogenetic Investigation of a West Nile Virus Strain Isolated from a Tick Sampled from Livestock in North Eastern Kenya. Parasites Vectors 2014, 7, 542. [Google Scholar] [CrossRef] [PubMed]
  71. Leu, J.H.; Yang, F.; Zhang, X.; Xu, X.; Kou, G.H.; Lo, C.F. Whispovirus. Curr. Top. Microbiol. Immunol. 2009, 328, 197–227. [Google Scholar] [PubMed]
  72. White, S.E.; Becnel, J.J.; Undeen, A.H.; Rotstein, M.J.; Moser, B.A.; Cockburn, A.; Fukuda, T. Epizootiology and Transmission of a Newly Discovered Baculovirus from the Mosquitoes Culex nigripalpus and C. quinquefasciatus. J. Gen. Virol. 2001, 82, 275–282. [Google Scholar] [CrossRef]
  73. Becnel, J.J. Prospects for the Mosquito Baculovirus Cuninpv as a Tool for Mosquito Control. J. Am. Mosq. Control Assoc. 2006, 22, 523–526. [Google Scholar] [CrossRef]
  74. Stephan, D.; Siddiqua, M.; Ta Hoang, A.; Engelmann, J.; Winter, S.; Maiss, E. Complete Nucleotide Sequence and Experimental Host Range of Okra Mosaic Virus. Virus Genes 2008, 36, 231–240. [Google Scholar] [CrossRef] [PubMed]
  75. Wang, L.; Lv, X.; Zhai, Y.; Fu, S.; Wang, D.; Rayner, S.; Tang, Q.; Liang, G. Genomic Characterization of a Novel Virus of the Family Tymoviridae Isolated from Mosquitoes. PLoS ONE 2012, 7, e39845. [Google Scholar] [CrossRef] [PubMed]
Figure 1. World geographical map outlining the geo-location of Kenya and China. Kenya is outlined in orange in the African continent while China is outlined in green in the Asian continent. The red stars represent sampling sites.
Figure 1. World geographical map outlining the geo-location of Kenya and China. Kenya is outlined in orange in the African continent while China is outlined in green in the Asian continent. The red stars represent sampling sites.
Viruses 10 00030 g001
Figure 2. Percentage proportions of viruses in different host groups from mosquitos in Kwale, Kenya, and the Hubei and Yunnan provinces of China.
Figure 2. Percentage proportions of viruses in different host groups from mosquitos in Kwale, Kenya, and the Hubei and Yunnan provinces of China.
Viruses 10 00030 g002
Figure 3. Stacked bar graph of all virus families (color coded) identified in the viromes of the Culex mosquito in Kwale county, Kenya and provinces of Hubei and Yunnan, China, 2014–2017. Only virus families with reads more than 30 were included in the figure.
Figure 3. Stacked bar graph of all virus families (color coded) identified in the viromes of the Culex mosquito in Kwale county, Kenya and provinces of Hubei and Yunnan, China, 2014–2017. Only virus families with reads more than 30 were included in the figure.
Viruses 10 00030 g003
Figure 4. Abundance of the 10 most common virus species in Culex mosquito viromes in Kwale, Hubei, and Yunnan.
Figure 4. Abundance of the 10 most common virus species in Culex mosquito viromes in Kwale, Hubei, and Yunnan.
Viruses 10 00030 g004
Table 1. Summary of the sampling sites, mosquito species and sequenced data.
Table 1. Summary of the sampling sites, mosquito species and sequenced data.
CountryLocationCoordinatesMosquito SpeciesNo. of Sequenced ReadsAverage Length of Reads
KenyaKwale4.46057° S, 39.47795° ECulex quinquefasciatus21,747,508150 bp
ChinaHubei30.8843° N, 112.5923° ECulex tritaeniorhynchus6,714,707125 bp
Yunnan24.9756° N, 101.4848° ECulex tritaeniorhynchus36,277,174398 bp

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/v10010030

AMA Style

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(1):30. https://doi.org/10.3390/v10010030

Chicago/Turabian Style

Atoni, Evans, Yujuan Wang, Samuel Karungu, Cecilia Waruhiu, Ali Zohaib, Vincent Obanda, Bernard Agwanda, Morris Mutua, Han Xia, and Zhiming Yuan. 2018. "Metagenomic Virome Analysis of Culex Mosquitoes from Kenya and China" Viruses 10, no. 1: 30. https://doi.org/10.3390/v10010030

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