Skip to main content

Evolutionary relationships of West Nile virus detected in mosquitoes from a migratory bird zone of Colombian Caribbean

Abstract

Background

West Nile virus (WNV) is a member of the genus Flavivirus, and it is transmitted between Culex sp. mosquitoes and avian hosts. Equids and humans are commonly infected with WNV as dead-end hosts, and the signs and symptoms of infection range from mild illness to neurologic symptoms as encephalitis, meningitis and sometimes death. Previous phylogenetic studies have classified WNV into six genetically distinct lineages and provided valuable insight on WNV dispersal patterns within the Americas and its emergence in different geographic areas. In this study, we isolated, sequenced and genetically characterized the NS5 and envelope genes for two WNV strains detected from Northern of Colombia. Herein we describe the evolutionary relationships with representative WNV-strains isolated in a variety of epidemic outbreaks and countries, to define the phylogeographic origin and possible implications in the epidemiology of this emergent virus in Colombia.

Findings

Fragments of the NS5 and Envelope genes were amplified with RT-PCR and sequenced to obtain 1186-nt and 1504-nt portions, respectively. Our sequences were aligned with 46 sequences from WNV-strains collected in the U.S., Mexico and Argentina for phylogenetic reconstruction using Bayesian methods. Sequence analyses identified unique non-synonymous substitutions in the envelope gene of the WNV strains we detected, and our sequences clustered together with those from the attenuated Texas – 2002 genotype.

Conclusions

A new strain closely related to attenuated strains collected in Texas during 2002 was identified from Colombia by phylogenetic analysis. This finding may explain the absence of human/equine cases of WNV-encephalitis or severe disease in Colombia and possibly other regions of South America. Follow-up studies are needed in ecosystems used by migratory birds areas and virological/entomological surveillance.

Findings

Introduction

West Nile virus (WNV) is a member of the Japanese encephalitis antigenic complex (JEV) within the Flavivirus genus and is transmitted by Culex spp. mosquitoes among birds. Other vertebrates such as mammals and reptiles also become infected [1]. Human infection with WNV causes mild to severe illness, sometimes affecting the nervous systems and provoking encephalitis, meningitis and death [2].

Since the first reports of WNV isolation from Africa, Europe, India, Russia, Israel, France, and its 1999 introduction into North America, WNV has extended its geographic distribution throughout the United States [3]. Subsequently, in the following years this arbovirus was detected in Canada, Mexico, Guatemala, Caribbean islands, and South America [46]. Serologic evidence for the natural circulation of WNV in Colombia has been observed in equids sampled from the department of Córdoba [7, 8] and other regions from the Caribbean [911]. Despite WNV being isolated for the first time from captive flamingoes in Santa Fé Zoo (Medellín, Colombia) [12], it is not clear why WNV has not been isolated from or been the cause of detectable disease in horses or humans in Colombia. Possible explanations include: circulation of WNV in remote enzootic cycles away from human settlements, limited vector competence of mosquito species, ornitophilic blood-feeding preferences, cross-protective immunity in humans from other flaviviruses (dengue, Saint Louis encephalitis viruses), or the circulation of WNV-strains with low or attenuated virulence [1012].

WNV occurs in four major lineages, but lineage 1 is epidemiologically relevant. Lineage 1 is subdivided into three clades (1a, 1b, 1c); clade 1a contains isolates from Africa, Europe, the Middle East, Russia and Americas [1, 3]. WNV-American strains have close relationships with three Old World isolates: IS98-ST1 (Israel - 1998), PaH001 (Tunez-1997) and goose-03 (Hungary-2003), and extensive studies have allowed detailed investigations of WNV microevolution in different areas over time and also the emergence of new genotypes [3]. In this sense, phylogenetic analysis has enabled the understanding of epidemiological patterns of emergence, dispersal routes, adaptation to new hosts/mosquitoes species, and spatio-temporal patterns of evolution [1, 3, 5, 12]. Evolutionary studies are necessary to identify “drivers” of emergence, molecular evolution of virulence and associations to ecological factors that allow the establishment of this arbovirus pathogen in human populations.

Between 2011 and 2013, a surveillance study was performed at one locality that is characteristic of a large migratory bird population in northern Colombia. We detected WNV in pools of mosquitoes and amplified two viral regions that were used to establish the phylogenetic relationships with strains isolated in the U.S., Mexico and Argentina, with the goal of defining the evolutionary relationships with genotypes previously isolated from outbreaks in the geographic areas mentioned above.

Material and methods

Samples

During virological surveillance for detection of emerging and re-emerging arboviruses between 2011–2013 in San Bernardo del Viento (Córdoba, Colombia) (9° 21´ 30.97” N, 75° 58´ 37.28” W) (Fig. 1), mosquitoes were collected using CDC-light/EVS traps that were baited with dry ice (CO2). All the insects sampled were separated into pools through morphological identification and triturated using minimum essential medium (MEM) supplemented with 10 % fetal bovine serum, 1 % penicillin, then clarified by centrifugation at 13,000 rpm for 30 min. Supernatant was used for RNA extraction and generic nested reverse transcriptase polymerase chain reaction (RT-PCR) for detection of flaviviruses (Table 1) [13]. Two mosquito pools, each containing 40 and 50 specimens (whose mosquitoes had a high percentage −99.57 % of similarity in their DNA Barcode sequences with Culex (Melanoconion) erraticus) were positive for a member of the Flavivirus genus. Sequencing of PCR-products, BLASTN and Neighbor-Joining phylogenetic analysis allowed identification of these sequences as West Nile virus.

Fig. 1
figure1

Map of the study area showing San Bernardo del Viento (Córdoba, Colombia)

Table 1 Primers and RT-PCR (generic/nested) conditions for arboviral detection in target groups

Molecular protocols

Total RNA was extracted from the PCR-positive homogenates with the RNeasy kit (Qiagen, Valencia) and used for RT-PCR using the One-Step RT-PCR Kit (Qiagen, Valencia, CA) as described previously [14, 15]. Two viral genes were amplified: NS5 (primers: FU1PMF-TACAACATGATGGGVAARAGWGARAA/cFD3PMR-ARCATGTCTTCYGTBGTCATCCA) and Envelope (primers: WN1101-GATGAATATGGAGGCGGTCA/WN1816A-CCGACGTCAACTTGACAGTG and WN1751-TGCATCAAGCTTTGGCTGGA/WN2504A TCTTGCCGGCTGATGTCTAT).

Sequence analysis

A total of 46 genome sequences of envelope/NS5 from WNV isolates, representative of American clade 1a-Cluster 4 geographic locations [16], were downloaded from GenBank (Table 2). Nucleotide sequences were aligned using the MAFFTv7.0 (http://mafft.cbrc.jp/alignment/server/) and then transferred in FASTA format to BioEdit (http://www.mbio.ncsu.edu/BioEdit/bioedit.html) for manual editing, keeping gaps consistent within the reading frame. The sequences of NS5 and envelope of the WNV-strains were aligned and evaluated in jModelTest v2.1.4 [17] using the Akaike criterion informative for identify appropriate substitution model nucleotide. One test was also performed for concatenated sequences. The XML file for Bayesian analysis was created in BEAUti v1.5.4 (http://www.molecularevolution.org/software/phylogenetics/beauti), describing model of sequences, invariants, gamma distribution, size of the chain run (20 million of generations), coalescent constant population and for accommodate the variation in substitution rate among branches, a random local clock model was chosen for this analysis [18].

Table 2 Sequences downloaded of Genbank and background information of WNV strain/isolates used in this study

Bayesian phylogenetic analysis was performed using the BEAST software package v2.1.3 [19], and estimation of the maximum clade credibility (MCC) phylogenetic tree was achieved using TreeAnnotator-v2.0.2. BEAST output was viewed with TRACERv-1.5 and evolutionary trees were generated in the FigTree-v1.3.1.

DNAsp-v5.0 [20] was used to establish polymorphic sites between NS5/envelope sequences characterized in our study and reference sequences of representative WNV strains.

Results and discussion

The two pools infected with WNV, corresponding probably to the mosquito species Culex (Melanoconion) erraticus (codes: CDCCA4-12 – CDCCZ2-21), were collected in September, 2012 from coastal mangroves in a migratory bird zone located in San Bernardo del Viento (Córdoba Department). Both isolates were RT-PCR amplified, including one sequence of 1504 nt in length from the envelope gene (Genbank accession number: KM212943 – KM212944) and another 1186 nt long from the NS5 gene (Genbank accession number: KM212941 – KM212942). jModelTest-v2.1.4 estimated the same model of substitution nucleotide, General-Time-Reversible (GTR) + gamma distribution (−lnL = 4792.34, AICc = 9798.9686) for both viral regions. A concatenated file with two NS5/envelope sequences (2690 nt) estimated the same model and was used for phylogenetic inference.

The phylogenetic tree inferred with Bayesian methods indicated that our two WNV-samples were closely related to strains Mosquito-v4369, Bird1519 and v4095, all belonging to southeastern coastal Texas genotype (Fig. 2). Our sequences differed from that of this WNV-genotype in three positions in the envelope gene: only two positions are unique (shaded positions: 547,641) and position 555 shared the same nucleotide with strains HNY1999, COL524 and COL9835. These substitutions were all synonymous. No unique mutations were found in NS5, but our sequences shared similar differences to the WNV-Texas genotype and other WNV genotypes (Table 3). The southeastern coastal Texas genotype includes several isolates collected in 2002 from Texas, and is considered to have an attenuated phenotype with a small plaques (sp) size, temperature sensitivity (ts), reduced replication in cell culture and reduced neuroinvasiveness that is dose-dependent [21]. The southeastern coastal Texas genotype has not been detected since 2002, suggesting its possible extinction [3]. In Colombia, previous work characterizing WNV-strains isolated from flamingoes (COL524/COL9835) showed a close genetic relationship with WNV strains isolated in Louisiana in 2001 and the NY99 strain (Fig. 2), but in vitro phenotypic characterization showed differences with the attenuated Texas-genotype. In fact, COL524/COL9835 has high virulence in chicken eggs and newborn/4-week-old Balb/c mice [12].

Fig. 2
figure2

Phylogenetic tree estimated by Bayesian analysis of 46 sequence strains (NS5-Envelope) of West Nile virus data under the GTR + G model of nucleotide substitution. The number accession of Genbank is followed by abbreviated name of WNV-strain. The branches in red belong to characterized previously by Osorio et al. (2012) and characterized by ours

Table 3 Polymorphic sites between envelope and NS5 sequences belonging to southeastern coastal Texas genotype and our samples from Northern of Colombia (Córdoba, San Bernardo del Viento)

Our results demonstrate genetic diversity of WNV strains circulating in Colombia. The presence of the attenuated coastal Texas genotype could explain, in part, the lack of human and equine cases detected. Previous work has shown WNV-seropositivity in horses from the Caribbean region with no disease reported [711]. Different WNV genotypes may converge along migratory bird flyways, which pass from the U.S. to Mexico, Colombia and Venezuela, Caribbean Islands, and other parts of South America south to Argentina. The high diversity of migratory birds in certain areas of Caribbean Colombia [3, 4, 12, 22] may be important for WNV maintenance.

Finally, our results indicate a close evolutionary relationship with the attenuated coastal Texas genotype requires further studies in cell culture and animal models to confirm the attenuated phenotype. Additional surveillance focused on avian and mosquito fauna is also needed to obtain more isolates of WNV in conserved Colombian ecosystems to further examine the genetic diversity of WNV and possible strain dissemination to another geographic areas.

References

  1. 1.

    Kramer L, Styer L, Ebel G. A Global Perspective on the epidemiology of West Nile virus. Annu Rev Entomol. 2008;53:61–81.

    Article  CAS  PubMed  Google Scholar 

  2. 2.

    Samuel M, Diamond M. Pathogenesis of West Nile Virus Infection: a Balance between Virulence, Innate and Adaptive Immunity, and Viral Evasion. J virol. 2006;80:9349–60.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. 3.

    Mann B, McMullen A, Swetnam D, Barret A. Molecular Epidemiology and Evolution of West Nile Virus in North America. Int J Environ Res Public Health. 2013;10:5111–29.

    Article  PubMed Central  PubMed  Google Scholar 

  4. 4.

    Komar N, Clarck G. West Nile virus activity in Latin America and Caribbean. Rev Panam Salud Publica. 2006;19:112–7.

    Article  PubMed  Google Scholar 

  5. 5.

    Morales MA, Barrandeguy M, Fabbri C, Garcia JB, Vissani A, Trono K, et al. West Nile virus isolation from equines in Argentina, 2006. Emerg Infect Dis. 2006;10:1559–61.

    Article  Google Scholar 

  6. 6.

    Ometto T, Durigon EL, de Araujo J, Aprelon R, de Aguiar DM, Cavalcante GT, et al. West Nile virus surveillance, Brazil, 2008–2010. Trans R Soc Trop Med Hyg. 2013;107:723–30.

    Article  PubMed  Google Scholar 

  7. 7.

    Mattar S, Edwards E, Laguado J, González M, Alvarez J, Komar N. West Nile virus infection in Colombian horses. Emerg Infect Dis. 2005;11:1497–8.

    Article  PubMed Central  PubMed  Google Scholar 

  8. 8.

    Jaramillo M, Peña J, Berrocal L, Komar N, González M, Ponce C, et al. Vigilancia centinela para el Virus del Oeste del Nilo en Culícidos y Aves Domésticas en el Departamento de Córdoba. Rev MVZ Cordoba. 2005;10:633–8.

    Google Scholar 

  9. 9.

    Mattar S, Parra M, Torres J. Limitaciones para el serodiagnóstico del virus del Nilo del oeste en zonas endémicas con co-circulación de Flavivirus en el Caribe colombiano. Colombia Médica. 2005;36:179–85.

    Google Scholar 

  10. 10.

    Góes-Rivillas Y, Taborda N, Díaz F, Góngora A, Rodas JD, Ruiz-Sáenz J, et al. Antibodies to West Nile virus in equines of Antioquia and Meta, Colombia, 2005–2008. Rev Colomb Cienc Pecu. 2010;23:462–70.

    Google Scholar 

  11. 11.

    Mattar S, Komar N, Young G, Alvarez J, Gonzalez M. Seroconversion for West Nile and St. Louis encephalitis viruses among sentinel horses in Colombia. Mem Inst Oswaldo Cruz. 2011;106:976–9.

    Article  PubMed  Google Scholar 

  12. 12.

    Osorio J, Ciuoderis K, Lopera J, Piedrahita L, Murphy D, Levasseur J, et al. Characterization of West Nile Viruses isolated from captive American flamingoes (Phoenicopterous ruber) in Medellin, Colombia. Am J Trop Med Hyg. 2012;87:565–72.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. 13.

    Sánchez-Seco M, Rosario D, Domingo C, Hernández L, Valdés K, Gúzman M, et al. Generic RT-nested-PCR for detection of flaviviruses using degenerated primers and internal control followed by sequencing for specific identification. J Virol Methods. 2005;1:101–9.

    Article  Google Scholar 

  14. 14.

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

    PubMed Central  CAS  PubMed  Google Scholar 

  15. 15.

    Beasley DWC, Davis CT, Guzman H, Vanlandingham DL, Travassos da Rosa APA, Parsons RE. Limited evolution of West Nile virus has occurred during its southwesterly spread in the United States. Virology. 2003;309:190–5.

    Article  CAS  PubMed  Google Scholar 

  16. 16.

    May F, Davis T, Tesh R, Barrett A. Phylogeography of West Nile Virus: from the Cradle of Evolution in Africa to Eurasia, Australia, and the Americas. J Virol. 2011;85:2964–74.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. 17.

    Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772.

    Article  CAS  PubMed  Google Scholar 

  18. 18.

    Drummond A, Suchard M. Bayesian random local clocks, or one rate to rule them all. BMC Biol. 2010;8:114.

    Article  PubMed Central  PubMed  Google Scholar 

  19. 19.

    Bouckaert R, Heled J, Kühnert D, Vaughan TG, Wu C-H, Xie D, et al. BEAST2: A software platform for Bayesian evolutionary analysis. PLOS Comput Biol. 2014;10, e1003537.

    Article  PubMed Central  PubMed  Google Scholar 

  20. 20.

    Librado P, Rozas J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.

    Article  CAS  PubMed  Google Scholar 

  21. 21.

    Davis T, Beasley D, Guzman H, Siirin M, Parsons R, Tesh R, et al. Emergence of attenuated West Nile virus variants in Texas, 2003. Virology. 2004;320:342–50.

    Article  Google Scholar 

  22. 22.

    Berrocal L, Peña J, González M, Mattar S. West Nile virus; ecology and epidemiology of an emerging pathogen in Colombia. Rev salud pública. 2006;8:218–28.

    Article  Google Scholar 

Download references

Acknowledgements

RHL was recipient of a doctoral fellowship from “Programa de Doctorados Nacionales – Colciencias” (Convocatoria 528). The authors acknowledge to “Departamento Administrativo de Ciencia, Tecnología e Innovación – Colciencias” by supporting research with grant 111549326198. Finally, JCGG was recipient of a Full-Time Professor Program (“Dedicación Exclusiva”) 2012–3 of the Medicine Faculty at Universidad de Antioquia.

To the entomologist Juan David Suaza who helped in all field work and entomological processing of the samples. Also to the direction and personnel of Laboratory of Biology and Insect Systematics, Sciences Faculty, Universidad Nacional de Colombia, sede Medellin, where all the molecular work was done. To Drs. Scott Weaver and Albert Auguste from National Galveston Laboratory - Medical Branch of Texas University (UTMB).

Author information

Affiliations

Authors

Corresponding authors

Correspondence to Richard Hoyos López or Juan Carlos Gallego-Gómez.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RHL performed the molecular experiments, phylogenetic analysis and wrote the first draft of the paper. JCGG conceived the study, SIU help with design and develop study. SIU and JCGG critically reviewed this manuscript. All authors read and approved the final manuscript.

Rights and permissions

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

López, R.H., Soto, S.U. & Gallego-Gómez, J.C. Evolutionary relationships of West Nile virus detected in mosquitoes from a migratory bird zone of Colombian Caribbean. Virol J 12, 80 (2015). https://doi.org/10.1186/s12985-015-0310-8

Download citation

Keywords

  • Reverse Transcriptase Polymerase Chain Reaction
  • West Nile Virus
  • Migratory Bird
  • Japanese Encephalitis
  • Envelope Gene