|Year : 2019 | Volume
| 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
|Date of Submission||08-Dec-2018|
|Date of Decision||20-Mar-2019|
|Date of Acceptance||02-Apr-2019|
|Date of Web Publication||30-May-2019|
Unit of Medical Statistics and Molecular Epidemiology, University Campus Bio-Medico of Rome, Rome
Source of Support: None, Conflict of Interest: None
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.
Keywords: Culex virus, Flaviviridae, Phylogenetic analysis, Evolution
|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-23
|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 2020 Jun 2];12:216-23. Available from: http://www.apjtm.org/text.asp?2019/12/5/216/259242
Eleonora Cella, Domenico Benvenuto. These authors contributed equally to this work.
| 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.
Insect-specific flaviviruses (ISFs) are distinct for the inability to be propagated in vertebrate cell lines and lack of vertebrate hosts. ISFs have been recognized as a major phylogenetic group along with mosquito-borne, tick-borne and no-known-vector flaviviruses,,, 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. 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,. 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,,. In particular, Culex (Cx.) theileri flavivirus (CTFV) seems to affect the flight activity of Culex mosquitos.
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 (http://www.ncbi.nlm.nih.gov/ genbank/). All the sequences were aligned and manually edited using the Bioedit program v7.2.5, as already described.
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. 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.
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). The Maximum Likelihood tree was reconstructed using MEGA7 with Tamura-Nei model as evolutionary model selected by ModelTest,,.
Adaptive Evolution Server (http://www.datamonkey.org/) 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. 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 http://www.mfa.gov.tr/turkey_myanmar-bilateraleconomic-and-commercial-relations.en.mfa. 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 and the programs ClustalO or Muscle,. Structural templates homologous to the sequences studied within this work have been identified with HHpred and SwissModel template search tools. Homology models were built with the program Modeller and validated with QMEAN and Prosa II  software tools. Three-dimensional structures were displayed and studied with PyMOL (The PyMOL Molecular Graphics System, Version 2.0 Schroedinger, LLC) or Chimera. Sequence conservation and mapping analysis has been carried out with the Consurf web tool.
| 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: Likelihood mapping of Culex theileri flavivirus dataset. The dots inside the triangles represents the likelihood of the possible unrooted topologies for each quartet. Numbers indicate the percentage of dots in the centre of the triangle corresponding to phylogenetic noise (star-like trees).|
Click here to view
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: Maximum likelihood tree of Culex theileri flavivirus. Along the branches an asterisk (*) indicating a statistical value from bootstrap (>70%).|
Click here to view
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. 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: Multiple sequence alignment between the fragments of Culex theileri flavirus RdRp domain of NS5 proteins and structurally solved homologous NS5 proteins from flaviviruses PDB codes of the RdRp domains are: 2HFZ (WNV); 4HDG (JEV); 4K6M (JEV-SA-14); 5JJS (DENV-3); 5WZ3 (ZIKV). Blue and green bars below the alignment indicate deletions and insertions in the CTFV RdRp sequences. Cyan arrows mark the CTFV RdRp sites under selective pressure. Labelled grey bars mark the conserved viral RdRp catalytic motifs. The red line shows the position of the “Middle finger” subdomain. Sequence position numbering refers to 4HDG structure.|
Click here to view
|Figure 4: Cartoon model of the chain A structure of the full-length Japanese encephalitis virus (PDB code 4K6M). The Methyltransferase (MTase) and RNA dependent RNA polymerase (RdRp) domains are labelled. Grey spheres indicate zinc ions, while sphere model in the MTase domain represents S-adenosyl-L-homocysteine. Dark red ribbon marks the CTFV RdRp fragments coverage of the full-length structure. Positions of the N-and|
Click here to view
|Figure 5: Structural superposition between the homology models of HE997074 (red cartoon) and KU958176 (cyan cartoon). Ovals enclose loops displaying conformational differences between the two models. Oval 1 corresponds to the “Middle” subdomain containing the position 567 [Figure 3] under selective pressure and the insertion Asn-His (positions 584585 in Figure 3). Oval 2 encloses the region of the large insertion in 4HDG (positions 636-644, Figure 3).|
Click here to view
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: Superposition between the homology model of CTFV RdRp NS5 with code HE997074 (red cartoon) and the structural template 4HDG (orange cartoon). Active site GTP is displayed as stick model. Grey and green spheres denote zinc and magnesium ions. Cyan and green sticks models indicate CTFV RdRp positions under selective pressure and insertions with respect to 4HDG. Orange sticks indicate the 4HDG insertion residues. Numbering refers to full-length 4HDG sequence. The two-residues insertion of CTFV RdRp between the 4HDG positions 557-558 and 584-585, respectively are indicated. Left figure is rotated by 180° along the vertical y axis with respect to that represented on the right panel.|
Click here to view
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: Cartoon depiction of the 4HDG structure coloured according to the conservation score calculated by ConSurf web tool. The colour scale ranges from cyan (high sequence variability) to dark red (high sequence conservation). White is the “neutrality” colour. GTP is displayed with stick models. Positions subject to selective pressure are represented as stick models and labelled accordingly. All of them occur in highly variable positions of the RdRp domain of Flavivirus NS5 proteins.|
Click here to view
| 4. Discussion|| |
Flavivirus have a great impact on human health,. 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.
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. 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),.
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,,.
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, several phylogenetically distinct ISFs have been repeatedly detected in various Culex and Aedes spp. Mosquitoes in different areas, suggesting independent introductions and multiple host-switching events. 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,.
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. 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.
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.
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
| References|| |
Mukhopadhyay S, Kuhn RJ, Rossmann MG. A structural perspective of the flavivirus life cycle. Nat Rev Microbiol
Calzolari 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
Moureau 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
Cook S, Holmes EC. Amultigene analysis of the phylogenetic relationships among the flaviviruses (Family: Flaviviridae) and the evolution of vector transmission. Arch Virol
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
|6.|Simmonds 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.
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
Crochu 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
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
Newman CM, Anderson TK, Goldberg TL. Decreased flight activity in Culex pipiens
(Diptera: Culicidae) naturally infected with Culex flavivirus. J Med Entomol
 Hall TA. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser
Schmidt HA, Petzold E, Vingron M, von Haeseler A. Molecular phylogenetics: Parallelized parameter estimation and quartet puzzling. J Parallel Distrib Comput
Swofford 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.
Posada D, Buckley TR. Model selection and model averaging in phylogenetics: advantages of akaike information criterion and Bayesian approaches over likelihood ratio tests. Syst Biol
Kumar S, Stecher G, Tamura K. MEGA7: Molecular evolutionary genetics analysis Version 7.0 for Bigger Datasets. Mol Biol Evol
Pond SL, Frost SD, Muse SV. HyPhy: Hypothesis testing using phylogenies. Bioinformatics
Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ. Jalview version 2: A multiple sequence alignment editor and analysis workbench. Bioinformatics
Sievers F, Higgins DG. Clustal omega for making accurate alignments of many protein sequences. Protein Sci
Edgar RC. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics
Zimmermann 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
Waterhouse 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
Eswar 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.
Benkert P, Kunzli M, Schwede T. QMEAN server for protein model quality estimation. Nucleic Acids Res
(Web Server issue): W510-514.
Wiederstein M, Sippl MJ. ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res
(Web Server issue): W407-410.
Meng EC, Pettersen EF, Couch GS, Huang CC, Ferrin TE. Tools for integrated sequence-structure analysis with UCSF Chimera. BMC Bioinformatics
Ashkenazy 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
Thompson AA, Peersen OB. Structural basis for proteolysis-dependent activation of the poliovirus RNA-dependent RNA polymerase. EMBO J
Ciccozzi 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
Ciccozzi 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
Stollar V, Thomas VL. An agent in the Aedes aegypti
cell line (Peleg) which causes fusion of Aedes albopictus
Calzolari M, Ze-Ze L, Ruzek D, Vazquez A, Jeffries C, Defilippo F, et al. Detection of mosquito-only flaviviruses in Europe. J Gen Virol
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
Parreira 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
Ergunay 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.
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
Enserink M. A mosquito goes global. Science
2008; 320: 864-866.
Scholte 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
Lu G, Gong P. A structural view of the RNA-dependent RNA polymerases from the flavivirus genus. Virus Res
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5], [Figure 6], [Figure 7]