Isolation and identification of a distinct strain of Culex Flavivirus from mosquitoes collected in Mainland China
Virology Journal volume 9, Article number: 73 (2012)
Culexflavivirus (CxFV) is an insect specific virus that has been isolated from Culexpipiens, Culexquinquefasciatus, Culextritaeniorhynchus and other Culex mosquitoes. It is a novel flavivirus isolated in Asia, North America, Central America and Africa. Phylogenetic analysis indicates that, based on the envelope gene (E gene) sequence, the worldwide CxFV strains can be divided into two genotypes.
A virus (SDDM06-11) was isolated from Culexpipiens collected in Shandong Province, China in 2006. The strain caused cytopathic effect (CPE) in Aedesalbopictus (C6/36) cells by 3 days post-infection and immunofluorescence assay (IFA) showed a reaction with Japanese encephalitis virus (JEV) polyclonal antibodies. Phylogenetic analysis of the E gene sequence showed CxFV formed two genotypes with the SDDM06-11 strain assigned to genotype 1. Analysis of the E gene nucleotide homology showed the virus possessed characteristic amino acids at specific sites. The nucleotide homology of the open reading frame (ORF) was 94.8%-95.1% between SDDM06-11 and isolates from Japan, Iowa and Texas, and 90.2%-90.5% between SDDM06-11 and isolates from Uganda and Mexico.
In this paper we report the first isolation and identification of an isolate of CxFV in mainland China. Phylogenetic analysis indicates the isolate belongs to genotype 1. Our findings provide insight into the occurrence of CxFV in Culex mosquito populations and its distribution on a global scale.
The genus Flavivirus in the family Flaviviridae is a group of single stranded positive sense RNA viruses that are primarily arthropod borne . Phylogenetic analysis classifies these viruses into mosquito-borne, tick-borne, no known vector and viruses that are specific to insects, but which do not infect vertebrates [2, 3]. This final group includes cell fusing agent virus (CFAV) , Kamiti River virus (KRV) and Culexflavivirus (CxFV).
CxFV is an insect-specific virus that can replicate in C6/36 cells but not in other mammalian cell lines. CxFV was first isolated in 2003 in Japan (in species Culexpipiens and Culextritaeniorhynchus) and then in Indonesia in 2004 (Culexquinquefasciatus). Subsequently the virus has been isolated in North America: Iowa and Texas in the United States [7, 8] (2007-2008, Culexquinquefasciatus, Culexrestuans); Central America: Guatemala (2006, Culexquinquefasciatus) , the Yucatan peninsula in Mexico [10, 11] (2007-2008, Culexquinquefasciatus, Culexinterrogator); Trinidad in the West Indies  (2008, Culexquinquefasciatus); and Uganda, Africa  (2008, Culexquinquefasciatus). E gene phylogenetic analysis of these isolates showed that they can be clearly divided into two genotypes [8, 9].
Relatively few flaviviruses have been isolated in China, the most common of which are Japanese encephalitis virus (JEV), dengue virus (DENV) and tick-borne encephalitis virus (TBEV) with presence of associated disease symptoms [13, 14]. In recent years, a novel insect-specific Flavivirus was isolated from Aedesvexans in rural corrals in Chaoyang city in Liaoning province in China , but the distribution and significance of these insect specific viruses remains unclear in China.
In this study, we report the isolation and identification of a distinct strain of CxFV from Culexpipiens mosquitoes collected in Shandong Province in China in 2006 and perform nucleotide sequencing analysis in order to compare with other CxFV isolates previously collected from around the world.
Mosquitoes used in this study were collected in August in 2006 in Dongming county, Shandong province, China. Mosquitoes were collected using light traps (Hubei Lucky Star Environment Protection Co., Ltd.) in fields, pig farms and dwellings. The traps were placed before sunset and mosquitoes were collected from the traps the following morning after sunrise. Collected mosquitoes were frozen for 30 min at -20°C then placed on an ice plate to determine mosquito species and exclude blood-fed and/or male mosquitoes . The mosquitoes were sorted according to species with 50 to 100 adults per pool, transported on dry ice and then stored in liquid nitrogen in the laboratory.
Cell culture and virus isolation
Mosquito cell line C6/36 and baby hamster kidney (BHK-21) cell line were used for virus isolation. The two cell lines were cultured with Eagle's minimum essential medium (MEM) (Sigma-Aldrich, USA) with 10% heat-inactivated fetal bovine serum (FBS) (Sigma), 2% non-essential amino acid (Sigma), 100 U/ml penicillin (Gibco, Invitrogen, USA), 100 μg/ml streptomycin (Gibco), and maintained at 28°C (for mosquito cells) and 37°C (for mammalian cells) with 5% CO2.
Pools of mosquitoes were homogenized using a Mixer Mill (Tissuelyser, Qiagen, Germany) and 3 mm stainless steel beads in 2 mL sterile plastic tubes and 1.5 mL of MEM containing 2% fetal bovine serum, 2 mM L-glutamine, 50 units/mL Penicillin and 100 ug/mL Streptomycin (Gibco Cat No 15140). Samples were centrifuged at 12,000 rpm for 20 minutes at 4°C. Then, 200 μl of the clarified homogenates were inoculated into single 5.5 cm2 Nunc tubes (Nunclon, Denmark) spread with a monolayer of C6/36 and BHK-21 cells respectively for 1 h at constant temperature. After discarding and refreshing with 500 μl medium, the cell cultures were incubated with the same conditions for 6-7 days. Cytopathic effects (CPE) were checked every 8 hours after incubation for 24 hours and observation over the next 6-7 days. At 70% CPE (3 days after post-infection), the culture supernatants were harvested and cellular debris was removed by centrifugation at 12,000 rpm. The supernatants were stored at -80°C until identification.
Indirect immunofluorescence assay
Virus isolate in cell culture was identified by indirect immunofluorescence assay (IFA). Infected and uninfected cell suspensions were applied to Teflon-coated 10-well slides, which were then air dried and fixed in acetone. The IFA was conducted with a panel of antiserum with fluorescein isothiocyanate-conjugated goat anti-mouse IgG as the second antibody (Sigma) by previously described procedures . The panel of polyclonal antibody (JEV-GSS) against JEV was prepared by our laboratory .
RT-PCR and nucleotide sequencing
Viral RNA was extracted from 140 μl of virus culture stocks using the QIAamp viral RNA extraction kit (Qiagen, Valencia, CA, USA) in accordance with the manufacturer's protocol and stored at -70°C until use. Briefly, purified RNA was used as template to finish the reverse transcription-polymerase chain reactions (RT-PCR) for cDNA synthesis using Ready-To-Go™ You-Prime First-Strand Beads (Amersham Biosciences, Piscatway, NJ, USA) and supplied random hexanucleotide primers. Flavivirus genus-specific primer sets  were used for identification. Then ten overlapping primers (Table 1) were designed from the Tokyo strain genome sequence (GenBank number: AB262759) to amplify the complete open reading frame (ORF) nucleotide sequence of the viral genomic RNA.
To avoid any potential contamination, JEV RNA and water were used as positive and negative controls respectively at all stages. RT-PCR products were purified using a commercial kit (ExoSAP-IT PCR Purification Kit, USB, USA). PCR products were sequenced using a PRIMS Ready Reaction Dyedeoxy Terminator Cycle Sequencing Kit on an ABI Prism 3100 Avant Genetic Analyser (Applied Biosystems, California, USA) using the PCR primer sets listed in Table 1. Dye terminator sequencing was used for both strands with double coverage and the results were inspected for accuracy and quality. After sequencing we got the ATGC software package (GENETYX Corp., Tokyo, Japan) was used to complete the sequence assembly.
Sequence data for the E and ORF gene of strain SDDM06-11 were deposited in GenBank with accession number JF938690 and JQ518484, respectively.
Sequence analysis and phylogenetic comparisons
Additional CxFV sequences were downloaded from GenBank, a list of the 49 CxFV isolates used in the analyses with origin and year of isolation is shown in Table 2. Nucleotide and amino acid sequence alignments were generated by CLUSTAL_X version 1.8 . The optimal nucleotide substitution model was estimated as the General Time Reversible with rate heterogeneity among sites and invariable sites using Akaike Information Criterion as implemented in jModeltest v.0.1.1 (http://darwin.uvigo.es/software/jmodeltest.html). The phylogenetic tree was computed with the GTR + I + G model using the maximum likelihood approach with the PHYML software package . The statistical significance of the phylogeny was estimated by the nonparametric bootstrap analysis with 100 replicates.
Neighbor Joining (NJ) and Maximum Likelihood (ML) trees were also estimated for the nucleotide alignments using PHYLIP, v3.6-alpha (http://atgc.lirmm.fr/phyml) and MEGA v5.1 with empirical base frequencies, and a gamma distribution to estimate rate variation among sites. Bootstrap values were determined for 100 replicates.
The trees were rooted using a CRFV (GenBank accession no. NC001564) sequence as the outgroup. Analysis of nucleotide and translated amino acid sequence identities for the E and ORF sequences was performed using the GeneDoc and Lasergene software packages (DNASTAR Inc, USA).
Virus isolation and identification
A total of 80 pools consisting of 4118 Culextritaeniorhynchus, 839 Culexpipiens and 27 Anopheles sinensis were collected from Dongming county, Shandong province, China, from August 9 to 11, 2006 and were processed for virus isolation. A virus isolate was identified that caused CPE in C6/36 cells and which was characterized by marked syncytia, aggregation and fusion by day 3 post-infection (Figure 1A and 1B); no CPE was observed in BHK-21 cells. Antigen IFA using JEV polyclonal antibody (JEV-GSS) identified the isolate as a flavivirus (Figure 2A and 2B). RT-PCR amplification using flavivirus primers FU1 and cFD2 . BLASTN analysis against the nr database at NCBI showed the highest scoring hits were all from CxFV with 95% similarity and an E-value of 0. The isolate was designated SDDM06-11.
Analysis of the E gene sequence
Table 3 shows a summary of the nucleotide and amino acid identities for the E gene according to genotype and geographical location amongst strains used in the phylogenetic analysis.
Specifically, the nucleotide identity between SDDM06-11 and genotype 1 strains (Japan, Indonesia and North America) (Iowa and Texas) was 94.4%-94.8%, and the corresponding identity of SDDM06-11 with genotype 2 (Uganda, Trinidad, Guatemala and Mexico) was 89.4%-90.2%. The identity between genotype 1 and genotype 2 was 89.1%-90.9%. The within genotype identity for genotype 1 was 94.4%-99.6% and within genotype 2 was 97.9%-99.2%. The E gene amino acid identity amongst all the CxFV strains was 97.2%-100%(Table 3).
Examination of the amino acid sequence for the E gene indicated that each genotype contained site-specific mutations (Table 4). A total of ten mutations were identified that were specific to genotype I and 2; these occurred at E-60 (Val), E-64 (Phe), E-117 (Gly), E-185 (His), E-236 (Val), E-326 (Arg), E-340 (Ile), E381 (His) E-407E (Leu), E-423 (Ile). Mutations that were unique to SDDM06-11 (i.e., that were not present in genotype 1 or 2) occurred at E-125 (Val) and E-414 (Val).
Analysis of the ORF gene sequence
The nucleotide and deduced amino acid sequences of the ORF of SDDM06-11 were compared with those of six other CxFV strains selected from the GenBank database. The nucleotide identity between SDDM06-11 and other genotype 1 sequences (Japan, Iowa and Texas) was 94.8%-95.1%, and between SDDM06-11 and genotype 2 sequences (Uganda and Mexico) was 90.2%-90.5%. The calculated amino acid identity for all the CxFV strains was 96.5%-98.2% (Table 5). The highest nucleotide identity occurred in the NS2a and NS2b genes, where the nucleotide identity between SDDM06-11 and genotype 1 strains was > 98% and > 96% respectively;. For genotype 2, the identity was > 94% and > 92% respectively. however, the amino acid homology was lower than ORF (Table 5).
The estimated ML phylogenetic tree based on the E gene nucleotide sequence is shown in Figure 3 and is consistent with the trees estimated for the using the NJ and ML methods using MEGA v5.1 (data not shown). The tree shows the CxFV strains are divided into two distinct genotypes. Genotype 1 is divided into two clades, with clade 1 incorporating sequences from Japan and Indonesia and clade 2 containing samples from the United States (Texas and Iowa) . Genotype 2 contains samples from Africa (Uganda), North America (Mexico), Central America (Guatemala) and the West Indies, (Trinidad) . The newly isolated Chinese strain SDDM06-11 is placed within the genotype 1 clade; however, it appears to have branched from the remaining isolates at a relatively early time point, suggesting that it belongs to a distinct cluster (Figure 3).
We report here the first isolation of a strain of CxFV in China. This strain is distinct from other isolates collected from around the world, containing signature amino acid mutations and forming a separate cluster in the estimated phylogenetic tree. This sample appears therefore to represent a new cluster within genotype 1; this is also supported by the analysis of the amino acid substitutions in the E protein; many of the mutations were shared by these genotypes, whereas SDDM06-11 contains several unique substitutions.
Investigation of the CPE of this strain also produced distinctive results. Previous studies of CPE produced by CxFV in C6/36 cells found that some genotype 1 strains could generate marked CPE in cells within 6 to 7 days but no such changes have been observed for genotype 2 strains [7, 9, 10]. For the SDDM06-11 from China the CPE was similar to observations for genotype 1 [6, 8], producing marked CPE but with a more rapid onset of effects (3-4 days) This difference may be a consequence of the detected amino acid changes in the E gene, but detailed experimental studies are necessary to gain further insight and understanding of the observed differences between the two genotypes.
The genotypes are formed with high bootstrap support and contain clear geographical subdivision, but on a global scale the geographical distribution is less clear with North America (Texas and Iowa), Japan and the Chinese isolate placed in genotype 1 and the Central America (Guatemala and Mexico), Trinidad and Africa strains placed in genotype 2. Similarly, the evolutionary relationship among the three genotypes is not apparent from the estimated phylogenetic tree. This is primarily because there are relatively few samples available for the analysis and only a partial evolutionary history of the virus can be recreated.
All CxFV samples to date have been isolated from members of the Culex genus, with the majority of samples collected from Culexpipiens and Culexquinquefasciatus. Genotype 2 primarily consists of Culexquinquefasciatus isolates, whereas genotype 1 is primarily Culexpipiens and contains smaller numbers of other Culex species [6–11]. The SDDM06-11 belonging to genotype 1 is also isolated from Culexpipiens.
An earlier theoretical study of species diversification derived an estimate for the relationship between the distribution of internal nodes in a tree predicted from a limited sample set and the number of different viruses from which the sample set were selected . When applied to the flaviviruses genus, it was estimated that are several thousand unsampled mosquito borne flaviviruses . Thus further collection and evaluation of CxFV samples is necessary to gain further insight into the evolutionary history of the virus and the host range. It is likely that this will reveal greater complexity and structure (i.e. additional genotypes and new strains) in the CxFV phylogeny and host population.
We have isolated and identified the CxFV strain in mainland China. Phylogenetic analysis places the isolate within genotype 1, but indicates it is quite distinct from other strains within this clade. This is further supported by analysis of the E gene and ORF sequence which shows the isolate has discrete difference in terms of nucleotide composition and amino acid mutations. Our findings suggest that CxFV may be widespread in Culex mosquito populations and could have a relevant role in the evolution of flaviviruses.
Dr. Wang Huanyu, Ph.D., is an associate professor at the State Key Laboratory for Infectious Disease Prevention and Control, the Institute for Viral Disease Control and Prevention, Chinese Center for Disease Control and Prevention. His current research focuses on arbovirus and viral encephalitis, especially in Flavivirus.
Fields BN, Knipe DM, Howley PM (Eds): Fields Virology 5th edition. Lippincott Williams & Wilkins; 2007:1154-1228.
Gaunt MW, Sall AA, de Lamballerie X, Falconar AK, Dzhivanian TI, Gould EA: Phylogenetic relationships of flaviviruses correlate with their epidemiology, disease association and biogeography. J Gen Virol 2001,82(Pt 8):1867-1876.
Kuno G, Chang GJ, Tsuchiya KR, Karabatsos N, Cropp CB: Phylogeny of the genus Flavivirus. J Virol 1998,72(1):73-83.
Stollar V, Thomas VL: An agent in the Aedesaegypti cell line (Peleg) which causes fusion of Aedesalbopictus cells. Virology 1975,64(2):367-377. 10.1016/0042-6822(75)90113-0
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(6):1095-1118. 10.1007/s00705-003-0019-7
Hoshino K, Isawa H, Tsuda Y, Yano K, Sasaki T, Yuda M, Takasaki T, Kobayashi M, Sawabe K: Genetic characterization of a new insect flavivirus isolated from Culexpipiens mosquito in Japan. Virology 2007,359(2):405-414. 10.1016/j.virol.2006.09.039
Blitvich BJ, Lin M, Dorman KS, Soto V, Hovav E, Tucker BJ, Staley M, Platt KB, Bartholomay LC: Genomic sequence and phylogenetic analysis of Culexflavivirus, an insect-specific flavivirus, isolated from Culexpipiens (Diptera: Culicidae) in Iowa. J Med Entomol 2009,46(4):934-941. 10.1603/033.046.0428
Kim DY, Guzman H, Bueno R Jr, Dennett JA, Auguste AJ, Carrington CV, Popov VL, Weaver SC, Beasley DW, Tesh RB: Characterization of CulexFlavivirus (Flaviviridae) strains isolated from mosquitoes in the United States and Trinidad. Virology 2009,386(1):154-159. 10.1016/j.virol.2008.12.034
Morales-Betoulle ME, Monzon Pineda ML, Sosa SM, Panella N, Lopez MR, Cordon-Rosales C, Komar N, Powers A, Johnson BW: Culexflavivirus isolates from mosquitoes in Guatemala. J Med Entomol 2008,45(6):1187-1190. 10.1603/0022-2585(2008)45[1187:CFIFMI]2.0.CO;2
Farfan-Ale JA, Lorono-Pino MA, Garcia-Rejon JE, Hovav E, Powers AM, Lin M, Dorman KS, Platt KB, Bartholomay LC, Soto V, 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(1):85-95.
Saiyasombat R, Dorman KS, Garcia-Rejon JE, Lorono-Pino MA, Farfan-Ale JA, Blitvich BJ: Isolation and sequence analysis of Culexflavivirus from Culex interrogator and Culexquinquefasciatus in the Yucatan Peninsula of Mexico. Arch Virol 2010,155(6):983-986. 10.1007/s00705-010-0665-5
Cook S, Moureau G, Harbach RE, Mukwaya L, Goodger K, Ssenfuka F, Gould E, Holmes EC, de Lamballerie X: Isolation of a novel species of flavivirus and a new strain of Culexflavivirus (Flaviviridae) from a natural mosquito population in Uganda. J Gen Virol 2009,90(Pt 11):2669-2678.
Wang H, Li Y, Liang X, Liang G: Japanese encephalitis in mainland china. Jpn J Infect Dis 2009,62(5):331-336.
Gao X, Nasci R, Liang G: The neglected arboviral infections in mainland China. PLoS Negl Trop Dis 2010,4(4):e624. 10.1371/journal.pntd.0000624
Wang ZS AS, Wang Y, Han R, Guo JQ: A new virus of flavivirus: Chaoyang virus isolated in Liaoning province. Chin J Public Health 2009,25(7):3.
Lu BL: Identification of mosquitoes. In Classification and Identification of Important Medical Insects of China. Edited by: Lu, BL, Wu, YH. Zhengzhou: Henan Science and Technology Publishing House; 2003:12-105. in Chinese
Li WJ, Wang JL, Li MH, Fu SH, Wang HY, Wang ZY, Jiang SY, Wang XW, Guo P, Zhao SC, et al.: Mosquitoes and mosquito-borne arboviruses in Qinghai-Tibet Plateau-Focus on the Qinghai area, China. Am J Trop Med Hyg 2010,82(4):705-711. 10.4269/ajtmh.2010.09-0649
Liang GD, He Y, Chen BQ, Zhao ZJ, Huang YJ, Chen L: Preparation of arbovirus group-specific PcAb and their use to identify newly isolated virus strains. Zhonghua Shi Yan He Lin Chuang Bing Du Xue Za Zhi 1993,7(4):374-376. in Chinese
Thompson JD, Gibson TJ, Higgins DG: Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics 2002., Chapter 2: Unit 2 3
Guindon S, Delsuc F, Dufayard JF, Gascuel O: Estimating maximum likelihood phylogenies with PhyML. Methods Mol Biol 2009, 537: 113-137. 10.1007/978-1-59745-251-9_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(6):881-888. 10.1080/10635150290102582
We thank the staff from Shandong Province Centers for Disease Control and Prevention and the two local Centers for Disease Control and Prevention (Heze city and Dongmingcounty) for assistance with the mosquito collection.
This work was supported by grants from the Ministry of Science and Technology of People's Republic of China (No. 2003BA712A08-01 and No. 2008ZX10004-010), the National Natural Science Foundation of China (31070145), Development Grant of State Key Laboratory for Infectious Disease Prevention and Control (2008SKLID105).
The authors declare that they have no competing interests.
Wang Huanyu1, Wang Haiyan2, Fu Shihong 1, Liu Guifang2, Liu Hong1, Gao Xiaoyan1, Song Lizhi2, Simon Rayner3, Xu Aiqiang2, Liang Guodong1*
WHYu carried out mosquito collection, virus isolation, nucleic acid detection and sequencing, participated in the sequence alignment, phylogenetic analysis and drafted the manuscript. WHYan, LGF, SLZ and XAQ participated in the mosquito collection. FSH, LH and GXY participated in the virus isolation, nucleic acid detection and phylogenetic analysis. SR participated in discussions and data analysis. LGD conceived the study and participated in its design and coordination. All authors read and approved the final manuscript.
About this article
Cite this article
Huanyu, W., Haiyan, W., Shihong, F. et al. Isolation and identification of a distinct strain of Culex Flavivirus from mosquitoes collected in Mainland China. Virol J 9, 73 (2012). https://doi.org/10.1186/1743-422X-9-73