Genome sequencing and genetic characterization of Culex Flavirirus (CxFV) provides new information about its genotypes
© The Author(s). 2016
Received: 13 April 2016
Accepted: 14 September 2016
Published: 23 September 2016
Culex Flavivirus (CxFV) is an insect-specific virus that is widely distributed and primarily infects mosquito species from the genus Culex. Its hosts include Culex tritaeniorhynchus, Culex quinquefasciatus, and Anopheles sinensis mosquitoes. Since its original identification, CxFV has been reported in several countries. Despite the increasing number of reports on CxFV, little is known about its genomic characteristics. It is unclear whether the phylogenetic relationships between the strains are influenced by host species and geographic location.
We characterized the Brazilian CxFV strain and performed a comprehensive genetic and phylogenetic characterization of CxFV based on all ORF sequences described so far. Our results revealed that the Brazilian strain is in a monophyletic clade with the Mexican strain. Overall, selective pressure indicates that the ORF is undergoing purifying selection.
The phylogenetic analysis revealed a strong association between climate and CxFV ancestry. Also, based on phylogeny and the genetic distance between the main branches of the tree, we propose the classification of the available sequences into two different genotypes. We also suggest the existence of two different subtypes within Genotype 1.
Flaviviruses are well known for causing important mosquito-borne human diseases. They are transmitted between arthropods and vertebrates, and are capable of replicating in both hosts. In 1975, a new member of this genus, incapable of replicating in vertebrate cells, was isolated from mosquito cell cultures and named Cell Fusing Agent Virus (CFAV) . More recently, the isolation of other insect-specific viruses has been reported [2–5]. One of these insect-specific viruses is Culex Flavivirus (CxFV), an insect-specific virus that was first identified in 2007 in Culex pipiens mosquitoes . Other hosts include Culex tritaeniorhynchus, Culex quinquefasciatus, and Anopheles sinensis mosquitoes [5, 7, 8].
Similar to other flaviviruses, CxFV encodes a polyprotein from a single-strand positive RNA open reading frame (ORF), flanked by 3′ and 5′-untranslated regions (UTR). The polyprotein of 3,364 aa is cleaved during and after translation into structural and non-structural proteins in the following order: C-prM(M)-E-NS1-NS2A-NS2B-NS3-NS4A-2K-NS4B-NS5 .
Since its identification, Culex Flavivirus has been reported in the United States, Mexico, Guatemala, Trinidad, Uganda, Indonesia, Japan, China, and Brazil, indicating that it is widely distributed [5, 6, 9–14]. Despite the increasing number of CxFV sequences deposited in GenBank, it is not clear if the phylogenetic relationships between the strains are influenced by the host species or geographic location. Moreover, little is known about the genetic profile of the CxFV genome.
Previously, we reported the first identification of CxFV in Brazil . Here, we perform a genetic and phylogenetic analysis of the complete ORF of the Brazilian CxFV strain. In addition, we perform the genetic characterization of CxFV based on the ORF sequences described so far, and we propose a novel genotype classification.
Culex Flavivirus (CxFV) strain
According to report published previously by this research group , The CxFV strain obtained in this study was derived from Culex quinquefasciatus mosquitoes, collected between April 2007 and January 2008, in the city of São José do Rio Preto (SJRP). The mosquitoes were pooled according to genus, gender and geographic location, for RNA extraction using Trizol protocol (Invitrogen).
The presence of Flavivirus was detected for each pool by Multiplex-Nested_RT-PCR . Ten positives pools were inoculated three times passages on C6/36. After the third passage pools were tested with primers specific for Culex Flavivirus (CxFV), Saint Louis Encephalitis Virus (SLEV), and West Nile Virus (WNV) . One CxFV isolated was used for sequence analysis in this manuscript.
Next generation sequencing
The CxFV RNA extracted from the C6/36 supernatant was quantified using a UV/Vis spectrophotometer (PicoDrop P200, Picodrop, Hinxton, Cambridgeshire, United Kingdom) and treated with Deoxyribonuclease (DNase I, Invitrogen, Carlsbad, CA, USA) to digest single and double-stranded DNA. cDNA was obtained using commercial SuperScript® III Reverse Transcriptase (Thermo Scientific, San Jose, CA, USA) according to the manufacturer’s protocol. A multiplex nested PCR was used once again to confirm CxFV positive sample [14, 15].
Double-stranded DNA was obtained using the SuperScript® Double-Stranded cDNA Synthesis Kit (Thermo Scientific, San Jose, CA, USA), purified with commercial magnetic beads (Agencourt AMPure XP, Beckman Coulter Inc., Indianapolis, IN, USA), and quantified. Sample were processed with the Nextera XT DNA Library Preparation Kit and sequenced using a commercial kit (MiSeq System, Illumina Inc., San Diego, CA, USA).
The assembly was performed by Geneious R6 (Biomatters, Auckland, New Zealand), using sequence of Toyama740/2005 CxFV strain (gb AB701772) as reference to assemble it.
Culex flavivirus sequences
All the analyses in this study were based on a 26-sequence dataset composed of all complete Culex Flavivirus ORF sequences (10,092 nucleotides) deposited in GenBank at the time of the study, plus the Brazilian sequence described here (KT726939). Their accession numbers are: AB701776.1, AB701775.1, AB701774.1, AB701773.1, AB701772.1, AB701771.1, AB701770.1, AB701769.1, AB701768.1, AB701767.1, AB701766.1, JQ308190.1, JQ308189.1, JQ308188.1, JQ308187.1, JQ308186.1, HQ678513.1, FJ663034.1, GQ165808.1, FJ502995.1, NC_008604.2, JQ518484.1, EU879060.1, AB377213.1, and AB262759.2. All Sequences were aligned using Muscle in the SeaView 4.4.2 package [16, 17].
Similarity analysis was carried out using the Sequence Identity Matrix tool from the BioEdit 184.108.40.206 package .
Mean genetic distances for the CxFV ORF and each specific region were calculated by the Kimura 2-parameters substitution model using Mega 6 .
The site-by-site genetic distances between CxFV_BR-RP01/2007 (KT726939) and the sequences identified in other countries were calculated in SimPlot 3.5.1 . The analysis was carried out using the Kimura 2-parameter distance model and a 200-bp sliding window. Sequences were grouped by geographic location.
Selective pressure analysis
Global and site-by-site ω were estimated in HyPhy . Global ω was calculated for the full ORF and for each genomic region using the HKY85_3x4 substitution model and global parameters. Site-specific selection was calculated Single Likelihood Ancestor Counting (SLAC) method, the MG94xHKY85 model, and a p-value of 0.05.
As HyPhy requires a starting tree for all analysis, a maximum-likelihood tree was constructed for the entire ORF and for specific regions by PhyML in the SeaView package . The GTR model was used, with optimized nucleotide equilibrium frequencies and site rate variation.
A Bayesian Markov Chain Monte Carlo (MCMC) analysis was used to estimate the phylogenetic tree using BEAST v1.8 . The analysis was performed under strict molecular clock using as tree prior a constant population size coalescent and MCMC was run for 10,000,000 steps and sampled every 1000 steps. Trees were summarized using Tree Annotator 1.8 available on Beast package and edited using FigTree 1.4 . BaTS 0.9 was used to statistically test for association between tree topology and geography, host species or climate .
The mean genetic distance between the sequences in specific branches of the phylogenetic tree was calculated in Mega 6 . Sequences were divided into groups according to their distribution in the phylogenetic tree and the mean base substitution per site between groups was conducted using the Kimura 2-parameters substitution model .
Culex Flavivirus (CxFV) identification
After collection, the mosquitoes were identified by gender and grouped into pools. A total of 83 pools of female Culex quinquefasciatus, ten mosquitoes per pool, were analyzed. Among the samples, 41 (49.4 %) were positive for CxFV in multiplex RT-PCR and ten were isolated in C6/36. After three passages in C6/36 cells isolation was confirmed from pool 73C by Multiplex-RT-PCR. The complete genome of the isolate from pool 73C was perform by Illumina Plataform.
Next generation sequencing
NGS by Illumina generated a total of 69,242 reads of coverage, with 60,405 of CxFV reads, that generated 10,706 nucleotides. The sequence obtained was analyzed on BLAST, for confirmation of its contents and was deposited on GenBank under the accession number KT726939 and is here referred as CxFV_BR-RP01/2007.
The sequence of the open reading frame (ORF) of the Brazilian strain was aligned with all available CxFV full genome sequences deposited on GenBank. The percentage of identity between the sequences was determined with the BioEdit package using CxFV_BR-RP01/2007 as the query sequence . The highest similarity was observed between the Brazilian sequence and the sequence from Mexico (98.5 %) followed by that of Uganda (98.2 %) (Additional file 1: Table S1). All other sequences were approximately 90 % identical to CxFV_BR-RP01/2007.
Selective pressure analysis
In order to test for association of phylogeny with either host, geographical location or climate the posterior set of trees estimated by BEAST were analyzed by BaTS, which quantifies and statistically tests the phylogeny-trait correlation. Evidence of strong association was detected between phylogeny and climate (p < 0.05). Also, Association index (AI) and parsimony score (PS) indicate association for both host and geographical location (p < 0.05). However monophyletic clade (MC) significance reveals that it only applies for Japan and China (location) and Culex quinquefasciatus (host) (Additional files 2–4: Tables S2–S4).
In order to estimate the evolutionary divergence between the main clades of the phylogenetic tree, the mean genetic distances between these groups were calculated. Our results show that Clade 1 and Clade 2 diverge 10.4 %, while the two main branches that compose Clade 1 diverge 5.1 %.
Insect-specific flaviviruses, including CxFV, are a highly divergent group within the Flavivirus genus that shares a common ancestor with all other members, including the disease-causing ones. Understanding the genetic characteristics of CxFV is important not only to understand its biology and evolutionary history, but also to clarify some aspects of the whole genus.
Culex Flavivirus primarily infects globally distributed mosquito species of the genus Culex, which are vectors for pathogenic flaviviruses like WNV, SLEV, and JEV [27–29]. The most represented hosts in this work, Culex pipiens and Culex quinquefasciatus, are members of the Culex pipiens complex that consists of closely related species that are difficult to distinguish morphologically . Studying the geographical distribution of the host is important to understand the evolutionary relationships between the identified CxFV strains. Two main monophyletic branches, Clade 1 and Clade 2, can be observed in the phylogenetic tree with high branch support. Regarding countries, all Asian strains are grouped in Clade 1 along with some sequences from USA. Clade 2 is composed by sequences from different countries with only one sequence each. When considering host species, Clade 2 is composed only of sequences derived from Culex quinquefasciatus while Clade 1 as sequences from different hosts. This scenario prevents a clear association by simply analyzing the ancestry. Statistical analysis however revealed significant association between phylogeny and climate (temperate x tropical). Clade 1 consists of virus from temperate climate regions while Clade 2 is composed from sequences from region of tropical climate. Climate association is related with host species distribution. Although Culex pipiens occurs in temperate regions and Culex quinquefasciatus in tropical and subtropical regions their range overlaps and it has been proven they can hybridize [30–32]. When testing for host-phylogeny association only Culex quinquefasciatus was significantly associated. Moreover, country association was only found for Japan and China. Considering that both country and host are related to climate these factors may also be associated with CxFV ancestry, however the limited number of sequences from some hosts and geographic regions is probably preventing a more consistent conclusion.
We analyzed the genetic distance between the main branches of the phylogenetic tree that was reconstructed based on all CxFV ORF sequences described so far, including the Brazilian strain described in this work. Clades 1 and 2 consist of strong groups with robust branch support, which are 10.4 % genetically distant. These observations sustain the idea that they could be two different genotypes of CxFV. For Japanese Encephalitis Virus, another Culex borne virus, genotypes are defined by 12 % difference in nucleotide composition of the highly divergent PrM gene [33, 34]. Genotypes of DENV, also a flavivirus, are defined by 6 % difference within a serotype based on the E region [35, 36]. In both cases, the limits were determined arbitrarily but supported by strong phylogenetic evidence such as ours. Other researchers have previously suggested that CxFV should be classified into different genotypes, based on the phylogenetic analysis of the E genomic region . According to those authors, Clade 1 described here would consist of Genotype 1 and Clade 2 of Genotype 2. Our analysis also evidenced that the two main branches, with strong branch support (1000), in Clade 1 diverge 5.1 %. These two branches could be considered different subtypes of CxFV Genotype 1 (1a and 1b). Similar results were obtained using the E region (Additional file 5: Table S5), suggesting it is a suitable region for CxFV genotyping and subtyping. We propose that groups of CxFV that are more than 9 % genetically distant could be classified as different genotypes and that differences higher than 4 % within these groups could be considered different subtypes. Although the cut-off values are arbitrarily defined, based on the genetic data available so far, they seem suitable.
Overall, selective pressure analysis indicates that the entire ORF, as well as specific regions, are undergoing purifying selection. Accordingly, site-by-site analysis found no positively selected sites and 106 sites under negative selection. This suggests that the virus is well adapted to its host and evolutionary forces are working against amino acid change. The infection of Culex mosquitoes by CxFV apparently does not result in disease. This observation is supported by the fact that the infection of C6/36 and AeAl-2, two Aedes albopictus cell lines, results only occasionally in moderate cytopathic effects [5, 6, 12]. If this is the case, any change in the genomic composition could disrupt the fine balance between what is advantageous for the virus and at the same time does not negatively affect the host, supporting the purifying selection.
Although the differences are not substantial, the NS2A region presents higher ω values and lower percentages of negatively selected sites. It is interesting, however, to notice that the same region also showed the lowest genetic distance. At first, these data might seem in disagreement. However, it suggests that NS2A displays low genetic variability in general, but a more relaxed pressure against amino acid change when compared to other regions. CxFV and other Flaviviruses from the insect-specific group present a frameshift signal that produces an alternate reading frame called fifo . The signal consists of a stem-loop structure in the RNA genome that is harbored inside the NS2A region . In this case, the selection is acting not only on the protein, but also on the RNA structure to avoid disturbing the frameshift signal. It is important to highlight that different amino acids might have similar chemical properties and a change would not necessarily disrupt the tertiary structure. This could explain why, despite the low genetic distance, this region is more permissive to amino acid change.
There are still many aspects of Culex Flavivirus and other insect-specific flaviviruses that needs to be clarified. Culex ssp. mosquitoes are distributed worldwide suggesting that CxFV should be present in more places than currently reported. Here we report an association between the ancestry of CxFV and climate. The identification of new strains from different countries and hosts is important for elucidating the ecology and evolutionary history of this virus. Based on the available sequences, we propose the classification of Culex Flavivirus into two genotypes (1 and 2) and the existence of two subtypes within genotype 1.
The authors acknowledge CNPq, CAPES and FAPESP for the financial support of MLN (2013/21719-3), CB (165802/2015-4; 2012/028560) and DV (PNPD/Microbiologia 1278383/2012).
CB carried out the genome analysis, analyzed the data and wrote the first draft of the manuscript. DCM isolated the virus and participate in the genome sequencing and analysis of the data. DV and LSU carried out the genome sequence and genome annotation, and analysis of the data. PR, JPAJ, MLN designed the study, analyzed the data and reviewed the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis 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.
- 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.View ArticlePubMedGoogle Scholar
- Huhtamo E, Moureau G, Cook S, Julkunen O, Putkuri N, Kurkela S, Uzcategui NY, Harbach RE, Gould EA, Vapalahti O, de Lamballerie X. Novel insect-specific flavivirus isolated from northern Europe. Virology. 2012;433:471–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Hobson-Peters J, Yam AW, Lu JW, Setoh YX, May FJ, Kurucz N, Walsh S, Prow NA, Davis SS, Weir R, et al. A new insect-specific flavivirus from northern Australia suppresses replication of West Nile virus and Murray Valley encephalitis virus in co-infected mosquito cells. PLoS One. 2013;8:e56534.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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:85–95.PubMedPubMed CentralGoogle Scholar
- 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 Culex pipiens mosquito in Japan. Virology. 2007;359:405–14.View ArticlePubMedGoogle Scholar
- Wang Z, Guo J, An S: JQ308188.1 - Culex flavivirus was isolated from Culex pipiens and Anopheles in liaoning province. In Nucleotide. GenBank: NCBI; 2011.Google Scholar
- Obara-Nagoya M, Yamauchi T, Watanabe M, Hasegawa S, Iwai-Itamochi M, Horimoto E, Takizawa T, Takashima I, Kariwa H. 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.View ArticlePubMedGoogle Scholar
- Blitvich BJ, Lin M, Dorman KS, Soto V, Hovav E, Tucker BJ, Staley M, Platt KB, Bartholomay LC. 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 Culex flavivirus (Flaviviridae) from a natural mosquito population in Uganda. J Gen Virol. 2009;90:2669–78.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim DY, Guzman H, Bueno Jr R, Dennett JA, Auguste AJ, Carrington CV, Popov VL, Weaver SC, Beasley DW, Tesh RB. Characterization of Culex Flavivirus (Flaviviridae) strains isolated from mosquitoes in the United States and Trinidad. Virology. 2009;386:154–9.View ArticlePubMedGoogle Scholar
- Morales-Betoulle ME, Monzon Pineda ML, Sosa SM, Panella N, Lopez MR, Cordon-Rosales C, Komar N, Powers A, Johnson BW. Culex flavivirus isolates from mosquitoes in Guatemala. J Med Entomol. 2008;45:1187–90.View ArticlePubMedGoogle Scholar
- 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. Arch Virol. 2010;155:983–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Machado DC, Mondini A, dos Santos Santana V, Yonamine PT, Chiaravalloti Neto F, Zanotto PM, Nogueira ML. First identification of Culex flavivirus (Flaviviridae) in Brazil. Intervirology. 2012;55:475–83.View ArticlePubMedGoogle Scholar
- de Morais Bronzoni RV, Baleotti FG, Ribeiro Nogueira RM, Nunes M, Moraes Figueiredo LT. Duplex reverse transcription-PCR followed by nested PCR assays for detection and identification of Brazilian alphaviruses and flaviviruses. J Clin Microbiol. 2005;43:696–702.View ArticlePubMedGoogle Scholar
- Gouy M, Guindon S, Gascuel O. SeaView version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010;27:221–4.View ArticlePubMedGoogle Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.Google Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999;73:152–60.PubMedPubMed CentralGoogle Scholar
- Pond SL, Frost SD, Muse SV. HyPhy: hypothesis testing using phylogenies. Bioinformatics. 2005;21:676–9.View ArticlePubMedGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.View ArticlePubMedGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.View ArticlePubMedPubMed CentralGoogle Scholar
- FigTree 1.4 [http://tree.bio.ed.ac.uk/software/figtree/].
- Parker J, Rambaut A, Pybus OG. Correlating viral phenotypes with phylogeny: accounting for phylogenetic uncertainty. Infect Genet Evol. 2008;8:239–46.View ArticlePubMedGoogle Scholar
- Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16:111–20.View ArticlePubMedGoogle Scholar
- WHO. West Nile virus. In: Fact sheet N°354. Geneva: Word Health Organization; 2011.Google Scholar
- WHO. Japanese encephalitis. In: Fact sheet No 386. Geneva: World Health Organization; 2014.Google Scholar
- CDC. Saint Louis Encephalitis. In: Technical Fact Sheet. Atlanta: Centers for Disease Control and Prevention; 2009.Google Scholar
- Farajollahi A, Fonseca DM, Kramer LD, Marm Kilpatrick A. “Bird biting” mosquitoes and human disease: a review of the role of Culex pipiens complex mosquitoes in epidemiology. Infect Genet Evol. 2011;11:1577–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Fonseca DM, Smith JL, Kim HC, Mogi M. Population genetics of the mosquito Culex pipiens pallens reveals sex-linked asymmetric introgression by Culex quinquefasciatus. Infect Genet Evol. 2009;9:1197–203.View ArticlePubMedPubMed CentralGoogle Scholar
- Fonseca DM, Smith JL, Wilkerson RC, Fleischer RC. Pathways of expansion and multiple introductions illustrated by large genetic differentiation among worldwide populations of the southern house mosquito. Am J Trop Med Hyg. 2006;74:284–9.PubMedGoogle Scholar
- Solomon T, Ni H, Beasley DW, Ekkelenkamp M, Cardosa MJ, Barrett AD. Origin and evolution of Japanese encephalitis virus in southeast Asia. J Virol. 2003;77:3091–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen WR, Tesh RB, Rico-Hesse R. Genetic variation of Japanese encephalitis virus in nature. J Gen Virol. 1990;71(Pt 12):2915–22.View ArticlePubMedGoogle Scholar
- Rico-Hesse R. Molecular evolution and distribution of dengue viruses type 1 and 2 in nature. Virology. 1990;174:479–93.View ArticlePubMedGoogle Scholar
- Rico-Hesse R. Microevolution and virulence of dengue viruses. Adv Virus Res. 2003;59:315–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Huanyu W, Haiyan W, Shihong F, Guifang L, Hong L, Xiaoyan G, Lizhi S, Rayner S, Aiqiang X, Guodong L. Isolation and identification of a distinct strain of Culex Flavivirus from mosquitoes collected in Mainland China. Virol J. 2012;9:73.View ArticlePubMedPubMed CentralGoogle Scholar
- Firth AE, Blitvich BJ, Wills NM, Miller CL, Atkins JF. Evidence for ribosomal frameshifting and a novel overlapping gene in the genomes of insect-specific flaviviruses. Virology. 2010;399:153–66.View ArticlePubMedPubMed CentralGoogle Scholar
- Blitvich BJ, Firth AE. Insect-specific flaviviruses: a systematic review of their discovery, host range, mode of transmission, superinfection exclusion potential and genomic organization. Viruses. 2015;7:1927–59.View ArticlePubMedPubMed CentralGoogle Scholar