Skip to main content

Co-circulation of Aedes flavivirus, Culex flavivirus, and Quang Binh virus in Shanghai, China

Abstract

Background

With increases in global travel and trade, the spread of arboviruses is undoubtedly alarming. Pathogen detection in field-caught mosquitoes can provide the earliest possible warning of transmission. Insect-specific flavivirus (ISFV) has been first detected in 1991 and documented worldwide in the latest ten years. Although infection with ISFVs is apparently limited to insects, an increase in the infection rate of mosquito-borne flaviviruses may be able to induce cytopathic effects in vertebrate cells during co-infection with other human pathogens. However, little is known whether ISFVs persist in most regions of China.

Methods

During the mosquito activity season in 2016, a surveillance program was carried out to detect ISFVs in mosquitoes in metropolitan Shanghai, China. The presence of ISFVs was randomly tested in different species of mosquitoes using RT-PCR-based and hemi-nested PCR assays, following by the sequencing of PCR products. Sequences from positive pooled samples were compared with those deposited in GenBank. Thereafter, sequences of representative insect flaviviruses were used for further phylogenetic and molecular evolutionary analyses.

Results

Our investigations showed: (1) the presence of Aedes flavivirus (AEFV) in 11/161 pooled samples (nine pools in Songjiang District, one pool in Huangpu District, and one pool in Qingpu District) of Aedes albopictus, (2) the presence of Quang Binh virus (QBV) in 10/195 pooled samples (all in Chongming District) of Culex tritaeniorhynchus; and (3) the presence of Culex flavivirus (CxFV) in 9/228 pooled samples (six pools in Pudong New Area, two pools in Huangpu District, and one pool in Chongming District) of Cx. pipiens. Furthermore, phylogenetic analyses of the gene sequences of envelope proteins indicated that Shanghai CxFV strains belonged to the Asia/USA genotype. The overall maximum likelihood estimation values (and 95% confidence interval) for CxFV, QBV, and AEFV in mosquitoes collected in Shanghai in 2016 were 1.34 (0.66–2.45), 1.65 (0.87–2.85), and 1.51 (0.77–2.70) per 1000, respectively.

Conclusions

This study reveals the presence and the geographical distribution of ISFVs, and determines the genetic variation and the infection rate of ISFVs in Shanghai, China. At least, three insect flaviviruses including ISFVs, AEFV, CxFV, and QBV, co-circulate in this area. To our knowledge, this is the first report of AEFV in China.

Multilingual abstracts

Please see Additional file 1 for translations of the abstract into the five official working languages of the United Nations.

Background

Transmitted by arthropods, arboviruses have become a significant cause of public health issues worldwide, with the potential of leading to unprecedented spread and causing epidemics. Most arboviruses of medical significance belong to three genera: Flavivirus (family Flaviviridae), Alphavirus (family Togaviridae), and Orthobunyavirus (family Peribunyaviridae) [1]. The genus Flavivirus contains close to 80 enveloped, single-stranded, positive-sense RNA viral species that exhibit a wide range of geographic distributions and diverse relationships with their hosts. Most known flaviviruses are associated with human diseases, and are transmitted to vertebrate hosts by mosquitoes or ticks. Their genomes, which vary in length from 9 to 13 kb, are made up of ten functional genes, of which three encode structural proteins (capsid, C; pre-membrane, prM; and envelop, E), and the remaining seven encode non-structural proteins (NS1, NS2a, NS2b, NS3, NS4a, NS4b, and NS5) [2]. According to known host types and phylogenetic relationships, members of the genus Flavivirus are divided into four groups. These include: 1) mosquito-borne flaviviruses—such as Dengue virus (DENV), Japanese encephalitis virus (JEV), Yellow fever virus, West Nile virus (WNV), and Zika virus—the majority of which are zoonoses that infect a range of vertebrate hosts via mosquitoes and can cause public health issues; 2) tick-borne flaviviruses, like the tick-borne encephalitis virus; 3) no-known vector vertebrate flaviviruses, which can replicate in vertebrate cells but with no known arthropod hosts; and 4) insect-specific flaviviruses (ISFVs), which are specific for mosquitoes and unable to replicate in mammalian cells [2,3,4]. The ISFVs are speculated to represent the genetic primordial form of the genus [5].

Based on the molecular structures and evolution rates of known flaviviruses, approximately 2000 unidentified flaviviruses are estimated to exist [6]. Cell fusing agent virus (CFAV), originally discovered and characterised by syncytium formation of cytopathic effect (CPE) in the Aedes aegypti cell line [7], is recognised as the first member of the ISFV group, and was isolated from several species of natural mosquitoes from Puerto Rico in 2006, including Ae. aegypti, Ae. albopictus, and Culex sp. [8]. The second member of this group, Kamiti River virus (KRV), isolated from Ae. macintoshi in Central Province, Kenyan in 2003, is considered to be the first ISFV isolated from nature [9]. Culex flavivirus (CxFV) was first isolated from Cx. pipiens in 2007 in Japan and Indonesia [5]. Moreover, CxFVs have subsequently been detected in Cx. pipiens, Cx. quinquefasciatus, Cx. restuans, Cx. tritaeniorhynchus, and Anopheles sinensis in Japan, the USA, China, Uganda, Mexico, and Brazil [1, 10,11,12,13,14,15,16]. Several novel flaviviruses have been isolated and characterised as ISFVs more recently. They include: Aedes flavivirus (AEFV), which was detected in Ae. flavopictus and Ae. albopictus mosquitoes from Japan [3]; Quang Binh virus (QBV), isolated from Cx. tritaeniorhynchus from Vietnam [17]; Nakiwogo virus (NAKV), isolated from Mansonia africana nigerrima from Uganda [12]; Calbertado virus, detected primarily in Cx. trasalis in Western Canada [18]; Spanish Culex flavivirus, detected in Cx. theileri and Cx. pipiens from Spain [19]; Ochlerotatua flavivirus (OcFV) detected in Ochlerotatua caspius from Italy, Portugal and Spain [19,20,21]; Hanko virus (HANKV), isolated from Oc. punctor and/or Oc. caspius from Finland [22]; and Yamadai flavivirus (YDFV), isolated from Cx. tritaeniorhynchus from Japan [23]. Additionally, flavivirus RNAs have been detected in sandflies [24], and one was also detected in chironomids, another family from the order Diptera [25].

Culex flavivirus was the first ISFV reported in China, where it was isolated in 2012 from Cx. pipiens in Dongming, Shandong Province [13]. Since then, an increasing number of CxFV detections have been documented in China [14, 16, 26]. Other ISFVs that have been isolated and characterised in China are the Quang Binh-like virus [27] and Yunnan Culex flavivirus (YNCxFV) [26]. Both of their nucleotide sequences were 83% identical to that of QBV, a species identity percentage that is just below the cut-off value of 84% proposed by Kuno [28], and were detected mainly in Cx. tritaeniorhynchus, the same host in which QBV has been detected [17, 26, 27]. However, it is unclear whether they are novel species of ISFVs, because, they were collected from the Yunnan Province, which borders Vietnam, the country where QBV was discovered [26]. In order to resolve this issue, neutralisation tests are recommended, when nucleotide sequence identity is higher than 80% [28].

The ISFV CPE can be weak, strain-dependent, or only visible after a number of blind passages [11, 29], indicating that a fine balance exists between virus escape mechanisms and the mosquito immune system. In contrast to most flaviviruses, ISFVs show peculiar characteristics. The sequences of some ISFVs can integrate into the genomes of mosquitoes both in the field and in laboratory cell lines, mainly occurring in Aedes mosquitoes [19, 30, 31]. Thus, when attempting to detect ISFVs, nucleic acid extracts should be treated with DNase to exclude interference by integrated DNA forms of the virus.

Insect-specific flaviviruses have been propagated experimentally only in mosquito cells, and are probably transmitted vertically because eggs, larvae, female, and male mosquitoes can be infected in nature [1, 8, 9, 11, 32, 33] and may be transmitted horizontally through the feeding process in the field [29]. It is apparent that ISFVs are unable to affect the health of birds, domestic animals, and human beings. However, they are carried by invasive mosquitoes and likely do not cross-react with known pathogenic flaviviruses [1]. Recent studies, both in the field and in the laboratory, have shown that co-infection with CxFV increases WNV infection rate [34, 35]. Alternatively, the possibility that CxFV can be attributed to a decrease in WNV morbidity in the Yucatan Peninsula of Mexico, through reducing the number of available competent vectors has been speculated [1]. Whether superinfection exclusion of pathogenic arbovirus occurs in invasive mosquitoes infected with ISFV requires further investigation. Immunological research has shown that, during co-infection by both human pathogens and ISFVs, ISFVs may induce vertebrate cell CPE when innate immunity fails due to the activity of human pathogens [36]. In contrast, it has also been proposed that ISFVs may modulate, or even suppress the immune responses of mosquitoes, thereby making the mosquitoes more susceptible to infection by a broad range of human pathogens [35].

Clearly, research on ISFVs is progressing rapidly and expanding on a worldwide scale. ISFVs in China are far more diverse than is currently reported. However, in most parts of China, ISFVs have not been studied. The purpose of this study is to identify the presence and the infection rate of ISFVs in field-caught mosquitoes in the municipality of Shanghai, China.

Methods

Survey area and mosquito collection

A surveillance program for vector pathogen, involving the collection of mosquitoes and detection of ISFVs studied the distribution and diversity of ISFVs in Shanghai from June to October 2016. Located in the Yangtze River Delta, Shanghai sits on the south edge of the estuary of the Yangtze in the middle portion of the East China coast. The region has a temperate climate, which is suitable for the reproduction of mosquitoes. Shanghai is at risk for transmission of mosquito-borne diseases, especially since it is the centre for economic trade and tourism, and has an abundance of migratory birds. The dominant mosquito species detected in this area are Ae. albopictus, Cx. pipiens, Cx. tritaeniorhynchus, and An. sinensis [37]. In this survey, six different districts, Huangpu District, Songjiang District, Jiading District, Qingpu District, Chongming Districtrict, and Pudong New Area, were chosen for the field research, covering different types of ecological environments, including urban areas, suburban areas, rural areas, and even conservation areas. During the surveillance period, two sampling methods were used to collect mosquitoes. CO2-baited traps were hung from sunset to sunrise in five monitoring sites of each district. Direct aspiration by humans was also used to bait mosquitoes after sunset for fifteen minutes in ten monitoring sites of each district. These two methods were carried out three times at even distributed days throughout the month.

ISFV identification

After collection, mosquitoes were identified using morphological characteristics according to the national key [38]. Some morphological confused specimens were determined by molecular methods, as reported previously [39, 40]. Mosquitoes were then pooled by species, sex, and date and location of collection, with one to 50 individuals per pool. Pooled mosquitoes were stored in 2 ml sterile plastic tubes containing 75% alcohol and frozen at − 20 °C. Sampled pools of mosquitoes were homogenised in a frozen block using a Mixer Mill (Jingxin, Shanghai, China) with one 3 mm and one 5 mm stainless-steel bead added. If fewer than 20 mosquitoes were present in the sample, 450 μl of TRIzol (Invitrogen, Carlsbad, CA) was added, and if more than 20 mosquitoes were present, 600 μl of TRIzol was added to the sample. The samples were then centrifuged at full speed (12 700 rpm) for 10 mins at 4 °C. Then, 350 μl of supernatant from each ground sample was added to the processing cartridge and placed into to the MagNA Pure 96 System (Roche, Basel, Switzerland) for automated RNA extraction. The MagNA Pure 96 Cellular RNA Large Volume Kit, including DNase to obtain DNA-free RNA was used for the RNA extraction according to the manufacturer’s instructions, with a final elution volume of 50 μl. First-strand cDNA was synthesised by RT-PCR using Takara PrimeScript RT reagent Kit with gDNA Eraser (TAKARA BIO, Shiga, Japan). After the RT reactions, PCR amplification targeting mosquito 18S rRNA was implemented with primers 18S417 and 18S920c [41] to verify the integrity of the RNA in each pool. If the control amplification was successful, the cDNA was amplified by hemi-nested PCR using a set of primers (cFD2 and MAMD, cFD2 and FS778) for the detection of the partial NS5 gene of flaviviruses as reported [42]. With the purpose to identify the genotypes of the CxFVs in the samples, we used CxFV-E forward and reverse primers (CxFV-E-F: 5′-ACTGGTGACGTTCAAGGCCATAAG-3′; CxFV-E-R: 5′-GCCGTGATCAGGTGCTGGTCATCG-3′), which target the E gene [43]. Amplified products were separated by agarose gel electrophoresis, then purified and sequenced in both directions by Sangon Biotech (Shanghai, China). Sequences were compared with those available in GenBank database using the BLAST program.

Phylogenetic analysis

Multiple sequence alignments were generated with the relevant ISFV genomes obtained from GenBank and this study using ClustalW2 [44] with default settings, and were manually adjusted if necessary. Neighbour-joining trees were established following the Kimura’s two-parameter (K2P) distance model [45] with 1000 bootstrap replications using MEGA v7.0 software [46]. Based on the Akaike Information Criterion, the best-fit model for the alignment was determined using Modeltest 3.7, in cooperation with PAUP* v4.0b10 [47]. Consequently, the construction of the maximum likelihood (ML) and Bayesian likelihood was completed under the GTR + I + G model for both the NS5 gene and CxFV E gene. The ML tree was performed by MEGA v7.0 software, with 1000 bootstraps. The Bayesian tree was built with MrBayes v3.2.1 [48], run for 10 million generations, with the first 25% generations discarded as burn-in. The trees were unrooted to provide the least biased topology, and visualised using Figtree v1.4.2 (http://tree.bio.ed.ac.uk/software/figtree/).

Calculation of infection rates

The sizes of the pools of collected mosquitoes varied considerably. For this reason, infection rates were calculated by Bias-corrected maximum likelihood estimation (MLE) and minimum infection rate (MIR) using the Excel add-in PooledInfRate v4 statistical software package [49]. The rates are expressed as the number of infected mosquitoes per 1000 collected mosquitoes.

Results

Detection of ISFV RNA in mosquito pools collected in Shanghai

A total of 3249 mosquitoes of the species Ae. albopictus, 3370 of the species Cx. pipiens, and 6045 of the species Cx. tritaeniorhynchus were collected in Shanghai during the active season for mosquitoes between June and October 2016. Of these, 161/440 pools of Ae. albopictus, 228/611 pools of Cx. pipiens, and 195/655 pools of Cx. tritaeniorhynchus were randomly tested for the presence of ISFVs. In order to ensure that pools only tested positive due to the amplification of viral RNA, but not due to the integration of ISFV sequences in the genome of the mosquitoes, RNA extraction protocols employed an additional DNAse step and nucleic acid extracts were treated with DNA Eraser before RT-PCR. The integrity of RNA extracted from the mosquito pools was checked by amplification of the mosquitoes’ 18S rRNA. This amplification was successful in 573 of the 576 pools. Sampled pools were screened for ISFVs using hemi-nested PCR. A fragment comprising 261 bp of the NS5 gene was amplified from pools positive for ISFVs. All the NS5 genes of the 30 ISFVs obtained in this study were detected from female mosquito pools. These sequences have already been deposited in GenBank. Details of their collection information and GenBank accession numbers are listed in Table 1.

Table 1 Summary of the insect-specific flaviviruses detected from different mosquito pools, captured in Shanghai, during June to October in 2016

A total of 44 nucleotide sequences from representative ISFVs that were retrieved from GenBank and the 30 sequences described here were aligned with a partial NS5 nucleotide sequence. This region was selected because the use of the NS5 gene as a target for diagnostic PCR protocols to detect flaviviruses makes partial sequences of this gene the most common. The topology of the NJ tree, ML tree, and Bayesian tree were almost identical for the major lineages, although node confidence values were slightly different among the three (Fig. 1). Hence, only the ML tree is presented with an unrooted phylogram to avoid assumptions regarding ancestry. As shown in Fig. 1, all the lineages, including individuals representing the same ISFV, formed a distinct clade with a high bootstrap value. The tree showed that the ISFVs can be divided into four main clusters, Culex-hosted ISFVs, including CxFV, Cx. theileri flavivirus, YDFV, YNCxFV and QBV; Mansonia-hosted ISFV, only so far containing NAKV; Ochlerotatua-hosted ISFV, comprising with OcFV and HANKV; and Aedes-hosted flaviviruses, including CFAV, AEFV, and KRV. The tree topology indicates that the Shanghai strains belong to the same group of ISFVs, share high sequence identity, and form independent clusters.

Fig. 1
figure 1

Molecular phylogeny of the 30 insect-specific flavivirus (ISFV) strains from Shanghai and other ISFVs based on partial NS5 gene (261 base pairs) sequences. The maximum likelihood tree was constructed by the GTR + I + R model. Sequences referable to the same host-genus are shown with the same colours in the tree. Bootstrap values (1000 replicates, not shown for less than 60%) of maximum likelihood, neighbour-joining, and Bayesian analyses are shown above the main lineages. The bar indicates 0.05 substitutions per site. AEFV, Aedes flavivirus; CFAV, cell fusing agent virus; CTFV, Culex theileri flavivirus; CxFV, Culex flavivirus; HANKV, Hanko virus; KRV, Kamiti River virus; NAKV, Nakiwogo virus; QBV, Quang Binh virus; OcFV, Ochlerotatua flavivirus; SCxFV, Spanish Culex flavivirus; YDFV, Yamadai flavivirus; YNCxFV, Yunnan Culex flavivirus. See Table 1 for information on the sequences of the Shanghai strains of ISFVs

Sequence analysis and phylogenetic characterisation of CxFV genotypes

The CxFVs were detected in nine pools of Cx. pipiens, distributed in Pudong New Area, Chongming District, and Huangpu District. Analysis of their partial NS5 sequence fragments revealed 96.17–97.30% similarity (BLASTn) with the sequences from CxFVs available in GenBank. The CxFV specific primers CxFV-E-F and CxFV-E-R that target a 1443 nt region of the CxFV E gene were used to determine the genotypes of the Shanghai CxFV strains. The amplification of E gene was successful in five of nine CxFV-positive pools. The E sequences of the five pools of mosquitoes showed a high level of identity (99.51–100.00%) with each other, and showed the highest similarity to sequence from the Liaoning DG1064 strain (JQ308188, 99.14–99.51% identity) at the nucleotide level, and lower homology to the Taiwan TW100322 strain (JX897904, 89.35–90.91% identity). The phylogenetic tree presented in Fig. 2 was constructed based on 38 E gene sequences of CxFVs. It was composed of two main monophyletic clades of high confidence. Clade 1 included sequences obtained from Asia and USA, detected primarily in Cx. pipiens. Clade 2 comprised CxFVs from Africa, the Caribbean, and Latin America that share the same host, Cx. quinquefasciatus. These two clades represented two genotypes. All the CxFV strains from Shanghai that were detected in this paper belong to the former Asia/USA clade. Genetic distance analyses based on E gene sequences showed that, at the nucleotide level, the K2P corrected genetic distances were 0.035 within Clade 1, 0.038 within Clade 2, and 0.113 between the two clades.

Fig. 2
figure 2

Molecular phylogenetic analysis of the five Culex flavivirus (CxFV) strains based on the nucleotide sequence of the envelope gene. The maximum likelihood tree was constructed by the GTR + I + R model. Bootstrap values (1000 replicates, not shown for less than 60%) of maximum likelihood, neighbour-joining, and Bayesian analyses are shown above the main lineages. The bar indicates 0.01 substitutions per site. Dark grey and light grey indicate CxFV genotype I and genotype II, respectively

Based on the short NS5 sequence fragments, the genetic identity between YNCxFV (KC464457) and CxFVs was 76.98–80.68%, while that between YNCxFV and QBV (NC_012671) was 85.88%. In the NS5 tree (Fig. 1), the YNCxFV sequence presented as a sister group to a group containing several CxFV sequences and the other composited sequences from QBVs.

Phylogenetic analysis of AEFV sequences

Eleven pools of Ae. albopictus tested positive for AEFVs. It is noteworthy that, of these 11 strains, nine were collected from three different streets in Songjiang District in June, and two were collected separately from Qingpu and Huangpu Districts. The results of the pairwise alignment showed that the new strains shared 98.85–100.00% nucleotide identity and were most similar to the AEFV Narita-21 strain (AB488408) from Japan, which shared 90.80–94.64% nucleotide identity with the partial NS5 gene.

Phylogenetic analysis of QBV sequences

Ten pools tested positive for QBVs. All QBVs described in this paper were detected from pools of Cx. tritaeniorhynchus that were collected in Chongming District of Shanghai. The NS5 gene sequences amplified from QBV-positive pools of Cx. tritaeniorhynchus share 89.77–91.60% identity with the QBV sequence (FJ644291) from Vietnam.

Infection rate of ISFVs in mosquitoes from Shanghai

The overall maximum likelihood estimation (MLE) values (and 95% confidence interval, 95% CI) of CxFV, QBV, and AEFV, expressed as the number of infected mosquitoes per 1000, for Cx. pipiens, Cx. tritaeniorhynchus, and Ae. albopictus were 1.34 (0.66–2.45), 1.65 (0.87–2.85), and 1.51 (0.77–2.70), respectively. These values indicate the differential spread of ISFVs among the different mosquito species. The infection rate of ISFVs also varied based on district. As shown in Table 2, the MLE-MIR of the AEFVs varies across the districts, with a mean value from 2.42 to 33.93 and an upper limit of 62.14 per 1000 Ae. albopictus mosquitoes.

Table 2 Maximal likelihood estimation (MLE) and minimum infection rate (MIR) of insect-specific flaviviruses during mosquito activity season of Shanghai in 2016

Discussion

The most common fragment of flavivirus genome available in the GenBank database is the widely used NS5 gene. Previous results have shown that the NS5 dataset and topologies are not significantly different from those inferred by the open reading frame; though the bootstrap supporting values at some nodes were lower [21, 28]. The phylogenetic tree of the ISFVs (Fig. 1) suggests divergence of four groups, reflecting sequences isolated from Aedes, Culex, Mansonia, and Ochlerotatua mosquitoes. The tree agrees with the hypothesis that each insect flavivirus is maintained in a host genus-associated manner; though the transmission patterns of a few members of this group, like CFAV [8] and CxFV [50], show taxonomic diversity, possibly resulting from frequent host-switching [4].

As evidenced from sequence identity and phylogenetic analysis, one AEFV strains in Ae. albopictus, and two CxFV strains in Cx. pipiens, captured in Huangpu District; nine AEFV strains in Ae. albopictus, captured in Songjiang District; ten QBV strains in Cx. tritaeniorhynus, and one CxFV strain in Cx. pipiens, captured in Chongming District; one AEFV strains in Ae. albopictus, captured in Qingpu District; and six CxFV strains in Cx. pipiens, were captured in Pudong New Area of Shanghai. No ISFV-positive mosquitoes were identified in Jiading District. All Shanghai strains of AEFVs, CxFVs, and QBVs were separately collected from mosquitoes belonging to the species Ae. albopictus, Cx. pipiens, and Cx. tritaeniorhynchus, indicating that ISFVs in Shanghai are probably host-related. Reports of ISFVs in China were organised and listed in Additional file 2. In this study, CxFV and QBV, previously detected in other regions of China [13,14,15,16, 26, 27], are also found in Shanghai. However, AEFV has not been found in China before, but this study is the first to report AEFV in China. At least, three insect flaviviruses including ISFVs, AEFV, CxFV, and QBV are detected in Shanghai, showing that they are co-circulating in this area. The presence of ISFVs in this area could be explain by the fact that Shanghai is an international and cosmopolite metropolis, where migratory vectors through frequent international trades and human migration could easily introduce insect flaviviruses.

The monthly MLE value (95% CI) of AEFV reached 40.18 (20.00–73.57) per 1000 mosquitoes, corresponding to approximately one in three of the homogenates tested. Other researches have also reported high infection rates of mosquitoes by ISFVs in the field [51, 52]. We continued to randomly test collected Ae. albopictus for AEFVs from July to October, but did not detect any further AEFV-positive pools in this district. This finding might indicate that the activity of these viruses in the local Ae. albopictus population is seasonal, an activity pattern that would support the hypothesis that the ability of ISFVs to transmit may vary seasonally, as proposed by Kim et al. [10]. Two other AEFV-positive pools were collected from Qingpu and Huangpu districts in August and September, respectively. The results of a pairwise alignment showed that these new strains were genetically closely related to the AEFV Narita-21 strain (AB488408) from Japan. This observation indicates that AEFV may have been introduced in Shanghai by tire-travelling tiger mosquitoes.

Culex flaviviruses are widely distributed in the nature and can infect various mosquito species [1, 10,11,12,13,14,15,16]. In China, CxFVs have been isolated and/or detected in mosquitoes from Shandong, Liaoning, Gansu, Henan, and Shanxi provinces [15, 16] since they were first reported in 2012 [13]. Unfortunately, the genomic information of the CxFVs isolated from Henan, Shandong, Shanxi, and Gansu provinces is not available in public databases for further phylogenetic analysis. All the CxFVs described here were isolated from Cx. pipiens; though CxFVs have occasionally also been isolated from other Culex species [10,11,12, 53], and have been detected in An. sinensis in Liaoning Province, China [26]. CxFVs can be divided into two genotypes, the Asia/USA genotype and the Africa/Caribbean/Latin American genotype, according to their E genes [50]. The E genes of the CxFVs in five of the nine CxFV positive pools collected in this study were successfully amplified. Phylogenetic analyses indicated that the CxFVs of Shanghai strains are more closely related to the subtropical lineage than to the tropical one. Based on cytology experiments, the Asia/USA genotype has been reported to cause the cell-fusion type of CPE in C6/36 cells [9, 13, 54]. The Africa/Caribbean/Latin American genotype, in contrast, induces relatively mild CPE, reduced cell density, and modification of cell shape [12], or no observed CPE [1], except for the isolate from Taiwan Province, China, which possesses a glycine residue at position 117 of its E protein. The presence of the glycine residue is a unique character shared by the Asia/USA genotype, which may be associated with the induction of CPE [11]. In this study, all the CxFVs of Shanghai strains had a unique glycine residue at position 117 of their E proteins, which is consistent with the Asia/US genotype.

The Quang Binh virus has been found in southwest Asia in Quang Binh, Vietnam [17] and in Yunnan Province in China [27]. The QBVs detected in this study were solely isolated in Cx. tritaeniorhynchus from the Chongming District. Chongming Island has important commercial harbours and a nature reserve in the Dongtan wetland that hosts numerous migratory birds. Migratory vectors through international trade, or wind-blown infected mosquitoes might introduce QBV, as other ISFVs [52].

It is unfortunate that virus isolation could not be done because the mosquitoes were preserved in 75% alcohol. The potential virus in the supernatant of mosquito homogenate cannot be cultivated in mosquito cell lines for further analyses. Thus, the actual infection rate of ISFVs in Shanghai may be underestimated in this study. This study might leverage deep genomic investigations and continuous ISFV mosquito surveillance for a better understanding of the ISFVs circulating in mosquitoes found in Shanghai. In the further, studies might be carried out to investigate any co-effect of ISFV and JEV, or DENV in field caught mosquitoes, as well as the impact on other arboviruses of public health importance such as JEV, DENV, and Zika virus.

Conclusions

This molecular survey has determined the presence, geographic distribution, genetic variation and infection rate of ISFVs in Shanghai, China. The results showed that, AEFV, CxFV, and QBV are co-circulating in Shanghai, with overall MLE values of 1.51, 1.34, and 1.65 per 1000 of Ae. albopictus, Cx. pipiens, and Cx. tritaeniorhynchus, respectively. Phylogenetic analysis of E gene sequences with those of reference strains has revealed that the Shanghai CxFVs belong to the Asia/USA genotype. Very importantly, this is the first report of AEFV in China.

Abbreviations

AEFV:

Aedes flavivirus

bp:

Base pair

CFAV:

Cell fusing agent virus

CI :

Confidence interval

CPE:

Cytopathic effect

CxFV:

Culex flavivirus

DENV:

Dengue virus

E gene:

Envelope gene

HANKV:

Hanko virus

ISFV:

Insect-specific flavivirus

JEV:

Japanese encephalitis virus

K2P:

Kimura’s 2-parameter

KRV:

Kamiti River virus

MIR:

Minimum infection rate

ML tree:

Maximum likelihood tree

MLE:

Maximal likelihood estimation

NAKV:

Nakiwogo virus

NS5:

Nonstructural 5 gene

OcFV:

Ochlerotatua flavivirus

QBV:

Quang Binh virus

WNV:

West Nile virus

YDFV:

Yamadai flavivirus

YNCxFV:

Yunnan Culex flavivirus

References

  1. Farfan-Ale JA, Loroño-Pino MA, Garcia-Rejon JE, Hovav E, Powers AM, Lin M, et al. Detection of RNA from a novel West Nile-like virus and high prevalence of an insect-specific flavivirus in mosquitoes in the Yucatan peninsula of Mexico. Am J Trop Med Hyg. 2009;80:85–95.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  2. Lobo FP, Mota BE, Pena SD, Azevedo V, Macedo AM, Tauch A, et al. Virus-host coevolution: common patterns of nucleotide motif usage in Flaviviridae and their hosts. PLoS One. 2009;4:e6282.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. 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–29.

    Article  PubMed  CAS  Google Scholar 

  4. Cook S, Moureau G, Kitchen A, Gould EA, de Lamballerie X, Holmes EC, et al. Molecular evolution of the insect-specific flaviviruses. J Gen Virol. 2012;93:223–34.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Hoshino K, Isawa H, Tsuda Y, Yano K, Sasaki T, Yuda M, et al. Genetic characterization of a new insect flavivirus isolated from Culex pipiens mosquito in Japan. Virology. 2007;359:405–14.

    Article  PubMed  CAS  Google Scholar 

  6. Pybus OG, Rambaut A, Holmes EC, Harvey PH. New inferences from tree shape: numbers of missing taxa and population growth rates. Syst Biol. 2002;51:881–8.

    Article  PubMed  CAS  Google Scholar 

  7. Stollar V, Thomas VL. An agent in the Aedes aegypti cell line (Peleg) which causes fusion of Aedes albopictus cells. Virology. 1975;64:367–77.

    Article  PubMed  CAS  Google Scholar 

  8. Cook S, Bennett S, Holmes E, de Chesse R, Moureau G, de Lamballerie X. Isolation of a new strain of the flavivirus cell fusing agent virus in a natural mosquito population from Puerto Rico. J Gen Virol. 2006;87:735–48.

    Article  PubMed  CAS  Google Scholar 

  9. Crabtree MB, Sang RC, Stollar V, Dunster LM, Miller BR. Genetic and phenotypic characterization of the newly described insect flavivirus, Kamiti River virus. Arch Virol. 2003;148:1095–118.

    Article  PubMed  CAS  Google Scholar 

  10. Kim DY, Guzman H Jr, Bueno R, Dennett JA, Auguste AJ, Carrington CVF, et al. Characterization of Culex flavivirus (Flaviviridae) strains isolated from mosquitoes in the United States and Trinidad. Virology. 2009;386:154–9.

    Article  PubMed  CAS  Google Scholar 

  11. Chen Y, Lin J, Fan Y, Tu W, Chang G, Chiou S. First detection of the Africa/Caribbean/Latin American subtype of Culex flavivirus in Asian country, Taiwan. Comp Immunol Microbiol Infect Dis. 2013;36:387–96.

    Article  PubMed  Google Scholar 

  12. Cook S, Moureau G, Harbach RE, Mukwaya L, Goodger K, Ssenfuka F, et al. Isolation of a novel species of flavivirus and a new strain of Culex flavivirus (Flaviviridae) from a natural mosquito population in Uganda. J Gen Virol. 2009;90:2669–78.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  13. Wang H, Wang H, Fu S, Liu G, Liu H, Gao X, et al. Isolation and identification of a distinct strain of Culex flavivirus from mosquitoes collected in mainland China. Virol J. 2012;9:73.

    Article  CAS  Google Scholar 

  14. Liang W, He X, Liu G, Zhang S, Fu S, Wang M, et al. Distribution and phylogenetic analysis of Culex flavivirus in mosquitoes in China. Arch Virol. 2015;160:2259–68.

    Article  PubMed  CAS  Google Scholar 

  15. Zha B, Yu D, Fu S, Li G, Lv Z, Jiang J, et al. Investigation of mosquitoes and arboviruses in Hexi corridor of Gansu province, China in 2011. Chin J Vector Biol Control. 2012;23:424–7. (in Chinese)

    CAS  Google Scholar 

  16. An S, Liu J, Ren Y, Wang Z, Han Y, Ding J, et al. Isolation of the Culex flavivirus from mosquitoes in Liaoning Province, China. Chin J Virol. 2012;28:511–6. (in Chinese)

    CAS  Google Scholar 

  17. Crabtree MB, Nga PT, Miller BR. Isolation and characterization of a new mosquito flavivirus, Quang Binh virus, from Vietnam. Arch Virol. 2009;154:857–60.

    Article  PubMed  CAS  Google Scholar 

  18. Tyler S, Bolling BG, Blair CD, Brault AC, Pabbaraju K, Armijos MV, et al. Distribution and phylogenetic comparisons of a novel mosquito flavivirus sequence present in Culex tarsalis mosquitoes from western Canada with viruses isolated in California and Colorado. Am J Trop Med Hyg. 2011;85:162–8.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Vazquez A, Sanchez-Seco MP, Palacios G, Molero F, Reyes N, Ruiz S, et al. Novel flaviviruses detected in different species of mosquitoes in Spain. Vector Borne Zoonotic Dis. 2012;12:223–9.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Ferreira DD, Cook S, Lopes A, de Matos AP, Esteves A, Abecasis A, et al. Characterization of an insect-specific flavivirus (OCFVPT) co-isolated from Ochlerotatua caspius collected in southern Portugal along with a putative new Negev-like virus. Virus Genes. 2013;47:532–45.

    Article  PubMed  CAS  Google Scholar 

  21. Calzolari M, Bonilauri P, Bellini R, Caimi M, Defilippo F, Maioli G, et al. Arboviral survey of mosquitoes in two northern Italian regions in. 2007 and 2008. Vector Bone Zoonotic Dis. 2010;10:875–84.

    Article  Google Scholar 

  22. Huhtamo E, Moureau G, Cook S, Julkunen O, Putkuri N, Kurkela S, et al. Novel insect-specific flavivirus isolated from northern Europe. Virology. 2012;433:471–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Kuwata R, Sugiyama H, Yonemitsu K, Nguyen D, Terada Y, Taniguchi M, et al. Isolation of Japanese encephalitis virus and a novel insect-specific flavivirus from mosquitoes collected in a cowshed in Japan. Arch Virol. 2015;160:2151–9.

    Article  PubMed  CAS  Google Scholar 

  24. Moureau G, Ninove L, Izri A, Cook S, de Lamballerie X, Charrel RN. Flavivirus RNA in phlebotomine sandflies. Vector Borne Zoonotic Dis. 2010;10:195–7.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Cook S, Chung BYW, Bass D, Moureau G, Tang S, McAlister E, et al. Novel virus discovery and genome reconstruction from field RNA samples reveals highly divergent viruses in dipteran hosts. PLoS One. 2013;8:e80720.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Zuo S, Zhao Q, Guo X, Zhou H, Cao W, Zhang J. Detection of Quang Binh virus from mosquitoes in China. Virus Res. 2014;180:31–8.

    Article  PubMed  CAS  Google Scholar 

  27. Feng Y, Li H, Zhu J, Zhang Y, Yang W, Fan J, et al. First isolation of Quang Binh-like virus from mosquitoes in China. Chin J Virol. 2014;30:57–61. (in Chinese)

    CAS  Google Scholar 

  28. Kuno G, Chang GJ, Tsuchiya KR, Karabatsos N, Cropp CB. Phylogeny of the genus Flavivirus. J Virol. 1998;72:73–83.

    PubMed  PubMed Central  CAS  Google Scholar 

  29. Calzolari M, Ze-Ze L, Vazquez A, Sanchez Seco MP, Amaro F, Dottori M. Insect-specific flaviviruses, a worldwide widespread group of viruses only detected in insects. Infect Genet Evol. 2016;40:381–8.

    Article  PubMed  Google Scholar 

  30. Roiz D, Vazquez A, Seco MP, Tenorio A, Rizzoli A. Detection of novel insect flavivirus sequences integrated in Aedes albopictus (Diptera: Culicidae) in northern Italy. Virol J. 2009;6:93.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Sanchez-Seco M, Vazquez A, Collao X, Hernandez L, Aranda C, Ruiz S. Surveillance of arboviruses in Spanish wetlands: detection of new flavi- and phleboviruses. Vector Borne Zoonotic Dis. 2009;10:203–6.

    Article  Google Scholar 

  32. Saiyasombat R, Bolling BG, Brault AC, Bartholomay LC, Blitvich BJ. Evidence of efficient transovarial transmission of Culex flavivirus by Culex pipiens (Diptera: Culicidae). J Med Entomol. 2011;48:1031–8.

    Article  PubMed  Google Scholar 

  33. Bolling B, Eisen L, Moore C, D Blair C. Insect-specific flaviviruses from Culex mosquitoes in Clorado, with evidence of vertical transmission. Am J Trop Med Hyg. 2011;85:169–77.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Kent RJ, Crabtree MB, Miller BR. Transmission of West Nile virus by Culex quinquefasciatus say infected with Culex flavivirus Izabal. PLoS Negl Trop Dis. 2010;4:e671.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  35. Newman CM, Cerutti F, Anderson TK, Hamer GL, Walker ED, Kitron UD, et al. Culex flavivirus and West Nile virus mosquito coinfection and positive ecological association in Chicago, United States. Vector Borne Zoonotic Dis. 2011;11:1099–105.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Tree MO, McKellar DR, Kieft KJ, Watson AM, Ryman KD, Conway MJ. Insect-specific flavivirus infection is restricted by innate immunity in the vertebrate host. Virology. 2016;497:81–91.

    Article  PubMed  CAS  Google Scholar 

  37. Zhou Z, Lv S, Zhang Y, Gu W, Guo Y, Jiang M, et al. Mosquito species, distribution and their pathogens in Shanghai, China. Chin J Vector Biol Control. 2015;26:28–32. (in Chinese)

    Google Scholar 

  38. Lu B, Xu J, Dong X. Fauna Sinica, Insecta, vol.9. Diptera: Culicidae II. 1st ed. Beijing: Science Press; 1997. (in Chinese)

    Google Scholar 

  39. Fang Y, Shi WQ, Zhang Y. Molecular phylogeny of Anopheles hyrcanus group (Diptera: Culicidae) based on mtDNA COI. Infect Dis Poverty. 2017;6:61.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Fang Y, Shi WQ, Zhang Y. Molecular phylogeny of Anopheles hyrcanus group members based on ITS2 rDNA. Parasite Vector. 2017;10:417.

    Article  Google Scholar 

  41. Hoffmann PR, Woodrow RJ, Calimlim PS, Sciulli RH, Effler PV, Miyamoto V, et al. West Nile virus surveillance: a simple method for verifying the integrity of RNA in mosquito (Diptera: Culicidae) pools. J Med Entomol. 2004;41:731–5.

    Article  PubMed  CAS  Google Scholar 

  42. Scaramozzino N, Crance JM, Jouan A, Debriel DA, Stoll F, Garin D. Comparison of flavivirus universal primer pairs and development of a rapid, highly sensitive heminested reverse transcription-PCR assay for detection of flaviviruses targeted to a conserved region of the NS5 gene sequences. J Clin Microbiol. 2001;39:1922–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Saiyasombat R, Dorman KS, Garcia-Rejon JE, Lorono-Pino MA, Farfan-Ale JA, Blitvich BJ. Isolation and sequence analysis of Culex flavivirus from Culex interrogator and Culex quinquefasciatus in the Yucatan peninsula of Mexico. Comp Immunol Microbiol Infect Dis. 2010;155:983–6.

    CAS  Google Scholar 

  44. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.

    Article  PubMed  CAS  Google Scholar 

  45. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    Article  PubMed  CAS  Google Scholar 

  46. Wilgenbusch JC, Swofford D. Inferring evolutionary trees with PAUP*. Curr Protoc Bioinformatics. 2003;6:6.4.

    Google Scholar 

  47. Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61(3):539–42.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Morariu VI, Srinivasan BV, Raykar VC, Duraiswami R, Davis LS. Automatic online tuning for fast Gaussian summation. Adv Neural Inf Proces Syst. 2008;21:1113–20.

    Google Scholar 

  49. Biggerstaff BJ. PooledInfRate, version 3.0: a Microsoft® excel® add-in to compute prevalence estimates from pooled samples. Fort Collins: CDC; 2006.

    Google Scholar 

  50. Bittar C, Machado DC, Vedovello D, Ullmann LS, Rahal P, Junior JPA, et al. Genome sequencing and genetic characterization of Culex flavirirus (CxFV) provides new information about its genotypes. Virol J. 2016;13:158.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  51. Farfan-Ale JA, Loroño-Pino MA, Garcia-Rejon JE, Soto V, Lin M, Staley M, et al. Detection of flaviviruses and orthobunyaviruses in mosquitoes in the Yucatan peninsula of Mexico in 2008. Vector Borne Zoonotic Dis. 2010;10:777–83.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Roiz D, Vazquez A, Rosso F, Arnoldi D, Girardi M, Cuevas L, et al. Detection of a new insect flavivirus and isolation of Aedes flavivirus in northern Italy. Parasite Vector. 2012;5:223.

    Article  Google Scholar 

  53. Obaranagoya M, Yamauchi T, Watanabe M, Hasegawa S, Iwaiitamochi M, Horimoto E, et al. Ecological and genetic analyses of the complete genomes of Culex flavivirus strains isolated from Culex tritaeniorhynchus and Culex pipiens (Diptera: Culicidae) group mosquitoes. J Med Entomol. 2013;50:300–9.

    Article  CAS  Google Scholar 

  54. Blitvich BJ, Lin M, Dorman KS, Soto V, Hovav E, Bradley JT, et al. Genomic sequence and phylogenetic analysis of Culex flavivirus, an insect-specific flavivirus, isolated from Culex pipiens (Diptera: Culicidae) in Iowa. J Med Entomol. 2009;46:934–41.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Acknowledgements

We are grateful to the staffs of the Centers for Disease Control and Prevention of Songjiang District, Chongming District, Qingpu District, Huangpu District, Jiading District, and Pudong New Area, Shanghai, China for their assistance in carrying out field work. We thank Dr. Kokouvi Kassegene for providing an excellent assistance with English editing during preparation of this manuscript.

Funding

This study was supported by the National Key Research and Development Program of China (No. 2016YFC1202000, 2016YFC1202002, 2016YFC1202003), and the Fourth Shanghai Public Health Three-Year Action Plan (No. GWIV-29).

Availability of data and materials

All data generated or analysed during this study are included in the published article and Additional file.

Author information

Authors and Affiliations

Authors

Contributions

YF, YZ, ZZ, WS, SX, YL, JW and QL designed the study; YF, ZZ, JW and WS conducted the field collections; YF, ZZ, YL and GL carried out laboratory experiments, performed phylogenetic analyses; YF wrote the original manuscript; SX, QL, GL and YZ revised the manuscript and gave approval for its publication. All the authors read and approved the final version of the manuscript.

Corresponding author

Correspondence to Yi Zhang.

Ethics declarations

Ethics approval and consent to participate

No specific permits were required for this study. The study did not involve endangered or protected species. Therefore, the local ethics committee deemed that approval was not necessary.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional files

Additional file 1:

Multilingual abstracts in the five official working languages of the United Nations. (PDF 271 kb)

Additional file 2:

Summary of the ISFVs documented in mainland China. (PDF 981 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fang, Y., Zhang, Y., Zhou, ZB. et al. Co-circulation of Aedes flavivirus, Culex flavivirus, and Quang Binh virus in Shanghai, China. Infect Dis Poverty 7, 75 (2018). https://doi.org/10.1186/s40249-018-0457-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40249-018-0457-9

Keywords