Phylogenetic analysis of West Nile Virus in Maricopa County, Arizona: Evidence for dynamic behavior of strains in two major lineages in the American Southwest

Temperature and time of host-seeking activity impact the efficacy of chemical control interventions targeting the West Nile virus vector, Culex tarsalis

West Nile virus (WNV) is the leading mosquito-borne disease causing-pathogen in the United States. Concerningly, there are no prophylactics or drug treatments for WNV and public health programs rely heavily on vector control efforts to lessen disease incidence. Insecticides can be effective in reducing vector numbers if implemented strategically, but can diminish in efficacy and promote insecticide resistance otherwise. Vector control programs which employ mass-fogging applications of insecticides, often conduct these methods during the late-night hours, when diel temperatures are coldest, and without a-priori knowledge on daily mosquito activity patterns. This study's aims were to 1) quantify the effect of temperature on the toxicity of two conventional insecticides used in fogging applications (malathion and deltamethrin) to Culex tarsalis, an important WNV vector, and 2) quantify the time of host-seeking of Cx. tarsalis and other local mosquito species in Maricopa County, Arizona. The temperature-toxicity relationship of insecticides was assessed using the WHO tube bioassay, and adult Cx. tarsalis, collected as larvae, were exposed to three different insecticide doses at three temperature regimes (15, 25, and 35°C; 80% RH). Time of host-seeking was assessed using collection bottle rotators with encephalitis vector survey traps baited with dry ice, first at 3h intervals during a full day, followed by 1h intervals during the night-time. Malathion became less toxic at cooler temperatures at all doses, while deltamethrin was less toxic at cooler temperatures at the low dose. Regarding time of host-seeking, Cx. tarsalis, Aedes vexans, and Culex quinquefasciatus were the most abundant vectors captured. During the 3-hour interval surveillance over a full day, Cx. tarsalis were most-active during post-midnight biting (00:00-06:00), accounting for 69.0% of all Cx. tarsalis, while pre-midnight biting (18:00-24:00) accounted for 30.0% of Cx. tarsalis. During the 1-hour interval surveillance overnight, Cx. tarsalis were most-active during pre-midnight hours (18:00-24:00), accounting for 50.2% of Cx. tarsalis captures, while post-midnight biting (00:00-06:00) accounted for 49.8% of Cx. tarsalis. Our results suggest that programs employing large-scale applications of insecticidal fogging should consider temperature-toxicity relationships coupled with time of host-seeking data to maximize the efficacy of vector control interventions in reducing mosquito-borne disease burden.

Correlates of Co-Infection with Coccidiosis and Avian Malaria in House Finches (Haemorhous mexicanus)

Pathogens have traditionally been studied in isolation within host systems; yet in natural settings they frequently coexist. This raises questions about the dynamics of co-infections and how host life-history traits might predict co-infection versus single infection. To address these questions, we investigated the presence of two parasites, a gut parasite (Isospora coccidians) and a blood parasite (Plasmodium spp.), in House Finches (Haemorhous mexicanus), a common passerine bird in North America. We then correlated these parasitic infections with various health and condition metrics, including hematological parameters, plasma carotenoids, lipid-soluble vitamins, blood glucose concentration, body condition, and prior disease history. Our study, based on 48 birds captured in Tempe, Arizona, US, in October 2021, revealed that co-infected birds exhibited elevated circulating lutein levels and a higher heterophil:lymphocyte ratio (H/L ratio) compared to those solely infected with coccidia Isospora spp. This suggests that co-infected birds experience heightened stress and may use lutein to bolster immunity against both pathogens, and that there are potentially toxic effects of lutein in co-infected birds compared to those infected solely with coccidia Isospora sp. Our findings underscore the synergistic impact of coparasitism, emphasizing the need for more co-infection studies to enhance our understanding of disease dynamics in nature, as well as its implications for wildlife health and conservation efforts.

The Mosquito Fauna of Arizona: Species Composition and Public Health Implications

Arizona is home to many mosquito species, some of which are known vectors of infectious diseases that harm both humans and animals. Here, we provide an overview of the 56 mosquito species that have been identified in the State to date, but also discuss their known feeding preference and the diseases they can (potentially) transmit to humans and animals. This list is unlikely to be complete for several reasons: (i) Arizona's mosquitoes are not systematically surveyed in many areas, (ii) surveillance efforts often target specific species of interest, and (iii) doubts have been raised by one or more scientists about the accuracy of some collection records, which has been noted in this article. There needs to be an integrated and multifaceted surveillance approach that involves entomologists and epidemiologists, but also social scientists, wildlife ecologists, ornithologists, representatives from the agricultural department, and irrigation and drainage districts. This will allow public health officials to (i) monitor changes in current mosquito species diversity and abundance, (ii) monitor the introduction of new or invasive species, (iii) identify locations or specific populations that are more at risk for mosquito-borne diseases, and (iv) effectively guide vector control.

A Spatially Resolved and Environmentally Informed Forecast Model of West Nile Virus in Coachella Valley, California

West Nile virus (WNV) is the most significant arbovirus in the United States in terms of both morbidity and mortality. West Nile exists in a complex transmission cycle between avian hosts and the arthropod vector, Culex spp. mosquitoes. Human spillover events occur when humans are bitten by an infected mosquito and predicting these rates of infection and therefore the risk to humans may be associated with fluctuations in environmental conditions. In this study, we evaluate the hydrological and meteorological drivers associated with mosquito biology and viral development to determine if these associations can be used to forecast seasonal mosquito infection rates with WNV in the Coachella Valley of California. We developed and tested a spatially resolved ensemble forecast model of the WNV mosquito infection rate in the Coachella Valley using 17 years of mosquito surveillance data and North American Land Data Assimilation System-2 environmental data. Our multi-model inference system indicated that the combination of a cooler and dryer winter, followed by a wetter and warmer spring, and a cooler than normal summer was most predictive of the prevalence of West Nile positive mosquitoes in the Coachella Valley. The ability to make accurate early season predictions of West Nile risk has the potential to allow local abatement districts and public health entities to implement early season interventions such as targeted adulticiding and public health messaging before human transmission occurs. Such early and targeted interventions could better mitigate the risk of WNV to humans.

Development of multiplex PCR based NGS protocol for whole genome sequencing of West Nile virus lineage 2 directly from biological samples using Oxford Nanopore platform

West Nile virus (WNV) can affect humans, birds, horses and another mammals, causing asymptomatic infection, mild febrile disease, neurological and systematic disease and death. In order to gain insight into the prevalence of WNV, a monitoring program has been established in the Republic of Serbia. Whole genome sequencing is essential for the molecular epizootiological analysis of virus entry and transmission routes, especially in high-risk regions. This paper describes the development of a multiplex PCR based NGS protocol for whole genome sequencing of WNV lineage 2 directly from biological samples using Oxford Nanopore (ONT) platform. The results obtained using this platform, confirmed by Sanger sequencing, indicate that this protocol can be applied to obtain whole sequences of the WNV genome, even when the virus concentration in the sample is medium, Ct value is approximately 30. The use of this protocol does not require prior virus isolation on cell culture nor the depletion of host nucleic acids.

St. Louis Encephalitis Virus in the Southwestern United States: A Phylogeographic Case for a Multi-Variant Introduction Event

Since the reemergence of St. Louis Encephalitis (SLE) Virus (SLEV) in the Southwest United States, identified during the 2015 outbreak in Arizona, SLEV has been seasonally detected within Culex spp. populations throughout the Southwest United States. Previous work revealed the 2015 outbreak was caused by an importation of SLEV genotype III, which had only been detected previously in Argentina. However, little is known about when the importation occurred or the transmission and genetic dynamics since its arrival into the Southwest. In this study, we sought to determine whether the annual detection of SLEV in the Southwest is due to enzootic cycling or new importations. To address this question, we analyzed 174 SLEV genomes (142 sequenced as part of this study) using Bayesian phylogenetic analyses to estimate the date of arrival into the American Southwest and characterize the underlying population structure of SLEV. Phylogenetic clustering showed that SLEV variants circulating in Maricopa and Riverside counties form two distinct populations with little evidence of inter-county transmission since the onset of the outbreak. Alternatively, it appears that in 2019, Yuma and Clark counties experienced annual importations of SLEV that originated in Riverside and Maricopa counties. Finally, the earliest representatives of SLEV genotype III in the Southwest form a polytomy that includes both California and Arizona samples. We propose that the initial outbreak most likely resulted from the importation of a population of SLEV genotype III variants, perhaps in multiple birds, possibly multiple species, migrating north in 2013, rather than a single variant introduced by one bird.

Phylogenetic Analysis of the 2020 West Nile Virus (WNV) Outbreak in Andalusia (Spain)

During recent decades West Nile Virus (WNV) outbreaks have continuously occurred in the Mediterranean area. In August 2020 a new WNV outbreak affected 71 people with meningoencephalitis in Andalusia and six more cases were detected in Extremadura (south-west of Spain), causing a total of eight deaths. The whole genomes of four viruses were obtained and phylogenetically analyzed in the context of recent outbreaks. The Andalusian viral samples belonged to lineage 1 and were relatively similar to those of previous outbreaks which occurred in the Mediterranean region. Here we present a detailed analysis of the outbreak, including an extensive phylogenetic study. As part on this effort, we implemented a local Nextstrain server, which has become a constituent piece of regional epidemiological surveillance, wherein forthcoming genomes of environmental samples or, eventually, future outbreaks, will be included.

Sample pooling methods for efficient pathogen screening: Practical implications

Due to the large number of negative tests, individually screening large populations for rare pathogens can be wasteful and expensive. Sample pooling methods improve the efficiency of large-scale pathogen screening campaigns by reducing the number of tests and reagents required to accurately categorize positive and negative individuals. Such methods rely on group testing theory which mainly focuses on minimizing the total number of tests; however, many other practical concerns and tradeoffs must be considered when choosing an appropriate method for a given set of circumstances. Here we use computational simulations to determine how several theoretical approaches compare in terms of (a) the number of tests, to minimize costs and save reagents, (b) the number of sequential steps, to reduce the time it takes to complete the assay, (c) the number of samples per pool, to avoid the limits of detection, (d) simplicity, to reduce the risk of human error, and (e) robustness, to poor estimates of the number of positive samples. We found that established methods often perform very well in one area but very poorly in others. Therefore, we introduce and validate a new method which performs fairly well across each of the above criteria making it a good general use approach.

Epidemiological hypothesis testing using a phylogeographic and phylodynamic framework

Computational analyses of pathogen genomes are increasingly used to unravel the dispersal history and transmission dynamics of epidemics. Here, we show how to go beyond historical reconstructions and use spatially-explicit phylogeographic and phylodynamic approaches to formally test epidemiological hypotheses. We illustrate our approach by focusing on the West Nile virus (WNV) spread in North America that has substantially impacted public, veterinary, and wildlife health. We apply an analytical workflow to a comprehensive WNV genome collection to test the impact of environmental factors on the dispersal of viral lineages and on viral population genetic diversity through time. We find that WNV lineages tend to disperse faster in areas with higher temperatures and we identify temporal variation in temperature as a main predictor of viral genetic diversity through time. By contrasting inference with simulation, we find no evidence for viral lineages to preferentially circulate within the same migratory bird flyway, suggesting a substantial role for non-migratory birds or mosquito dispersal along the longitudinal gradient.

Twenty years of West Nile virus spread and evolution in the Americas visualized by Nextstrain

It has been 20 years since West Nile virus first emerged in the Americas, and since then, little progress has been made to control outbreaks caused by this virus. After its first detection in New York in 1999, West Nile virus quickly spread across the continent, causing an epidemic of human disease and massive bird die-offs. Now the virus has become endemic to the United States, where an estimated 7 million human infections have occurred, making it the leading mosquito-borne virus infection and the most common cause of viral encephalitis in the country. To bring new attention to one of the most important mosquito-borne viruses in the Americas, we provide an interactive review using Nextstrain: a visualization tool for real-time tracking of pathogen evolution (nextstrain.org/WNV/NA). Nextstrain utilizes a growing database of more than 2,000 West Nile virus genomes and harnesses the power of phylogenetics for students, educators, public health workers, and researchers to visualize key aspects of virus spread and evolution. Using Nextstrain, we use virus genomics to investigate the emergence of West Nile virus in the U S, followed by its rapid spread, evolution in a new environment, establishment of endemic transmission, and subsequent international spread. For each figure, we include a link to Nextstrain to allow the readers to directly interact with and explore the underlying data in new ways. We also provide a brief online narrative that parallels this review to further explain the data and highlight key epidemiological and evolutionary features (nextstrain.org/narratives/twenty-years-of-WNV). Mirroring the dynamic nature of outbreaks, the Nextstrain links provided within this paper are constantly updated as new West Nile virus genomes are shared publicly, helping to stay current with the research. Overall, our review showcases how genomics can track West Nile virus spread and evolution, as well as potentially uncover novel targeted control measures to help alleviate its public health burden.

Towards Translational Epidemiology: Next-Generation Sequencing and Phylogeography as Epidemiological Mainstays

Next-generation sequencing, coupled with the development of user-friendly software, has achieved a level of accessibility that is revolutionizing the way we approach epidemiological investigations. We can sequence pathogen genomes and conduct phylogenetic analyses to assess transmission, identify from which country or city a pathogen originated, or which contaminated potluck item resulted in widespread foodborne illness. However, until recently, these types of studies have been rarities, limited to specific investigations usually conducted over the short term. Given the feasibility and realized public health benefits of ascertaining pathogen relationships, federal, state, and county agencies are building their sequencing capacities, either through acquisition of equipment or collaborative activities. In this perspective, I detail research projects that our group collaborates on with county and state public health agencies, where the objective is to identify pathogen source locations with the longer-term goal of implementing proactive interventions.

An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar

How viruses evolve within hosts can dictate infection outcomes; however, reconstructing this process is challenging. We evaluate our multiplexed amplicon approach, PrimalSeq, to demonstrate how virus concentration, sequencing coverage, primer mismatches, and replicates influence the accuracy of measuring intrahost virus diversity. We develop an experimental protocol and computational tool, iVar, for using PrimalSeq to measure virus diversity using Illumina and compare the results to Oxford Nanopore sequencing. We demonstrate the utility of PrimalSeq by measuring Zika and West Nile virus diversity from varied sample types and show that the accumulation of genetic diversity is influenced by experimental and biological systems.