Asian Pacific Journal of Tropical Medicine

: 2019  |  Volume : 12  |  Issue : 5  |  Page : 216--223

Phylogeny of Culex theileri virus flavivirus in Spain, Myanmar, Portugal and Turkey

Eleonora Cella1, Domenico Benvenuto1, Daniele Donati2, Francesco Garilli3, Silvia Angeletti4, Stefano Pascarella5, Massimo Ciccozzi1,  
1 Unit of Medical Statistics and Molecular Epidemiology, University Campus Bio-Medico of Rome, Rome, Italy
2 School of Nursing, Faculty of Medicine, Department of Biomedicine and Prevention, Tor Vergata University, Via Montpellier 1, 00133, Rome, Italy
3 Faculty of Medicine, University Campus Bio-Medico of Rome, Rome, Italy
4 Unit of Clinical Laboratory Science, University Campus Bio-Medico of Rome, Italy
5 Department of Biochemical Sciences “A. Rossi Fanelli”, University of Rome “La Sapienza”, 00185 Roma, Italy

Correspondence Address:
Massimo Ciccozzi
Unit of Medical Statistics and Molecular Epidemiology, University Campus Bio-Medico of Rome, Rome


Objective: To study the genetic diversity of Culex theileri flavivirus and the spread of this virus among Spain, Portugal and Turkey. Methods: A database consisting of 55 sequences of the NS5/3’UTR region of Culex theileri flavivirus group downloaded from GenBank were aligned and manual edited with Bioedit. ModelTest v. 3.7 was used to select the simplest evolutionary model that adequately fitted the sequence data. Maximum likelihood analysis was performed using MEGA7. The phylogenetic signal of the dataset was investigated by the likelihood mapping analysis. Results: The phylogenetic tree showed three clusters. Myanmar sequences clusterd together with Turkish sequences, Spain and Portugal strains grouped together and two Turkish sequences grouped separately. Selective pressure analysis showed a moderate percentage of sites (22.5%) under pervasive negative selection and only 1% under pervasive positive selection. The sites subject to selective pressure in CTFV RdRp NS5 fragments have been located onto the predicted three-dimensional structure. Conclusions: Phylogenetic and evolutionary analysis can be an important tool for understanding the evolutionary impact of the probable contemporary existence between non-pathogenic and pathogenic flaviviruses among these vectors.

How to cite this article:
Cella E, Benvenuto D, Donati D, Garilli F, Angeletti S, Pascarella S, Ciccozzi M. Phylogeny of Culex theileri virus flavivirus in Spain, Myanmar, Portugal and Turkey.Asian Pac J Trop Med 2019;12:216-223

How to cite this URL:
Cella E, Benvenuto D, Donati D, Garilli F, Angeletti S, Pascarella S, Ciccozzi M. Phylogeny of Culex theileri virus flavivirus in Spain, Myanmar, Portugal and Turkey. Asian Pac J Trop Med [serial online] 2019 [cited 2019 Jun 19 ];12:216-223
Available from:

Full Text

 1. Introduction

The genus flavivirus comprises more than 70 RNA virus species that include Yellow fever virus, Dengue virus, Japanese encephalitis virus, and the Tick-borne encephalitis virus complex. Many of these arthropod-borne viruses represent dangerous threats to human health and have been subjected to intensive research to unravel their molecular and virological properties[1].

Insect-specific flaviviruses (ISFs) are distinct for the inability to be propagated in vertebrate cell lines and lack of vertebrate hosts[2]. ISFs have been recognized as a major phylogenetic group along with mosquito-borne, tick-borne and no-known-vector flaviviruses[3],[4],[5], ISFs share many similar properties with widely known pathogenic arthropod-borne flaviviruses such as dengue virus (DENV), West Nile virus (WNV) and tick-borne encephalitis virus (TBEV). They are enveloped, single-stranded RNA viruses, with a single open reading frame (ORF) that encodes the viral polyprotein, which is then cleaved by viral and host proteases to form three structural (C, preM and E), and seven nonstructural (NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5) proteins[6]. ISFs seem to have a global distribution and have been detected in field-collected mosquitoes from Europe, East Asia, Africa, the American Continent and Australia[2],[7]. The ISFs possess potentially distinct and unique features among flaviviruses, which are likely to have an impact on the circulation of their related vector-borne pathogens. Their widespread distribution without obvious reservoirs or amplification hosts requires exploration of potential mechanisms, such as vertical transmission, contributing to their dispersion[7],[8],[9]. In particular, Culex (Cx.) theileri flavivirus (CTFV) seems to affect the flight activity of Culex mosquitos[10].

The aim of this study was to investigate the evolutionary patterns of Cx. theileri virus using all the available sequences from National Center for Biotechnology Information (NCBI). It has performed a phylogenetic analysis including phylogenetic tree and homology model analysis.

 2. Materials and methods

2.1. Sequence dataset

A dataset were built based on previous investigations about Cx. theileri flavivirus species made on the NCBI website, the previous investigation found a clade consisting of 55 sequences of NS5/3’ UTR region of Cx. theileri flavivirus which has been selected and downloaded from GenBank ( genbank/). All the sequences were aligned and manually edited using the Bioedit program v7.2.5, as already described[11].

2.2. Phylogenetic analysis

The phylogenetic signal was investigated by means of the likelihood mapping analysis of 10 000 random quartets by using TreePuzzle program as already described[12]. In this analysis, groups of four randomly chosen sequences (quartets) were evaluated using Maximum Likelihood. For each quartet, the three possible unrooted trees were reconstructed under the selected substitution model. The likelihoods of each tree were then plotted on a triangular surface, so that fully resolved trees fall into the corners and the unresolved quartets in the center of the triangle (indicating a star-like signal). When using this strategy, if more than 30% of the dots fall into the center of the triangle, the data are considered unreliable for the purposes of phylogenetic inference[12].

The evolutionary model was chosen, as the best-fitting nucleotide substitution model in accordance with the results of the hierarchical likelihood ratio test implemented in MODELTEST software (version 3.7)[13]. The Maximum Likelihood tree was reconstructed using MEGA7 with Tamura-Nei model as evolutionary model selected by ModelTest[13],[14],[15].

Adaptive Evolution Server ( was used to find eventual sites of positive or negative selection using FEL (Fixed Effects Likelihood) which uses a maximum-likelihood (ML) approach to infer nonsynoymous (dN) and synonymous (dS) substitution rates on a per-site basis for a given coding alignment and corresponding phylogeny[16]. Sites were considered to have been subjected to statistically significant positive or negative selection based on P value < 0.05 (reference sequence used, accession number: KX652375.1).

Possible spreading ways from Spain, Portugal, Myanmar and Turkey have been analyzed using the website http://portugalturkey. com/ and for a better understanding of possible spreading paths of the vector of the virus.

2.3. Homology modelling

Sequence alignments and analyses were obtained through the Jalview editor[17] and the programs ClustalO[18] or Muscle[19],[20]. Structural templates homologous to the sequences studied within this work have been identified with HHpred[20] and SwissModel[21] template search tools. Homology models were built with the program Modeller[22] and validated with QMEAN[23] and Prosa II [24] software tools. Three-dimensional structures were displayed and studied with PyMOL (The PyMOL Molecular Graphics System, Version 2.0 Schroedinger, LLC) or Chimera[25]. Sequence conservation and mapping analysis has been carried out with the Consurf web tool[26].

 3. Results

3.1. Phylogenetic analysis

The percentage of dots falling in the central area of the triangle was 4.6%, the dataset did not show more than 30% of noise and contained sufficient phylogenetic signal [Figure 1].{Figure 1}

The phylogenetic tree showed three clusters [Figure 2]. The cluster I included 28 strains isolated from Myanmar and Turkey. Cluster II included only two sequences from Turkey. The cluster III included 25 strains isolated from Spain and Portugal.{Figure 2}

Clade I was composed by two subclusters: 18 sequences all isolated in 2014 in Myanmar grouped in cluster I a and cluster I b included 10 sequences of Cx. theileri flavivirus isolated from 2014 to 2015 in Turkey.

Cluster II include only two Turkish sequences isolated in 2014 and 2015. Cluster III consisted of two different clusters (IIIa and IIIb): III a included Spanish isolates (in 2006) with the sequences isolated in Portugal in 2012 and IIIb strains isolated in Spain and Portugal respectively from 2003 to 2006 and in 2012 [Figure 2].

FEL analyzed 297 sites and found evidence of pervasive positive selection in 2 sites with P<0.05, it has been evidence of pervasive negative selection in 67 sites (22.5%) with P<0.05.

3.2. Homology modelling

A multiple sequence alignment among the RNA-dependent RNA polymerase (RdRp) fragments of CTFV NS5 proteins and the structurally solved NS5 sequences from flavivirus has been obtained [Figure 3]. According to the alignment, the CTFV RdRp fragments here considered cover the region 488-798 of the full-length Japanese encephalitis virus NS5 (PDB code 4K6M). This portion [Figure 4] contains the “Palm”, “Middle” and part of the “Thumb” subdomains of the human right hand architecture characteristic of the NS5 proteins[27]. To map the structural variability of the variant CTFV RdRp fragments and their sites under selection pressure, homology modelling has been applied. Among the available structural templates, the PDB structure with code 4HDG corresponding to the crystal structure of the RNA-dependent RNA polymerase domain from Japanese encephalitis Virus NS5 in complex with GTP, has been selected. The decision has been taken in consideration of the suitable quality of the structure in terms of resolution (2.38 Å) and R-value (0.276), for the presence of the ligand GTP, and for the high sequence identity to the targets (about 53%). The sequence heterogeneity of the set of CTFV RdRp sample suggested to build models for at least two variants each representing a sequence subgroup [Figure 3]. To this purpose, the fragments denoted by the codes HE997074 and KU958176 in Figure 3 have been selected. Models created from the two fragments do not show significant differences except for conformations of exposed loops [Figure 5]; for that reason, the HE997074 model was used as reference.{Figure 3}{Figure 4}{Figure 5}

The sequence alignment in [Figure 3] hints at the presence of several structural differences between the CTFV RdRp fragment set and the corresponding NS5 structures from flaviviruses. In particular, two deletions in the CTFV fragment set correspond to Gln residues in positions 528, and position intervals 636-644 of 4K6M and 4HDG sequences [Figure 6]. In addition, two 2-residue insertions are between sequence positions 33-34 and 62-63 of the CTFV RdRp as reported in [Figure 3], which correspond to positions 557-558 and 584-585 of 4K6M and 4HDG sequences [Figure 5].{Figure 6}

The sites subject to selective pressure in CTFV RdRp NS5 fragments have been located onto the predicted three-dimensional structure [Figure 6]. The two sites are at columns 44 and 134 of the alignment displayed in [Figure 3] which correspond to positions 567and 655 of 4K6M and 4HDG structures.

Position 44 (567 of 4K6M) displays Ala or Ile residues in the flaviviruses sequences whereas it shows Ala in the structurally solved NS5 proteins. It is located within the -helix encompassed by positions 562-576 in 4K6M and it is buried with respect to the protein surface.

Position 134 in CTFV RdRp NS5 fragments (655 of 4K6M) displays Asn and, in two cases, His. The corresponding position in the structurally solved NS5 proteins has Asn or Lys. The position is within the “Palm” of the RdRp domain, exposed to the solvent, at a kink of a helix encompassed by 4K6M positions 642-662. In CTFV RdRp the N-terminal segment of the helix is deleted [Figure 3], [Figure 4] and replaced by a predicted loop.

ConSurf analysis conducted on the NS5 RdRp domain suggests that the position 44 is at a site conserved within homologous NS5 sequences from different flaviviruses [Figure 7]. The conservation may reflect the stabilizing role of the hydrophobic interaction taking place between the Ala/Ile and other hydrophobic side chains from the surrounding protein environment (not shown). The position 134 sits at a site with generally variable residues within homologous NS5 sequences [Figure 7].{Figure 7}

 4. Discussion

Flavivirus have a great impact on human health[28],[29]. Many viruses assigned to the flavivirus genus are exclusively replicated and detected in insects and named ISF and in last years, an increasing number of newly discovered ISFs probably due to the bio surveillance efforts regarding arthropod-borne viruses and new sequencing techniques were observed[30].

Studies on phylogenetic relationships of ISFs in relationship with others accepted flaviviruses are not clear, and they have previously been suggested to form a sort of divergent outgroup suggesting the possibility to represent an ancestral lineage of flaviviruses

Despite statistically supported separation of the strains originating from Portugal and Turkey based on CTFV NS5 RdRp gene region, current information indicates that these strains constitute local variants of CTFV. Although the most comprehensive characterization of the CTFV was performed from specimens collected in Portugal during 2009 and 2010, there is evidence for the circulation of several related and potentially identical strains[7]. Partial NS5 sequences closely related to CTFV have been detected in Spain and Portugal from 2006 to 2010 with different acronyms (CxthFV) or names (Spanish Culex flavivirus, SCxFV)[31],[32].

Conversely, some variations, especially found on the putative viral capsid, envelope and non-structural proteins NS1, NS2a, NS3, NS4b and NS5, are detected in strains from either location, which might be reflecting environmental adaptive processes.

Moreover, presence of partial CFTV NS5 sequences could also be detected in Culex pipiens pools in Turkey as well in Portugal and Spain on various occasions between 2006 and 2015[32],[33],[34].

Thus, the host range of CFTV may not be restricted to Culex theileri. Although a grouping of the ISFs according to their host mosquito species seemed plausible in initial analyses[35], several phylogenetically distinct ISFs have been repeatedly detected in various Culex and Aedes spp. Mosquitoes in different areas[7], suggesting independent introductions and multiple host-switching events[3]. As we can see from the data there is a consistent phylogenetic relationship among the presence of ISF in Spain, Portugal, Myanmar and Turkey, probably the infected Cx. theileri has spread by land from Spain to Portugal and from Turkey to Myanmar, vectored by mosquitos through commercial trading as so as happened with Aedes albopictus from Africa to Mediterranean areas in the recent years. Recently, passive surveillance activities to monitor mosquito species was performed in Europe since international trading based on tires is a well documented route of transmission of Aedes albopictus worldwide[36],[37].

Selective pressure analysis showed a moderate percentage of sites (22.5%) under pervasive negative selection and only two sites (1%) under pervasive positive selection. The first position under positive selective pressure has two amino acidic variants (Ile and Ala), interestingly the Turkish and Myanmar strain have the same residue (Ile) and Spanish and Portuguese strains only Ala residue. Similarly, in the other position under selective pressure the Turkish strains have only His residue and Myanmar, Spanish and Portuguese strains only Asp residue. These results suggest a different distribution and different spread of CTFV strains between these countries.

Protein sequence comparison highlighted several sequence differences between the CTFV RdRp fragments and the reference structurally solved RdRp domains from flaviviruses. Although it is very difficult if not impossible to predict the functional and structural impact of these variations, it may be argued that since insertions and deletions occur mostly in exposed loops, they may alter chain flexibility and interaction with other molecular components. Of particular interest can be considered the insertion in the CTFV RdRp “Middle finger” loop (positions 584-585 of 4HDG), which contributes to the interaction of the domain with the connected MTase module. It has been suggested that the interplay is relevant for NS5 catalytic process[38]. One of the positions (site 567 of 4HDG) under selective pressure occurs in a buried conserved site within the helix encompassed by positions 562-576 of 4K6M. Conservation of the hydrophobic characteristic (only Ala or Ile occur in the CTFV sequences) strongly suggests that this site is important for providing stability through van der Waals contacts with other hydrophobic side chains from the adjacent helix. The second position under selection, (site 655 of 4HDG) is exposed and close to the large deletions in CTFV RdRp fragment which shortens the -helix encompassed by the 4HDG positions 636-644.

The CTFV described in this phylogenetic and evolutionary analysis have been isolated from mosquitoes collected in Portugal during the summertime (2009-2010), so as from mosquitoes collected between 2001 and 2005 in the southeast and northeast of Spain. Turkey and Myanmar had virus isolation in 2014 and 2015. This was not unexpected because Spain and Portugal are neighbor countries in geological barriers between them, and in Turkey maybe for trading, between these three countries, as previously suggested[38, 39], indicating that this virus has a wider distribution and can be not restricted to a single mosquito genus or species. The different sites under selective pressure, underlining the different environmental pressure on the mosquito vectors living differently in Turkey and Iberian Peninsula, this pressure can be important, in causing a change of the vector for the virus as happened for other flaviviruses so as Chikungunya virus[28].

It has been found evidence of consistent commercial trades between Portugal and Turkey due to the development of trade agreements, these bilateral trades have started in 2014 that is the year of the first Cx. theileri flavivirus outbreak, moreover among the goods exported from Portugal to Turkey there are several potential hidings for mosquitoes (such as cables, cabinets, motor cars) and this is probably the way of spreading of Cx. theileri from Portugal to Turkey. Similarly there is a bilateral Economic and Commercial Relations between Turkey and Myanmar suggesting the possible way of virus spreading[39].

Cx. theileri are also considered as competent West Nile Virus (WNV) vectors among the mosquito species distributed throughout Europe. An interaction of CFTV and WNV in the mosquito host might be possible, referring to the nucleotide mutations for selective environmental pressure. The genome evolution probably can support the hypothesis that insect-specific flavivirus may adapt to a variety of vectors switching in different host. Phylogenetic and evolutionary analysis can be an important tool for understanding the evolutionary impact of the probable contemporary existence between non-pathogenic and pathogenic flaviviruses among these vectors. The molecular epidemiological surveillance is important to improve the knowledge about the interaction between pathogenic versus non-pathogenic viruses and host vector to control the pathogenic arboviruses epizootic in Mediterranean areas.

Conflict of interest statement

We declare that we have no conflict of interest


1Mukhopadhyay S, Kuhn RJ, Rossmann MG. A structural perspective of the flavivirus life cycle. Nat Rev Microbiol 2005; 3: 13-22.
2Calzolari M, Zé-Zé L, Vázquez A, Sánchez 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-388.
3Moureau G, Cook S, Lemey P, Nougairede A, Forrester NL, Khasnatinov M, et al. New insights into flavivirus evolution, taxonomy and biogeographic history, extended by analysis of canonical and alternative coding sequences. PLoS One 2015; 10: e0117849.
4Cook S, Holmes EC. Amultigene analysis of the phylogenetic relationships among the flaviviruses (Family: Flaviviridae) and the evolution of vector transmission. Arch Virol 2006; 151: 309-325.
5Cook 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-234.
6Simmonds P, Becher P, Collett MS, Gould EA, Heinz FX, Meyers G, et al. Family Flaviviridae virus taxonomy: Classification and nomenclature of viruses. ninth report of the international committee on taxonomy of viruses. San Diego: Elsevier; 2011, p. 1003-1020.
7Blitvich 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-1959.
8Crochu S, Cook S, Attoui H, Charrel RN, De Chesse R, Belhouchet M, et al. Sequences of flavivirus-related RNA viruses persist in DNA form integrated in the genome of Aedes spp. mosquitoes. J Gen Virol 2004; 85: 1971-1980.
9Cook 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-2678.
10Newman CM, Anderson TK, Goldberg TL. Decreased flight activity in Culex pipiens (Diptera: Culicidae) naturally infected with Culex flavivirus. J Med Entomol 2016; 53: 233-236.
11[11] 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-98.
12Schmidt HA, Petzold E, Vingron M, von Haeseler A. Molecular phylogenetics: Parallelized parameter estimation and quartet puzzling. J Parallel Distrib Comput 2003; 63: 719-727.
13Swofford D, Sullivan J. Lemey P, Salemi M, Vandamme AM. Phylogeny inference based on parsimony and other methods with PAUP. In: The phylogenetic handbook–A practical approach to phylogenetic analysis and hypothesis. 2nd ed. New York: Cambridge University Press; 2009, p. 267-312.
14Posada D, Buckley TR. Model selection and model averaging in phylogenetics: advantages of akaike information criterion and Bayesian approaches over likelihood ratio tests. Syst Biol 2004; 53: 793-808.
15Kumar S, Stecher G, Tamura K. MEGA7: Molecular evolutionary genetics analysis Version 7.0 for Bigger Datasets. Mol Biol Evol 2016; 33(7): 1870-1874.
16Pond SL, Frost SD, Muse SV. HyPhy: Hypothesis testing using phylogenies. Bioinformatics 2005; 21(5): 676-679.
17Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ. Jalview version 2: A multiple sequence alignment editor and analysis workbench. Bioinformatics 2009; 25(9): 1189-1191.
18Sievers F, Higgins DG. Clustal omega for making accurate alignments of many protein sequences. Protein Sci 2018; 27(1): 135-145.
19Edgar RC. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics 2004; 5: 113.
20Zimmermann L, Stephens A, Nam SZ, Rau D, Kubler J, Lozajic M, et al. A completely reimplemented mpi bioinformatics toolkit with a new HHpred server at its core. J Mol Biol 2018; 430(15): 2237-2243.
21Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res 2008; 46(W1): W296-W303.
22Eswar N, Webb B, Marti-Renom MA, Madhusudhan MS, Eramian D, Shen MY, et al. Comparative protein structure modeling using Modeller. Curr Protoc Bioinformatics 2006; Chapter 5:Unit-5.6. doi: 10.1002/0471250953.bi0506s15.
23Benkert P, Kunzli M, Schwede T. QMEAN server for protein model quality estimation. Nucleic Acids Res 2009; 37(Web Server issue): W510-514.
24Wiederstein M, Sippl MJ. ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res 2007; 35(Web Server issue): W407-410.
25Meng EC, Pettersen EF, Couch GS, Huang CC, Ferrin TE. Tools for integrated sequence-structure analysis with UCSF Chimera. BMC Bioinformatics 2006; 7: 339.
26Ashkenazy H, Abadi S, Martz E, Chay O, Mayrose I, Pupko T, et al. ConSurf 2016: An improved methodology to estimate and visualize evolutionary conservation in macromolecules. Nucleic Acids Res 2016; 44(W1): W344-350.
27Thompson AA, Peersen OB. Structural basis for proteolysis-dependent activation of the poliovirus RNA-dependent RNA polymerase. EMBO J 2004; 23(17): 3462-3471.
28Ciccozzi M, Lo Presti A, Cella E, Giovanetti M, Lai A, El-Sawaf G, et al. Phylogeny of dengue and chikungunya viruses in Al Hudayda governorate, Yemen. Infect Genet Evol 2014; 27: 395-401.
29Ciccozzi M, Peletto S, Cella E, Giovanetti M, Lai A, Gabanelli E, et al Epidemiological history and phylogeography of West Nile virus lineage 2. Infect Genet Evol 2013; 17: 46-50.
30Stollar V, Thomas VL. An agent in the Aedes aegypti cell line (Peleg) which causes fusion of Aedes albopictus cells. Virology 1975; 64: 367-377.
31Calzolari M, Ze-Ze L, Ruzek D, Vazquez A, Jeffries C, Defilippo F, et al. Detection of mosquito-only flaviviruses in Europe. J Gen Virol 2012; 93: 1215-1225.
32Vazquez 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-229.
33Parreira R, Cook S, Lopes A, de Matos AP, de Almeida AP, Piedade J, et al. Genetic characterization of an insect-specific flavivirus isolated from Culex theileri mosquitoes collected in southern Portugal. Virus Res 2012; 167: 152-161.
34Ergunay K, Gunay F, Sarikaya Y, Kar S, Orsten S, Oter K, et al. Active biosurveillance reveals first records of Chikungunya virus from Turkey: Implications for public health in Turkey and Eastern Europe. 26th European Congress of Clinical Microbiology and Infectious Diseases (ECCMID) Amsterdam, O149; April 9-12, 2016.
35Hoshino 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-414.
36Enserink M. A mosquito goes global. Science 2008; 320: 864-866.
37Scholte EJ, Den Hartog W, Braks M, Reusken C, Dik M, Hessels A. First report of a North American invasive mosquito species Ochlerotatus atropalpus (Coquillett) in the Netherlands, 2009. Euro Surveill 2009; 14: pii19400.
38Lu G, Gong P. A structural view of the RNA-dependent RNA polymerases from the flavivirus genus. Virus Res 2017; 234: 34-43.
39Turkish Embassy in Naypyidaw, Ministry of Economy, EIU, World Bank. Turkey-Myanmar bilateral economic and commercial relations. [Online] Avaliable at: [Cited on 5 April 2019].