Factors associated with mosquito pool positivity and the characterization of the West Nile viruses found within Louisiana during 2007
© Christofferson et al. 2010
Received: 20 April 2010
Accepted: 25 June 2010
Published: 25 June 2010
West Nile virus (WNV) is an arbovirus of public health importance in the genus Flavivirus, a group of positive sense RNA viruses. The NS3 gene has a high level of substitutions and is phylogenetically informative. Likewise, substitutions in the envelope region have been postulated to enable viruses to subvert immune responses. Analysis of these genes among isolates from positive mosquitoes collected in Louisiana illustrates the variation present in the regions and provides improved insight to a phylogenetic model. Employing a GIS eco-regionalization method, we hypothesized that WNV pool positivity was correlated with regional environmental characteristics. Further, we postulated that the phylogenetic delineations would be associated with variations in regional environmental conditions.
Type of regional land cover was a significant effect (p < 0.0001) in the positive pool prediction, indicating that there is an ecological component driving WNV activity. Additionally, month of collection was significant (p < 0.0001); and thus there is a temporal component that contributes to the probability of getting a positive mosquito pool. All virus isolates are of the WNV 2002 lineage. There appears to be some diversity within both forested and wetland areas; and the possibility of a distinct clade in the wetland samples.
The phylogenetic analysis shows that there has been no reversion in Louisiana from the 2002 lineage which replaced the originally introduced strain. Our pool positivity model serves as a basis for future testing, and could direct mosquito control and surveillance efforts. Understanding how land cover and regional ecology effects mosquito pool positivity will greatly help focus mosquito abatement efforts. This would especially help in areas where abatement programs are limited due to either funding or man power. Moreover, understanding how regional environments drive phylogenetic variation will lead to a greater understanding of the interactions between ecology and disease prevalence.
West Nile Virus (WNV) is the most widely distributed arbovirus in the world, occurring on all continents save Antarctica . Its lack of vector specificity compared to other arboviruses has allowed it to use a wide variety of mosquito species in its enzootic cycles . WNV was introduced into the United States in 1999 and from its entry point of New York City it spread across the continental United States. Phylogenetic evidence traced this strain to a similar strain isolated in Israel in 1998 . In 2001 a new WNV strain appeared. In 2002, this genotype became the dominant WN02 strain that was significantly associated with an increase in numbers of human morbidity and mortality cases in the US. In fact, the number of deaths from 1999-2001 were significantly less than the number of deaths in 2002 alone, though whether this association is due to direct virulence in humans or an indirect result of the virulence in birds remains unclear [1, 3, 4].
WNV is a member of the genus Flavivirus, a group of positive sense RNA viruses. The genome is composed of a single open reading frame that produces ten viral proteins: three structural proteins (capsid C, membrane prM/M, envelope E) and seven non-structural proteins (NS1, NS2A, NS2B, NS3 NS4A, NS4B, and NS5) . The NS3 gene plays an important role in the replication of the virus, encoding a protein with four functions: a serine protease, a nucleoside triphosphatase, an RNA 5'triphosphatase, and a helicase . Phylogenetic analyses of WNV have most commonly utilized differences in the envelope protein, but the capsid, prM protein, and non-structural proteins have also been informative [4, 6–12]. Analyses done on complete genomes have given similar results to trees made from prM and envelope proteins . To determine the genetic variability in Louisiana, the envelope coding and NS3 coding regions were analyzed. The NS3 gene has a high level of substitution and is phylogenetically informative  and mutations in the envelope region have been postulated to enable viruses to subvert immune responses . Analysis of these genes would illustrate the variation present in Louisiana as well as provide improved insight for our phylogenetic model.
Many phylogenetic studies have been geographically focused [6, 8, 14, 15]. Geographic Information Systems (GIS) based "region" classifications were successfully used to model WNV transmission risk in humans in northeast Ohio where local environmental features to model transmission . Employing a similar eco-regionalization method, we hypothesized that WNV pool positivity was correlated with regional environmental characteristics. Further, we postulated that the phylogenetic delineations would be associated with variations in regional environmental conditions.
Positive Pool Predictors
Before building the WNV model, we tested several factors via multinomial regression to determine whether these were significant predictors of the mosquito species found, regardless of WNV positivity. Interestingly, no factors (land cover, month) were significant for the 2007 data. Therefore, it was determined that no noteworthy correlation existed between these variables and thus no colinearity issues would arise by inclusion of all variables in the model.
Species in Positive Mosquito Pools
Number of positive pools
Genetic Distances for 2007 Samples
Origin and classification of sequenced isolates
Accession Number (E/N)
Differences of Least Square Means
In some locations, WNV positivity in mosquito pools serves as a predictor for human cases . Our findings that the month of collection and eco-region were significant predictors of positive pools suggest that ecological and temporal factors influence WNV activity and can assist the public health sector predict or prevent cases of human WNV infection. Mosquito abatement programs are operated on a parish-wide basis, so any useful model would ideally work for the parish as a whole. Therefore, looking at the ecology of the parish as a whole- and even grouping parishes into ecological regions or types- will assist in determining a robust model as we did in our pool positivity model, can be used at the very least, on a parish level. Ideally, a model would serve the state as a whole, focusing on regional activity.
Phylogenetic analysis will also help us better understand how the changes in the genome may lead to a change in virulence. It is interesting that the statistical model predicts positive pools in forested areas while the phylogenetic study shows more diversity in the wetland areas. It is also interesting that there is more genetic diversity within Louisiana than without. That is, the diversity when comparing between the wetlands and forests within Louisiana is greater than the diversity when comparing the Louisiana isolates and those outside of the state. This suggests that WNV evolution is not as constrained in Louisiana as compared to other locations as has been suggested by others . Further, some component of eco-regions could have an important role in diversification. This could be due to a difference in bird populations found in each of these ecologies, or it could be a function of extrinsic, microhabitat conditions, such as temperature or humidity which could exert a selection pressure on the viruses.
The phylogenetic analysis shows that there has been no reversion in Louisiana from the 2002 lineage which replaced the originally introduced strain. Our statistical model serves as a basis for future testing, directing mosquito control efforts and surveillance programs. Though we believe our findings to be a significant start to a potentially long term project, it was not without its pitfalls. For example, there is a confounding factor of spatial clustering of the ecologies tested. However, if the spatial component was the only source of diversity, thus discounting the ecological diversity, we would not expect that the phylogenetic signature of the Southern Louisiana strains should be markedly different in topology than the Northern Louisiana strains. The fact that the wetland strains formed a distinct monophyletic group as compared to the forest strains, which were characterized by a lack of phylogenetic structure, suggests the diversity seen within these regions is the result of some other, perhaps ecological, characteristic.
Additionally, the specificity of the land cover classification needs to be more precise. This land cover generalization to parish is a good start to determine if further investigation is garnered, and here we show that it is. Understanding how land cover and regional ecology effects mosquito pool positivity will help focus mosquito abatement efforts. This would especially help in areas where abatement programs are limited due to either funding or man power. Moreover, understanding how regional environments drive phylogenetic variation will lead to a greater understanding of the interactions between ecology and disease prevalence.
According to the mosquito surveillance data, there are two primary species of mosquitoes that serve as possible WNV vectors to people in Louisiana: Culex quinquefasciatus and Aedes albopictus. There were twenty six species from eight genera captured and submitted for testing. Of these, fifteen species were found to be positive for WNV (Table 1).
Mosquito Trap Sites and Field Collection of Mosquitoes
Data was provided by parish mosquito control departments. Each parish operates independently with its own trapping protocols and methods. Not all parishes actively sampled throughout the year due to considerations of the local mosquito activity levels. We therefore analyzed the months that all target parishes had in common (June-November), which captured the majority of the WNV transmission period .
Mosquitoes were sexed and the females were pooled according to genus and species; that is, one pool consisted of a single species. The pools were then homogenized and submitted to the Louisiana Animal Disease Diagnostic Laboratory (LADDL) at the Louisiana State University (LSU) School of Veterinary Medicine (SVM). LADDL is the state testing facility for mosquito pools for all parishes, so criteria for positive pools is the same across all parishes.
Virus Detection and Sequencing
Pools were obtained from the LADDL at the LSU SVM. Viral RNA was extracted from 140 μl of the supernatant from the mosquito pool homogenate using the QIAmp Viral RNA Extraction Kit following manufacturer's instructions (Qiagen, Valencia, CA). One microliter of the extracted viral RNA suspension was used as template for the reverse-trascriptase polymerase chain reaction (RT-PCR) using Superscript™ III RT-PCR kit (Invitrogen, Carlsbad, CA) with the previously described protocol . Upon confirmation of the presence of amplified viral DNA by gel electrophoresis, the remaining sample was cleaned using Qiagen PCR Cleanup kit following manufacturer's instructions. Cleaned DNA was then sent to the Gene Probes and Expression Systems Laboratories of the Division of Biotechnology and Molecular Medicine at the Louisiana State University School of Veterinary Medicine for sequencing. Sequencing was performed on a Beckman Coulter 8800 (Pasadena, CA) using the manufacturer's reagents and methods.
SAS version 9.1.3 was used to code the data as a binary response, where a mosquito pool that was positive for was coded as '1,' while the negative pools were coded as '0.' The probabilities reported are the probabilities of the event = 1 (WNV positive). A confidence level of 95% was used for all tests; a stepwise selection process was invoked to cull out non-significant effects from the model. PROC GLIMMIX with a binary distribution specified was used to construct the model and obtain the odds ratios, as well as to obtain differences in least squares means for the effects of parish and month. PROC GLIMMIX is a useful alternative to PROC LOGISTIC when wanting to obtain differences in least squares and/or modeling random effects, which is not easily done in PROC LOGISTIC. Predicted means and odds ratios are the same between the two procedures. Confidence intervals and tests for significance will vary when random effects are included the model in PROC GLIMMIX, but there were no random effects modelled here .
The complete envelope (E) and non-structural protein 3 (NS3) genes were aligned separately for all successfully recovered isolates as well as 24 reference strains, which represented the Israel 1998 strain and isolates spanning the contiguous United States and one Mexican isolate . Alignments and the creation of the E-NS3 contigs were done using Vector NTI software and exported to GeneDoc for trimming. Alignments were imported into MEGA 4 and converted to MEGA format. Bootstrap analysis (n = 1000) using maximum parsimony was performed and a tree produced. As the topologies of the E and NS3 phylogenies were the same, our tree represents a contiguous sequence of these two genes. The topology tree was collapsed with a node confidence of 70%. Genetic distances and means were also obtained using Mega 4 . Not all samples in the statistical model were available for the phylogenetic analysis.
We would like to thank the mosquito abatement programs for East Baton Rouge, St. Tammany, Iberville, Ouachita, and Caddo parishes. Thanks also should be extended to Dan Chisenhall for logistical support.
Appendix: Accession Numbers
The GenBank Accession numbers of the sequences reported in this paper are:
[GenBank:HM538823, GenBank:HM538811, GenBank:HM538824, GenBank:HM538812, GenBank:HM538829, GenBank:HM538807, GenBank:HM538825, GenBank:HM538810, GenBank:HM538827, GenBank:HM538806, GenBank:HM538826, GenBank:HM538808, GenBank:HM538828, GenBank:HM538809, GenBank:HM538821, GenBank:HM538813, GenBank:HM538817, GenBank:HM538805, GenBank:HM538822, GenBank:HM538816, GenBank:HM538818, GenBank:HM538803, GenBank:HM538819, GenBank:HM538804, GenBank:HM538820, GenBank:HM538814.]
- Kramer LD, Styer LM, Ebel GD: A global perspective on the epidemiology of West Nile virus. Annu Rev Entomol 2008, 53:61–81.PubMedView Article
- Lanciotti RS, Roehrig JT, Deubel V, Smith J, Parker M, Steele K, Crise B, Volpe KE, Crabtree MB, Scherret JH, et al.: Origin of the West Nile virus responsible for an outbreak of encephalitis in the northeastern United States. Science 1999, 286:2333–2337.PubMedView Article
- Moudy RM, Meola MA, Morin LL, Ebel GD, Kramer LD: A newly emergent genotype of West Nile virus is transmitted earlier and more efficiently by Culex mosquitoes. Am J Trop Med Hyg 2007, 77:365–370.PubMed
- Davis CT, Ebel GD, Lanciotti RS, Brault AC, Guzman H, Siirin M, Lambert A, Parsons RE, Beasley DW, Novak RJ, et al.: Phylogenetic analysis of North American West Nile virus isolates, 2001–2004: evidence for the emergence of a dominant genotype. Virology 2005, 342:252–265.PubMedView Article
- Brinton MA: The molecular biology of West Nile Virus: a new invader of the western hemisphere. Annu Rev Microbiol 2002, 56:371–402.PubMedView Article
- Bertolotti L, Kitron UD, Walker ED, Ruiz MO, Brawn JD, Loss SR, Hamer GL, Goldberg TL: Fine-scale genetic variation and evolution of West Nile Virus in a transmission "hot spot" in suburban Chicago, USA. Virology 2008, 374:381–389.PubMedView Article
- Jerzak G, Bernard KA, Kramer LD, Ebel GD: Genetic variation in West Nile virus from naturally infected mosquitoes and birds suggests quasispecies structure and strong purifying selection. J Gen Virol 2005, 86:2175–2183.PubMedView Article
- Ebel GD, Dupuis AP, Ngo K, Nicholas D, Kauffman E, Jones SA, Young D, Maffei J, Shi PY, Bernard K, Kramer LD: Partial genetic characterization of West Nile virus strains, New York State, 2000. Emerg Infect Dis 2001, 7:650–653.PubMedView Article
- Davis CT, Galbraith SE, Zhang S, Whiteman MC, Li L, Kinney RM, Barrett AD: A combination of naturally occurring mutations in North American West Nile virus nonstructural protein genes and in the 3' untranslated region alters virus phenotype. J Virol 2007, 81:6111–6116.PubMedView Article
- Beasley DW, Davis CT, Guzman H, Vanlandingham DL, Travassos da Rosa AP, Parsons RE, Higgs S, Tesh RB, Barrett AD: Limited evolution of West Nile virus has occurred during its southwesterly spread in the United States. Virology 2003, 309:190–195.PubMedView Article
- Gray RR, Veras NM, Santos LA, Salemi M: Evolutionary characterization of the West Nile Virus complete genome. Mol Phylogenet Evol 2010,56(1):195–200.PubMedView Article
- Lanciotti RS, Ebel GD, Deubel V, Kerst AJ, Murri S, Meyer R, Bowen M, McKinney N, Morrill WE, Crabtree MB, et al.: Complete genome sequences and phylogenetic analysis of West Nile virus strains isolated from the United States, Europe, and the Middle East. Virology 2002, 298:96–105.PubMedView Article
- Pfeffer M, Foster JE, Edwards EA, Brown MB, Komar N, Brown CR: Phylogenetic analysis of Buggy Creek virus: evidence for multiple clades in the Western Great Plains, United States of America. Appl Environ Microbiol 2006, 72:6886–6893.PubMedView Article
- Andreadis TG, Anderson JF, Vossbrinck CR, Main AJ: Epidemiology of West Nile virus in Connecticut: a five-year analysis of mosquito data 1999–2003. Vector Borne Zoonotic Dis 2004, 4:360–378.PubMedView Article
- Ebel GD, Carricaburu J, Young D, Bernard KA, Kramer LD: Genetic and phenotypic variation of West Nile virus in New York, 2000–2003. Am J Trop Med Hyg 2004, 71:493–500.PubMed
- LaBeaud AD, Gorman AM, Koonce J, Kippes C, McLeod J, Lynch J, Gallagher T, King CH, Mandalakas AM: Rapid GIS-based profiling of West Nile virus transmission: defining environmental factors associated with an urban-suburban outbreak in Northeast Ohio, USA. Geospat Health 2008, 2:215–225.PubMed
- Harper RA: Forest Inventory & Analysis Factsheet Louisiana 2005. Forest Inventory & Analysis Factsheet Louisiana 2005 2006.
- Gleiser RM, Mackay AJ, Roy A, Yates MM, Vaeth RH, Faget GM, Folsom AE, Augustine WF Jr, Wells RA, Perich MJ: West Nile virus surveillance in East Baton Rouge Parish, Louisiana. J Am Mosq Control Assoc 2007, 23:29–36.PubMedView Article
- Mackay AJ, Roy A, Yates MM, Foil LD: West Nile virus detection in mosquitoes in East Baton Rouge Parish, Louisiana, from November 2002 to October 2004. J Am Mosq Control Assoc 2008, 24:28–35.PubMedView Article
- Chisenhall DM, Mores CN: Diversification of West Nile virus in a subtropical region. Virol J 2009, 6:106.PubMedView Article
- Richards SL, Mores CN, Lord CC, Tabachnick WJ: Impact of extrinsic incubation temperature and virus exposure on vector competence of Culex pipiens quinquefasciatus Say (Diptera: Culicidae) for West Nile virus. Vector Borne Zoonotic Dis 2007, 7:629–636.PubMedView Article
- Chisenhall DM, Vitek CJ, Richards SL, Mores CN: A method to increase efficiency in testing pooled field-collected mosquitoes. J Am Mosq Control Assoc 2008, 24:311–314.PubMedView Article
- Christofferson RC: Old vs. New: A comparison of PROC LOGISTIC and PROC GLIMMIX. SAS Global Forum; San Antonio, TX; 2008.
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 24:1596–1599.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.