West Nile virus (WNV) —a mosquito-borne arbovirus—entered the USA through New York City in 1999 and spread to the contiguous USA within three years while transitioning from epidemic outbreaks to endemic transmission. The virus is transmitted by vector competent mosquitoes and maintained in the avian populations. WNV spatial distribution is mainly determined by the movement of residential and migratory avian populations. We developed an individual-level heterogeneous network framework across the USA with the goal of understanding the long-range spatial distribution of WNV. To this end, we proposed three distance dispersal kernels model: 1) exponential—short-range dispersal, 2) power-law—long-range dispersal in all directions, and 3) power-law biased by flyway direction —long-range dispersal only along established migratory routes. To select the appropriate dispersal kernel we used the human case data and adopted a model selection framework based on approximate Bayesian computation with sequential Monte Carlo sampling (ABC-SMC). From estimated parameters, we find that the power-law biased by flyway direction kernel is the best kernel to fit WNV human case data, supporting the hypothesis of long-range WNV transmission is mainly along the migratory bird flyways. Through extensive simulation from 2014 to 2016, we proposed and tested hypothetical mitigation strategies and found that mosquito population reduction in the infected states and neighboring states is potentially cost-effective.
West Nile disease (WND) is a vector-borne zoonosis which may result from infection by West Nile virus (WNV), a member of the family Flaviviridae, genus Flavivirus. This virus is the most common cause of arboviral disease in the United States [1]. From 1999 to 2017, more than 48 thousands WNV disease cases were reported to the Centers for Disease Control and Prevention (CDC) and more than two thousands of these reported cases resulted in death [2]. WNV is maintained in an enzootic transmission cycle between competent mosquitoes and birds. Birds are the reservoir and amplifying host for this virus. The US Centers for Diseases Control and Prevention (CDC) has identified WNV infection in more than three hundred species of birds. Infected bird movement is likely a key factor that affects the geographic spread of WNV, especially given the different habitats and routes of various species. Although many bird species may be infected with WNV, the American robin is considered an important amplifier of WNV and maybe a driver geographic spread because WNV-infected American robins have low mortality and high viremia [3,4]. Members of the Culex genus of mosquito are the principal vectors of this virus in the United States [5]. Humans, horses, and other mammals can be infected with WNV. However, these infections result in relatively low virus titers (viremia) therefore the infected animals and people are considered dead-end hosts (not capable of infecting feeding mosquitoes). Therefore, they do not have any epidemiological impact on WNV transmission or geographic spread [6]. To understand the transmission dynamics of WNV, several mathematical models have been developed [3,7–10]. These models predict the threshold conditions for WNV spreading in different scenarios. However, most of these models do not consider the spatial dynamics of WNV. Space or geographic spread has a significant role in WNV disease dynamics and modeling of WNV spatial spreading is complex because of the interactions of multiple potential mosquito vectors, avian amplifiers, and mammalian hosts. Liu et al. [9] developed a patchy model to analyze the spatial spreading of WNV, where patches are geographical space. They assumed patches are identical, spatial dispersal of birds and mosquitoes are symmetric within patches, and movement of birds and mosquitoes are only one-dimensional. According to this investigation, long-range dispersal of infected bird populations determines the spatial spread of WNV, not the dispersal of infected mosquito populations. Other investigators proposed a reaction-diffusion model [10], where they have spatially extended the non-spatial model of Wonham et al. [8] to mathematically estimate the spread of WNV. Here, diffusion terms in the reaction-diffusion partial differential equations represent vector mosquito and host bird population movements. They identified traveling wave solutions in their model and calculated the rate of spatial spread of infection. Durand et al. [11] developed a discrete time deterministic meta-population model in order to analyze the circulation of WNV between Southern Europe and West Africa. Another spatial model proposed by Maidana and Yang [12] used a system of partial differential reaction-diffusion equations. They also calculated the speed of disease dissemination by investigating the traveling wave solution of their model. They concluded, mosquito movements do not play an important role in disease dissemination. In addition, they included vertical transmission in their model and determined that vertical transmission is not an important factor for the spatial spread of WNV. Most WNV spread models are mathematical deterministic compartmental models. However WNV spread is highly stochastic because of the demography and movement of hosts and vectors varies between different locations. The major weaknesses of these models are the number and complexity of the compartments required to account for the many host and vector populations. In turn, the number of compartments increases the number of unknown parameters. Approximation of these parameters in any biological system is very challenging and prone to estimation errors which can create inaccuracies in the model outputs. We developed an individual-based heterogeneous network framework to understand WNV geographic spread. To build the network framework, we used the American Robin population density across the contiguous United States. The demographic characteristics of avian host populations and vector populations are not homogenous geographically, so we used a heterogeneous network framework. The transmission intensity of WNV depends on the abundance of WNV-infected vector mosquitoes in a given location. Mosquito population numbers fluctuate with local weather and season throughout the year, therefore we used a temperature dependent transmission rate. Although dead-end hosts cannot spread WNV to mosquitoes, we have quantified WNV case data only for humans, which we used to estimate unknown parameters. To understand the WNV spatial distribution, we proposed distance dispersal kernels, which describes the probability of dispersal with respect to distances. In this framework, we proposed three types of distance dispersal kernels: 1) exponential, 2) power-law, and 3) power-law biased by flyway. Then we compared the three distance kernels using approximate Bayesian computation based on sequential Monte Carlo sampling (ABC-SMC) method [13–18]. After conducting an extensive simulation for 2014-2016, we observed that an adapted fat-tailed or power-law kernel, which has long-distance links in specified directions can best describe the WNV human case data. We tested this network framework for the best kernel with the human case data and found that simulated results for more than 41 states of 49 states are consistent with the reported WNV cases. Our results support previous work on WNV spreading [3], which also modeled WNV spreading with migratory birds. We validate our work computationally from human incidence data. We proposed several theoretical mitigation strategies to control WNV and calculated their estimated costs. From the analysis of mitigation strategies, we suggest that potentially effective mitigation policies would include the application of mitigation control in areas with active transmission and in immediate neighboring states.
The study area of this research was the contiguous United States where WNV is considered endemic. We modeled WNV case distributions for 2014-2016. We used three data sets each year to develop our model. The first dataset contained the average monthly temperatures. Mosquito vector abundance correlated with temperature. Temperature data was from the National Centers for Environmental Information [19]. The second dataset contains American Robin population data from eBird [20]. This is a database for bird abundance and distribution, which is formed by the Cornell Lab of Ornithology and National Audubon Society. We used total observation of American Robin in each state of the USA for each month. The robin data set was used to train the network model. The American Robin is abundant throughout the United States and is a preferred food source for many WNV-competent mosquito species [21]. Based on host feeding patterns of the Culex genus of mosquitoes, robins are the most common WNV amplifying host [22–24]. Other important susceptible birds, such as American crow were not used because although they are an indicator species (high crow mortality), they are unlikely to spread virus geographically as they are mostly a residential species. In addition, as an indication of epidemic start point, we used WNV human incidence data. Many species of birds have long-distance migration during the spring and fall. Therefore the network does not focus on one long-distance migrating bird species but aggregates all species along the known flyways. To estimate model parameters we used human case data for WNV from CDC [2], which is the third dataset. To explore WNV long-distance spatial distribution in the USA, we used an individual-based heterogeneous network framework. In this framework, birds are on the individual level, a node represents an individual bird and connection between nodes is the possibility of virus dispersal from one infected bird to another susceptible bird by mosquito vectors. Links or connections are formed by movement of birds or movement of vectors. If there is no link between nodes then infected birds and insects are not moving virus between nodes. All virus transmission occurs by local competent vector mosquitoes. There is some evidence of bird-to-bird transmission, but it likely does not contribute to or maintain outbreaks. We split the bird population into four compartments; susceptible, exposed, infected, and recovered. Although, in the literature most mathematical models do not consider the exposed avian class when modeling WNV [8,12,25,26]. Birds transmit virus to mosquitoes when a susceptible mosquito vector takes an infected blood meal, then the mosquito becomes infectious after the extrinsic incubation period (EIP), or the time needed for the virus to spreads from the mosquito mid gut to the salivary glands; usually this process takes 7 to 14 days [3,27]. In addition, an infected bird can infect many mosquitoes simultaneously and also an infected mosquito can bite many susceptible or infected birds. Therefore, there is some delay in the system, to represent this delay we added the exposed class. We estimated exposed period from data by using the approximate Bayesian computation with sequential Monte Carlo sampling (ABC-SMC) method. After the exposed period, birds entered the infected compartment and an infected bird transitions to recovered after 4-5 days. To simulate this model, we used generalized epidemic mean-field (GEMF) framework developed by the Network Science and Engineering (NetSE) group at Kansas State University [28]. In GEMF, each node stays in a different state and the joint state of all nodes follows a Markov process [28–30]. The node level description of this Markov process is: P r [ x i (t + Δ t) = 1 | x i (t) = 0, X (t) ] = β (T) Y i Δ t (1) P r [ x i (t + Δ t) = 2 | x i (t) = 1, X (t) ] = λ Δ t + o (Δ t) (2) P r [ x i (t + Δ t) = 3 | x i (t) = 2, X (t) ] = δ Δ t + o (Δ t) (3) Here, X (t) is the joint state of all individual nodes at time t. xi (t) is a node state, xi (t) = C means node i is in C compartment at time t, C = 0,1, 2,3 corresponds to susceptible, exposed, infected, and recovered compartment. Yi is the number of infected neighbors of node i, β (T) is the transmission rate from one infected bird to one susceptible bird, which is a function of temperature, λ is the rate for exposed to infectious state, and finally, a node recovers from infectious state at a rate δ. For the spatial dynamic characteristics of WNV transmission, we built a network framework, which has 49 sub-networks one for each adjoining states of the contiguous United States plus the District of Columbia. The number of nodes in each sub-network is proportional to the size of the avian population in that state [20]. We considered the mosquito season June-October for the simulation period. Although the mosquito season is not the same for all states, mosquitoes are active from June to September in all of the states at these times [33]. The network for the avian population is (V, E). Here, V is the set of nodes, which is the union of nodes of all sub-network, V = SN1 ∪ SN2 ∪ SN3 ∪ ………‥ ∪ SN49, here SNi is a set of nodes in the sub-network i and E is the set of links among individual nodes. To build sub-networks, we used the total number of observations of American Robin for states per month in the simulation time period. | S N i | = max m j = m 1: m 2 (O B S m j i) * S c + N 0, here, O B S m j i is the total number of observations of American Robins in state i in month mj, N0 is the error term and N0 ∼ N (5,2) for this model. m1 is the first month after May and m2 is the last month before October when the average monthly temperature is greater than T0. Sc is the scaling constant. In each sub-network, we assumed that nodes are connected through Erdos-Renyi (n, p) random network topology [37]. In this network topology, we created links randomly among nodes with a probability p. Here, n is the number of nodes in a sub-network and p is the probability to form an edge. We set the probability p = R * log (n) /n, here R is a constant (R ≥ 2), as this value is more than the threshold value for the connectedness of an Erdos-Renyi graph [38], so nodes of a sub-network are locally connected. We will refer these networks as a local network in the subsequent sections of this paper. To build connections among sub-networks, we considered long-distance dispersal kernels [7,39], which describe the probability of dispersal with respect to distances. Dispersal kernels provide a simple model of dispersal to model dispersal events. For long-distance events, we used three types of kernel models; 1) Exponential, 2) power-law, and 3) power-law-flyway, which is a power-law kernel biased by flyway. The dispersal phenomenon in this work is not conserved because of long-distance movement of migratory birds and seasonality within bird populations. Some long-distance migratory birds can disperse outside the contiguous United States or outside the network nodes, which are discrete points. The connection probability between two nodes does not represent the probability that a single dispersal event happens rather it represents the probability of contact and subsequent pathogen transmission between them. A simple caricature of the network is shown in Fig 1. There are three sub-networks, A, B, and C. The links, which formed local networks are shown by solid lines. These links are introduced by Erdos-Renyi (n, p) network topology. Dashed lines are inter-links among sub-networks. These links established by using long-distance dispersal kernels. In this framework, we adopted approximate Bayesian computation based on a sequential Monte Carlo sampling (ABC-SMC) method for parameter estimation and model selection [13–18]. The role of mosquito populations in WNV transmission is expressed by disease transmission rate β. This framework used different transmission rates in different parts of the network corresponding to the local mosquito abundance. Using this heterogeneous feature in the framework, we evaluated theoretical mosquito population management measures to reduce the outbreak size or transmission rates in the state level. Some states such as Kansas, do not have statewide mosquito surveillance or management, but in these theoretical scenarios, it is assumed they can develop or benefit from effective statewide mosquito management programs. The framework will simply estimate how much the mosquito abundance is reduced or maintained based on the theoretical outcomes of coordinated control. Furthermore, we realize mosquito control is generally conducted on a county or municipal level, but the human case data is only available on a state level. Therefore the recommendations are for the lowest resolution of the data, which is state level but applies to counties and municipalities as well. If vector management is increased in a sub-network, then transmission rates will be changed by, β r = β R F, here βr is the reduced transmission rate and RF is the reduction factor. Then management costs will be Cost = RF * NSc, here NSc is the number of states where control measures were applied. We considered supplemental management measures with the existing management measures. We used two types of mitigation strategies across the United States, 1) dynamic infected place tracing strategy and 2) static ranked based strategy. In the infected place tracing, we traced the infected states, then plan the mitigation strategies according to them. For this type of mitigation strategies, we considered three cases; 1) case-1: only infected: applied control only in the infected states; 2) case-2: infected & first neighbors: applied control in the infected states with its first neighboring states (whose distance is less than 500km), and 3) case-3: infected & first neighbors & second neighbors: applied control in the infected states with its first neighboring states, and also with its second neighboring states (whose distance is in 500 − 1000km). For infected tracing control measure, we kept track of infected places monthly. If SNi sub-network is infected for month t, then control measures were applied for the month t + 1 based on these three cases. In the static ranked based mitigation strategy, we ranked the states by different variables (for example, temperature, size of the avian population etc.). For this strategy, we considered three cases; 1) temp. : states ranked by temperature, 2) pop. : states ranked by avian population size, and 3) temp. & pop. : states ranked by temperature and avian population size both, then we applied management measures in the top 30% of the states.
In this spatial-temporal individual-based heterogeneous network framework, we used three distance kernel models. The fundamental basic WNV epidemic model is the same for all the three network kernels. In the entire network, there are 49 sub-networks representing the 48 adjoining contiguous states plus the District of Columbia. All sub-network nodes are locally connected. The topology of the local network is Erdos-Renyi. The total nodes for the year 2015 was |V| = 7657 and the scaling constant is Sc = 0. 02. Here, E = El ∪ Edd; |El| is the number of total intra-links for all local networks, which is around 167000-170000 and |Edd| is the number of total inter-links among sub-networks. The description of sub-networks is provided in Table B in the S3 Text. We started the epidemic from states with the highest human incidence prior to June. We started the epidemic for the year 2015 by adding two infected nodes, one in sub-network SN4 (California) and another in sub-network SN42 (Texas). Connections among sub-networks are developed by distance dispersal kernels. Parameters for these kernels are estimated from the ABC-SMC method. To test the performance of this framework, we used estimated parameters from Table 1 for power-law kernel influenced by flyway. We set the parameters value; Kpl = 2. 3147, β0 = 0. 0059day-1, λ = 0. 0721day-1, and δ = 0. 2031day-1. The simulation period for the avian population model is from week-23 to week-44. The output of avian population was used as the input of zoonotic spillover compartment. Then we compared the output of zoonotic spillover compartment with human case data for week 24 to week 45. We considered a one-week lag between WNV incidence in birds and WNV incidence in humans. In humans, WNV-infected individuals (approximately 20%) develop a mild febrile illness after 3–6 days [56]. Peak of reporting of dead birds is one week prior than the reporting peak of human incidence [57]. In Fig 5, the mean simulated human case from the 49 sub-networks is compared with the weekly human case data for 2015 for the contiguous USA. The absolute errors between them are shown here. From this whisker plot, we can see that the median of the absolute error for the states is close to zero. In Fig 5, the largest outlier is California (marked by black circles). These outliers result from a mismatch between the simulated peak human incidence time and the observed human incidence peak time possibly because the very long state (north to south) has weather which is very different in southern California (warmer and drier) than northern California (cooler and wetter) causing a difference between peak mosquito seasons in the southern and northern parts. We compared the total yearly incidence of human WNV from this model with the state level reported case data. The results are shown in Fig 6. For 2015, we found that the case data for 42 of 49 locations were within the simulation results. The states where human cases were different from the simulation results were over-reported states (Nevada) and under-reported states (Louisiana, Mississippi, Nebraska, North Dakota, South Dakota, and Washington). The possible reason for this mismatch are reporting error or overwintering of virus in birds or mosquitoes or another bird species (not robins) is the key reservoir species for that state To build a disease prevalence map, we grouped the states in four categories; 1) higher prevalence —incidence is more than 100,2) intermediate prevalence—incidence is in between 50-99,3) moderate prevalence—incidence is in between 25-49 and 4) low prevalence—incidence is less than 25. To group the states, we used the median of the simulation results. The disease prevalence map from the model are presented in Fig 7a and from observed data are presented in Fig 7b. Among 49 locations, 40 locations are in the same prevalence group in both maps. We applied mitigation strategies on the power-law-flyway kernel network model to find the optimal mitigation plan. Fig 8a shows the number of infected states or epidemic size for dynamic infected places tracing. Epidemic size decreased faster with increased reduction factor for case-2 (infected & first neighbors) and case-3 (infected & first neighbors & second neighbors) than case-1 (only infected). The number of states where control measures were applied is displayed in Fig 9, which is proportional to cost. Therefore, the cost was minimal for case-2 than other two cases for RF > 2. From the cost analysis, we concluded that, although the cost for case-1 is less at the beginning of the yearly outbreak, we need to apply management only in the infected places, however by the end of the year the total cost for case-2 will smaller because of the smaller epidemic size. The results of the static ranked based mitigation strategy measure are presented in Fig 8b. We observed that, before RF = 4. 5, number of infected states for temp. & pop. dropped earlier than others. Number of infected states or epidemic size was smaller for temp. than pop. after RF > 3, infected population of a sub-network are more positively correlated with temperature. The NSc is always the same for these three cases. For all mitigation strategies, minimum epidemic size could be 2, as we started the epidemic from two states.
We proposed an individual-based heterogeneous network framework and tested three dispersal kernels to understand the spatial spread patterns of WNV human case data across the contiguous United States. This framework requires fewer parameters and has more flexibility to represent the spatial-temporal dynamics of WNV. Adding parameters can make the framework more realistic, for example, more competent bird species, landscape features for habitat preferences of host and vector species, daylight conditions [32], pathogen invasion from outside of USA, variable susceptibility among different hosts and vectors, WNV strain variability, mosquito and virus overwintering, vertical transmission, human movement characteristics etc. However, inclusion of too many factors increases model complexity which makes model optimization difficult given the availability of limited observational data. On the other hand, a simple model may insufficient to represent WNV spatial dynamics. Computational models need to be developed and parameters calculated with sufficient detail to be biologically accurate if they are used to evaluate epidemic management measures. However, for most biological systems, reliable parameter information is unknown. Unknown parameters or inaccurate assumptions add uncertainty to the model. Our framework has only four parameters to estimate (network Parameter K, transmission rate β, transition rate from exposed to infectious state, λ, and human spillover, η). This framework has compartments only for the avian population (susceptible, exposed, infected, and recovered), and is not species specific. We reduced the compartments for vector population by implementing them implicitly through transmission rate between infected nodes and susceptible nodes. The presented framework and dispersal kernel network model has an intermediate complexity that approximate Bayesian computation based on sequential Monte Carlo sampling (ABC-SMC) method successfully calibrated and estimated the parameters with the available data. If more data becomes available, it is possible to add them in this model for improved performance of the model. Furthermore, this framework is flexible and therefore can represent various hosts and vectors including with population seasonality, which plays an important role in WNV dynamics. For host population seasonality, we added a node property Activity, this property allows us to control active host populations in the network in a specific time period. We added vector seasonality in this framework with a temperature dependent transmission rate. This framework proposed one exponential and two fat-tailed distance kernel models for long-distance transmission of WNV with each model having increasing complexity and similarities to natural avian movement. WNV spatial distribution is very complex because WNV can infect more than 300 bird species, some of which are residential birds and short-distance migrators which disperse less than 500 km distances (short connections) whereas some species are long-distance migratory birds creating long connections. The long-distance migratory birds are the long-distance dispersal (LDD) agents for WNV. Previous studies tried to analyze spreading of WNV using a traveling wave with constant velocity, however, WNV spread more rapidly across the North America than would be expected from the assumption of constant velocity traveling wave [58]. Likely this is because traveling wave models unlike distance dispersal kernel models for WNV spreading do not capture the long-distance migrating birds which can have various migratory ranges and distances. Distance dispersal kernels have more flexibility to represent the different bird migration distances and can account for accelerating invasions. However, exponential kernels produce short-connections and therefore like traveling waves are limited to constant expansion, unlike fat-tailed power-law kernels which can generate accelerating invasions by creating the long-distance connections from migratory birds [59]. However, a general fat-tailed power-law kernel makes long-distance links in every direction which does not follow the incidence of WNV. Instead, a power-law-flyway kernel can be used to produce the long connections in the direction of flyways and short links in other directions. Bayesian inference was used to test which of the three kernel models best described WNV distribution on the network for three most recent years (2014- 2016). The power-law-flyway kernel best described the distribution of WNV cases because the long-range WNV transmission was concentrated mainly along the migratory bird flyways. The general power-law kernel overestimated the incidence data in some states because it was creating long-distance links in all directions. The performance for the power-law-flyway dispersal kernel model was evaluated for the three most recent years (2014-2016) when WNV was endemic in the USA. The observed case data for the 49 locations were within the range of the simulated results for 41 states for 2014 (Fig B in S1 Text), 42 states for 2015 (Fig 6), and 45 states for 2016 (Fig D in S1 Text). For all three years, the simulated results were similar to the observed data, except in Colorado, Louisiana, Mississippi, Nevada, Nebraska, North Dakota, and Washington. Nevada was over-reported for 2015 and all others were under-reported. The power law flyway dispersal kernel network model reported more WNV human incidence in Nevada than reported cases, one possible reason for over-reporting cases in Nevada has rural areas, which tend to under report human cases, whereas mosquito control districts and health departments, focused in urban areas, must test birds and mosquitoes, which explains why CDC reported WNV infected mosquitoes in 25% of counties in Nevada. The under-reported states had more human cases than predicted by the model. Under-reporting by the power-law-flyway kernel network model is likely because overwintering of the virus in some states (for example, Louisiana, Mississippi etc.), which was not considered. The overwintering infected Culex mosquitoes can stay in hibernacula such as sewers, houses, caves, and other warm areas in urban, suburban, and rural areas and initiate the outbreak in the spring. Furthermore, there may be under-reporting of cases by the model if robins are not the main reservoir species in a state, which would be predicted between gulf coast states (Louisiana and Mississippi) and northern states such as North and South Dakota and Washington. Stochastic simulations are useful tools to select the optimal future mitigation strategy after outbreaks of invasive species and pathogens. The foot-and-mouth disease (FMD) epidemics in 2001 in the United Kingdom developed by Keeling et al. [60], and mitigation strategies for pandemic influenza in the United States [61] are two well developed models with similarities to the current model. These models explore possible control measures such as culling, vaccination etc. for FMD [60], and vaccination, quarantine etc for influenza [61]. Most of these strategies can be examined with the network framework however, avian culling or vaccination for WNV control is not feasible. Vector control (or mosquito control) is a viable mitigation strategy for WNV, which is not considered by the other two models (FMD and influenza). To be applicable to any pathogens and inclusive of new mitigation methods, the mitigation strategies are non-specific and the predicted effectiveness of the mitigation methods can be adjusted to other methods. In the planning of the mitigation strategies, there is a trade-off between control measures effectiveness and their cost both monetary and loss of life. A stochastic simulation tool can decide the optimal mitigation strategy by dealing with this trade-off. Mitigation strategies for WNV were tested using the power-law-flyway dispersal kernel network model. The mosquito management measures are not specific to larvae or adults, rather simply generally accepted best practices to reduce mosquito abundance for the purpose of reducing pathogen transmission. The mitigation strategy analysis proposes supplemental measures in addition to the existing mosquito management in each state because the states had yearly reported WNV cases despite the existing management methods. To reduce WNV spread, a theoretical policy would be management in neighboring regions and not exclusively in the infected places. Although this approach can cost more at the beginning of the epidemic season however at the end, it can reduce total cost by decreasing the size of the epidemic. If management measures are applied only in the infected states, it is not possible to control the epidemic because of long-distance migratory birds. This is statewide management in a unified effort. We acknowledge that states do not conduct mosquito management in this way, but to test the spillover it was necessary to do the simulation in this way because only state-level data was available. The infected place tracing mitigation technique has been used to control other diseases (for example, FMD, influenza etc.), although their host population and control measure means are different, however, the main concept behind the mitigation techniques are similar. The findings from this research to control WNV epidemic can be useful to select optimal mitigation strategies for other pathogens. This research showed that the inclusion of directional long-distance dispersal of migratory birds improves model representations of the spatial patterns of WNV spread in the United States. The simulation of our framework in the context of long-distance directional dispersal suggested that cooperation and communication can facilitate early treatment and reduced outbreak sizes because of reduced WNV dispersal by American robins.
The underlying pattern of West Nile virus (WNV) geographic spread across the United States is not completely clear, which is a necessary step for continental or state level mitigation strategies to reduce WNV transmission. We report a network model that explains the geographic spread of WNV in the United States. West Nile virus is a mosquito-borne pathogen that infects many avian species with different movement ranges. From our research, we found that migration patterns and routes play an essential role in the WNV spatial distribution. The virus spreads in all directions at short distances because of local birds and short-distance migratory birds. However, the virus also disperses long distances along the avian migratory routes. Our model is designed to be flexible and therefore can be used to explore spreading patterns of other infectious diseases in other geographic locations.
invertebrates
medicine and health sciences
pathology and laboratory medicine
statistics
pathogens
microbiology
vertebrates
social sciences
animals
simulation and modeling
viruses
animal behavior
rna viruses
mathematics
network analysis
infectious disease control
insect vectors
zoology
research and analysis methods
infectious diseases
computer and information sciences
birds
medical microbiology
behavior
mathematical and statistical techniques
microbial pathogens
monte carlo method
disease vectors
insects
arthropoda
mosquitoes
psychology
animal migration
eukaryota
west nile virus
flaviviruses
viral pathogens
biology and life sciences
species interactions
physical sciences
amniotes
statistical methods
organisms
A spatio-temporal individual-based network framework for West Nile virus in the USA: Spreading pattern of West Nile virus