Mosquito-borne diseases are a global health priority disproportionately affecting low-income populations in tropical and sub-tropical countries. These pathogens live in mosquitoes and hosts that interact in spatially heterogeneous environments where hosts move between regions of varying transmission intensity. Although there is increasing interest in the implications of spatial processes for mosquito-borne disease dynamics, most of our understanding derives from models that assume spatially homogeneous transmission. Spatial variation in contact rates can influence transmission and the risk of epidemics, yet the interaction between spatial heterogeneity and movement of hosts remains relatively unexplored. Here we explore, analytically and through numerical simulations, how human mobility connects spatially heterogeneous mosquito populations, thereby influencing disease persistence (determined by the basic reproduction number R0), prevalence and their relationship. We show that, when local transmission rates are highly heterogeneous, R0 declines asymptotically as human mobility increases, but infection prevalence peaks at low to intermediate rates of movement and decreases asymptotically after this peak. Movement can reduce heterogeneity in exposure to mosquito biting. As a result, if biting intensity is high but uneven, infection prevalence increases with mobility despite reductions in R0. This increase in prevalence decreases with further increase in mobility because individuals do not spend enough time in high transmission patches, hence decreasing the number of new infections and overall prevalence. These results provide a better basis for understanding the interplay between spatial transmission heterogeneity and human mobility, and their combined influence on prevalence and R0.
Citation: Acevedo MA, Prosper O, Lopiano K, Ruktanonchai N, Caughlin TT, Martcheva M, et al. (2015) Spatial Heterogeneity, Host Movement and Mosquito-Borne Disease Transmission. PLoS ONE 10(6): e0127552. https://doi.org/10.1371/journal.pone.0127552
Academic Editor: James M. McCaw, The University of Melbourne, AUSTRALIA
Received: May 8, 2014; Accepted: April 16, 2015; Published: June 1, 2015
Copyright: © 2015 Acevedo et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: All relevant data are available via Github (https://github.com/maacevedo/Spatial-heterogeneity-host-movement-and-vector-borne-disease-transmission) and within the supporting information files.
Funding: This work was supported by National Science Foundation (NSF) Quantitative Spatial Ecology, Evolution, and Environment (QSE3) Integrative Graduate Education and Research Traineeship Program Grant 0801544, and NSF Doctoral Dissertation Improvement Grant (DEB-1110441). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
More than half of the world’s population is infected with some kind of vector-borne pathogen [1–3], resulting in an enormous burden on human health, life, and economies . Vector-borne diseases are most common in tropical and sub-tropical regions; however, their geographic distributions are shifting because of vector control, economic development, urbanization, climate change, land-use change, human mobility, and vector range expansion [5–9].
Mathematical models continue to play an important role in the scientific understanding of vector-borne disease dynamics and informing decisions regarding control [10–14] and elimination [15–17], owing to their ability to summarize complex spatio-temporal dynamics. Although there is increasing interest in the implications of spatial processes for vector-borne disease dynamics [18–22], most models that describe these dynamics assume spatially homogeneous transmission, and do not incorporate host movement [23–25]. Yet, heterogeneous transmission may be the rule in nature [26–28], where spatially heterogeneous transmission may arise due to spatial variation in vector habitat, vector control, temperature, and rainfall, influencing vector reproduction, vector survival and encounters between vectors and hosts [29, 30].
Movement of hosts among patches with different transmission rates links the pathogen transmission dynamics of these regions . In the resulting disease transmission systems some patches may have environmental conditions that promote disease transmission and persistence (i.e., hotspots), while other patches may not be able to sustain the disease without immigration of infectious hosts from hotspots . Control strategies often focus on decreasing vectorial capacity in hotspots [33, 34] with some successes, such as malaria elimination from Puerto Rico , and some failures [36, 37], such as malaria control efforts in Burkina Faso . An often overlooked factor when defining sites for control efforts is a patch’s connectivity to places of high transmission. For example, malaria cases during the 1998 outbreak in the city Pochutla, Mexico were likely caused by human movement into the city from nearby high transmission rural areas, despite active vector control in Pochutla . Understanding the interaction between connectivity—defined by the rate of movement of hosts among patches—and spatial heterogeneity in transmission via mathematical models has the potential to better inform control and eradication strategies of mosquito-borne diseases in real-world settings [37, 40].
In this study, we ask, how host movement and spatial variation in transmission intensity influense malaria long-term persistence and prevalence. First, we show analytically that transmission intensity is an increasing function of spatial heterogeneity in a two-patch system, where the patches are connected by host movement. Second, we apply a multi-patch adaptation of the Ross-Macdonald modeling framework for malaria dynamics to explore the implications of spatial heterogeneity in transmission intensity and human movement for disease prevalence and persistence. The mosquitoes that transmit malaria typically move over much smaller spatial scales than their human hosts. Thus, we assume that mosquito populations are isolated in space. The varying size of mosquito populations across a landscape introduces spatial heterogeneity in transmission intensity. This heterogeneity, coupled with the fact that humans commonly move among areas with varying degrees of malaria transmission, makes malaria an ideal case study.
Materials and Methods
The Ross-Macdonald modeling approach describes a set of simplifying assumptions that describe mosquito-borne disease transmission in terms of epidemiological and entomological processes . Although it was originally developed to describe malaria dynamics, the modeling framework is simple enough to have broad applicability to other mosquito-borne infections. One of the most important contributions of the Ross-Macdonald model is the identification of the threshold parameter for invasion R0, or the basic reproductive number. Threshold quantities, such as R0, often form the basis of planning for malaria elimination. In some cases R0 also determines the long-term persistence of the infection. Here, we define persistence to mean uniform strong persistence of the disease; that is whether the disease will remain endemic in the population, and bounded below by some positive value, over the long term. Mathematically, a disease is uniformly strongly persistent if there exists some ϵ > 0 such that limsupt → ∞ I(t) ≥ ϵ for any I(0) > 0, where I(t) is the number of infected individuals at time t [42, 43].
To extend the Ross-Macdonald model to a landscape composed of i = 1, …, Q patches we need to account for the rate of immigration and emigration of humans among the Q patches. The full mathematical derivation of the multi-patch extension (Eq 1) from the original Ross-Macdonald model can be found in S1 Text.
For each patch i, the rates of change in the proportion of infected mosquitoes, the number of infected hosts, and the total number of humans are calculated as where Ni describes the total size of the human population in patch i, Ii represents the number of infected hosts in patch i, zi represents the proportion of infected mosquitoes in patch i, and kji represents the rate of movement of human hosts from patch i to patch j. Note that 1/kji describes the amount of time (days in this particular parameterization) an individual spends in patch i before moving to patch j. For simplicity, we assumed that the rate of host movement was symmetric between any two patches, and equal amongst all patches, such that k = kij = kji. We further assumed that the initial human population densities for each patch were equal. This constraint on the initial condition, along with the assumption of symmetric movement, causes the population size of each patch to remain constant, that is, dNi/dt = 0 for all i. We also assumed that the only parameter that varies among patches is the ratio of mosquitoes to humans, mi. The rate ai at which mosquitoes bite humans, the probability ci a mosquito becomes infected given it has bitten an infected human, the probability bi a susceptible human is infected given an infectious mosquito bite, the mosquito death rate gi, the human recovery rate ri, and the extrinsic incubation period (the incubation period for the parasite within the mosquito) ni, are all assumed constant across the landscape. Consequently, for all i = 1, …, Q, ai = a, bi = b, ci = c, gi = g, ri = r, and ni = n.
In this model there is no immunity conferred after infection. Furthermore, although host demography (births and deaths) can play an important role in transient disease dynamics, because our focus is the relationship between equilibrium prevalence and R0 under the assumption of constant patch population sizes, we omit host demography. Choosing constant birth rates Λ = μN and natural host mortality rates μ in each patch yields identical R0 and equilibria to our model, with the exception that r is replaced by r + μ. Thus, including host demography in this way would result in a slight decrease in R0 and prevalence by decreasing the infectious period. How host demography influences the relationship between R0 and prevalence when patch population sizes are not constant, and moreover, when host demography is heterogeneous, is an interesting question that remains to be explored. These simplifying assumptions yield the following system of 2Q equations, (1)
Differences in the ratio of mosquitoes to humans, mi results in a network of heterogeneous transmission, where each patch in the network is characterized by a different transmission intensity. The basic reproduction number for an isolated patch (i.e., one not connected to the network through human movement) is defined by , where αi: = mi abe−gn and β: = ac, and is a measure of local transmission intensity. Furthermore, R0,i is a threshold quantity determining whether disease will persist in patch i in the absence of connectivity. In particular, if R0,i > 1, malaria will persist in patch i, while if R0,i ≤ 1, it will go extinct in the absence of connectivity with other patches. R0,i (local transmission) increases with the ratio of mosquitoes to humans mi, and if more transmission occurs, more people are infected at equilibrium. These results, however, do not necessarily hold in a network where hosts move among patches . Indeed, movement can cause the disease to persist in a patch where it would otherwise die out [20, 44].
To address this limitation of the isolated patch reproduction number, we used the next generation approach [45, 46] to calculate R0 for the whole landscape. This approach requires the construction of a matrix K = FV−1, where J = F−V is the Jacobian of the 2Q-dimensional system evaluated at the disease-free equilibrium, F is nonnegative, and V is a nonsingular M-matrix. F contains terms related to new infection events, and V contains terms of the Jacobian related to either recovery or migration events. This choice satisfies the conditions for the theory to hold, and the important consequence of this approach is that the spectral radius of the next generation matrix ρ(K) is less than one if and only if the disease-free equilibrium is locally asymptotically stable. Defining R0 = (ρ(K))2, we have that the disease-free equilibrium is locally asymptotically stable when R0 < 1 and unstable when R0 > 1. We proved (see S2 Text) that System (1) exhibits uniform weak persistence of the disease when R0 > 1; that is, when R0 > 1, there exists an ϵ > 0 such that , for any initial condition for which . Furthermore, because our model is an autonomous ordinary differential equation, uniform weak persistence implies uniform strong persistence. Consequently, when R0 > 1, there exists an ϵ > 0 such that , for any initial condition for which [42, 43]. A generalization of our multi-patch system (see System (8) in ) exhibits a unique endemic equilibrium when R0 > 1 which is globally asymptotically stable. Likewise, the disease-free equilibrium for their model is globally asymptotically stable when R0 ≤ 1. In fact, Auger et al.  proved this result even when migration is neither constant across the landscape, nor symmetric.
Because R0,i defines a threshold for disease persistence in an isolated patch and R0 defines a threshold for disease persistence in the connected network, we use these two quantities as surrogates for local patch persistence when patches are isolated, and persistence in the connected network as a whole, respectively. Prevalence, on the other hand, was calculated as the total proportion of infected hosts in the landscape at equilibrium.
Heterogeneity in transmission intensity was quantified using the coefficient of variation (CV) of the ratio of mosquito to humans (m) such that (2) where describes the average ratio of mosquito to humans in the landscape and represents the standard deviation associated with this average. This coefficient of variation is a simple measure commonly used in landscape ecology to quantify landscape heterogeneity .
We analyze two cases: (1) a simple two-patch system (Q = 2) where we study analytically the relationship between spatial heterogeneity, R0 and prevalence. Then, (2) we address a similar question in a multi-patch system (Q = 10) where each patch is characterized by their unique transmission intensity (see below).
We use an analytical approach (see S3 Text) to study the relationship between R0, prevalence, and spatial heterogeneity in the special case where the network is composed of two connected patches (Q = 2). Transmission heterogeneity in the system is created by choosing different values for m1 and m2, the ratio of mosquitoes to humans in the two patches, and quantified by the coefficient of variation, CV. We define to be the average of m1 and m2, and study the behavior of R0 and prevalence as CV increases.
To study the implications of spatial heterogeneity in transmission intensity, in the presence of host movement, for disease prevalence and persistence, we generated a landscape composed of Q = 10 discrete patches connected by movement (Fig 1). We used this landscape to simulate a spatially homogeneous configuration in transmission intensity and four heterogeneous configurations (Fig 1). As with the two-patch analysis, the variation in transmission intensity was attained by varying the ratio of mosquitoes to humans mi, while keeping all other parameters constant (Table 1). The ratio of mosquitoes to humans in each patch was drawn from a normal distribution such that in the homogeneous configuration mi = 60, and in the four heterogeneous configurations , , , and . This resulted in the same mean transmission intensity in each of the landscape configurations (), although the range (min R0,i, max R0,i varied among the five configurations: [2.17, 2.17], [1.04, 3.33], [0.03, 4.66], [0.03, 5.96], and [0.03, 6.83] from the homogeneous landscape to the most heterogeneous configuration, respectively (Fig 1). This resembles, in part, variation in malaria transmissibility reported in South America and Africa . To determine how host movement affected persistence and prevalence, and how their relationship depended upon variation in patch transmissibility, we varied the rate of host movement between all patches (k) from 0 to 0.2 (days−1) in 1 × 10−2 increments. This rate was equal among all patches. Given that population size was also equal among patches we are evaluating the simple case where population size is constant and movement is symmetric among patches. We replicated this simulation 100 times for each configuration.
Nodes represent patches characterized by their randomly generated R0,i, and links represent host movement. Each configuration represents a particular scenario of spatial heterogeneity in transmission intensity, which increases with increasing coefficient of variation (CV).
To evaluate the effect of heterogeneity in transmission intensity on disease dynamics, we first proved analytically for the two-patch model that the network reproduction number R0, and the total disease prevalence limt → ∞(I1(t)/N+I2(t)/N) increase with variance , even if , and consequently the average transmission intensity (R01+R02)/2 between the two regions, remains constant (see Theorems 0.0.2 and 0.0.4 in S3 Text). Because CV is proportional to the square root of the variance V, this implies that disease persistence and prevalence increase with CV. However, the influence of heterogeneity on R0 becomes less profound as connectivity between the two patches increases (see Proposition 0.0.3 in S3 Text).
Spatial heterogeneity in transmission intensity increased long-term persistence of infection (R0) in the multi-patch system (Fig 2). Yet, increasing host movement-rate decreased R0 in the spatially heterogeneous scenarios. Spatial homogeneity resulted in the lowest R0 of all landscape configurations (Fig 2), which is consistent with our conclusions derived analytically from the two-patch system (see above). R0 in this homogeneous case was also independent of movement because the system was effectively a one patch system. In contrast, in all heterogeneous configurations, increasing host movement-rate resulted in a decrease in R0 that approached an asymptote. The value of this asymptote increased with increasing spatial heterogeneity (Fig 2), which is also consistent with our analytic results for the two-patch case.
Lines represent means and shaded areas 95% confidence intervals. Spatial heterogeneity in transmission intensity increases with the coefficient of variation (CV). (c) Box-plots shows the distribution of patch-specific transmission intensities R0,i in 100 simulations for each level of spatial heterogeneity. Note how variance increases with CV, while the average remains similar among configurations.
Similarly, spatial heterogeneity in transmission intensity increased disease prevalence in the multi-patch system. Spatial homogeneity in transmission intensity resulted in the lowest prevalence of all landscape configurations (Fig 2). Maximum prevalence and the asymptotic prevalence with increasing spatial heterogeneity in transmission intensity, which again, agrees with our conclusions derived for the two-patch case. Disease prevalence initially increased with increasing movement, was maximized at relatively low movement rates and later decreased. The movement rate, k, that maximized prevalence increased with increasing heterogeneity and occurred at movement rates corresponding to once every 0.5 to 1.5 years. This suggests that the rate of movement required to maximize disease prevalence increases with increasing spatial heterogeneity in transmission intensity. Note that, in the simulations, mean R0,i remained the same for all scenarios while variance increased with increasing coefficient of variation, as expected (Fig 2). In all heterogeneous configurations prevalence and R0 followed a non-monotonic relationship in the presence of host movement (Fig 3).
We have explored the way that disease prevalence and R0 — two important measures of mosquito-borne pathogen transmission — display a complex non-monotonic relationship as a result of spatial heterogeneity in mosquito density and human mobility. Heterogeneity in mosquito density and mosquito bionomic patterns affecting vectorial capacity drive spatially heterogeneous biting patterns, while human mobility connects isolated areas that can have very different mosquito populations. We illustrated these patterns analytically in a two-patch system, and numerically in a multi-patch extension of the Ross-Macdonald modeling framework. We showed that prevalence was maximized at low rates of movement, whereas R0 always decreased with increasing movement rates. These results suggest that the relationship between R0 and prevalence is intimately intertwined with the interaction between host movement and the degree of spatial heterogeneity in a region.
Transmission heterogeneity generally promotes persistence in host-parasite systems [18, 58–61]. This heterogeneity may have a spatial component arising from spatial variation in factors affecting mosquito ecology such as habitat distribution or host finding ability [25, 61]. Our results showed that disease persistence decreased with increasing rates of movement even in highly spatially heterogeneous landscapes with multiple transmission hotspots (Figs 1 and 2). At low rates of movement, transmission was highly heterogeneous, with high rates of transmission in some patches and low in others. R0 was higher in this scenario, because our calculation of R0 describes the average number of potential infections that arise from an average infected host in the system and thus its magnitude is being influenced by conditions in high transmission patches (Fig 4). Transmission becomes more homogeneous with increasing rate of movement resulting in individual patch transmissibility more similar to the overall average (Fig 4). A similar result was found in a study of the metapopulation dynamics of Schistosomiasis (bilharzia) , where increased social connectivity sometimes reduced large-scale disease persistence because as mobility increases infectious individuals spent less time in areas of high transmission distributing infection away from hotspots. Thus, acknowledging host movement patterns is required to better understand disease persistence in heterogeneous landscapes.
The proportion of infected hosts in the low transmission patch increase with increasing rate of movement because it is receiving infected immigrants from other patches with high transmission. The proportion of infected hosts in the high transmission patch decrease with increasing rate of movement because of increasing emigration of infected hosts to other patches.
Results from our numerical simulations support previous theoretical and empirical work showing that disease prevalence is generally maximized at low to intermediate levels of movement [31, 63, 64]. Our results add to this body of theory by showing that the amount of movement required to achieve peak prevalence increases with increasing spatial transmission heterogeneity. At very low rates of movement, individuals spend most of their time in a single patch. In transmission hotspots most hosts are already infected at equilibrium and most bites do not yield new infections. A relatively small increase in movement will significantly increase the number of hosts exposed to very intense transmission (Fig 4). Therefore, as connectivity increases, the number of infectious bites in high transmission patches decrease, yet, this decrease is offset by the increase in the number of susceptibles that visit these patches. As connectivity continues to increase, hosts spend less time in high transmission patches resulting in a decrease in the number of hosts that become infected in high transmission patches. This causes the number of infectious bites in high transmission patches to decline, ultimately causing fewer people to be infected, and prevalence decreases. The different behaviors of prevalence and R0 in the presence of spatial heterogeneity and mobility suggest a role for models including mobility and spatial scale in the estimation of prevalence based on R0 estimates, because the assumed positive relationship between the two is disrupted .
Reproduction numbers (R0) are useful to understand the intensity of transmission in a region and are often used to design and evaluate control measures of mosquito-borne diseases. The estimation of R0 can be done using several different methods, including estimating number of infectious bites on a person per year [1, 61, 65, 66]. Generally, depending on the assumptions about superinfections and density dependence among parasites, R0 is proportional to the inverse of the fraction of uninfected individuals at equilibrium (i.e. R0 and prevalence are positively correlated) [67, 68]. Yet, this relationship between prevalence and R0 has been shown to be disrupted by heterogeneous biting [18, 58, 61, 67–69]. Our analysis of the two-patch system illustrated that increasing heterogeneity increases both prevalence and R0, but the multi-patch numerical simulations show this effect is diminished as connectivity increases suggesting that the human“activity space” — or how humans spend time between areas of varying mosquito densities — is also an important determinant of the relationship between R0 and prevalence . For example, assuming that transmission intensity across two regions is the average of the transmission intensity in each region will underestimate the disease burden, particularly at low to intermediate levels of connectivity. Therefore our results emphasize the necessity for reasonable estimates of host movement rates, because individual patch transmission intensities do not uniquely determine overall transmission intensity and prevalence.
Our findings have important practical implications for mosquito-borne disease control in heterogeneous landscapes in the presence of symmetric host movement. Our results show that the dynamics of spatially heterogeneous system are driven primarily by the characteristics of areas with the highest potential for transmission by mosquitoes, which supports the idea that hotspots should be targeted for control efforts. If control strategies are untargeted these high transmission areas may represent residual areas where the disease persists with the potential to re-colonize others [32, 71, 72], or maintain transmission throughout the system. This is shown by the persistence of malaria in many landscape scenarios, despite R0,i < 1 in many patches (Fig 2a and 2c). Thus, controlling malaria transmission in areas with heterogeneous transmission requires a combination of interventions that include mosquito control, the reduction of human infectious reservoirs, and vaccination targeted towards high transmission areas .
Finally, human movement between areas often changes over time, and predicting how these changes will affect transmission and prevalence requires understanding the effect of connectivity on prevalence and the initial degree of movement. If human movement is very low initially, an increase in movement is likely to increase endemic prevalence, while an initially high human movement will likely result in a decrease in endemicity if movement increases further. Therefore, knowing the degree of connectivity between areas and how connectivity changes over time is also important to management and elimination planning . Recent studies are beginning to analyze human movement in relation to mosquito-borne pathogen transmission [70, 73–75], and these show great promise for improving models of mosquito-borne pathogen transmission across geographic scales.
S1 Text. Multi-patch model derivation.
Derivation of a multi-patch extension of the Ross-Macdonald model in Eq (1) from a single-patch model.
S2 Text. Theorem 0.0.1.
Mathematical proof showing that system of equations in (1) exhibit uniform weak persistence.
This study greatly benefited from insightful discussions with A. Tatem and C. Cosner. We also thank three anonymous reviewers for their insightful and helpful comments on pervious version of this manuscript. Funding was provided by the National Science Foundation (NSF) Quantitative Spatial Ecology, Evolution, and Environment (QSE3) Integrative Graduate Education and Research Traineeship Program Grant 0801544 at the University of Florida. MAA was also supported by an NSF Doctoral Dissertation Improvement Grant (DEB-1110441).
Conceived and designed the experiments: MAA OP KL NR TC MM CO DLS. Performed the experiments: MAA OP KL NR. Analyzed the data: MAA OP KL NR. Contributed reagents/materials/analysis tools: MAA OP KL NR TC MM CO DLS. Wrote the paper: MAA OP KL NR TC MM CO DLS.
- 1. Gething PW, Patil AP, Smith DL, Guerra CA, Elyazar IR, Johnston GL, et al. A new world malaria map: Plasmodium falciparum endemicity in 2010. Malaria Journal. 2011;10:378. pmid:22185615
- 2. Gething PW, Elyazar IRF, Moyes CL, Smith DL, Battle KE, Guerra CA, et al. A long neglected world malaria map: Plasmodium vivax endemicity in 2010. PLoS Neglected Tropical Diseases. 2012;6(9):e1814. pmid:22970336
- 3. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496:504–507. pmid:23563266
- 4. Varmus H, Klausner R, Zerhouni E, Acharva T, Daar AS, Singer PA. Public Health. Grand challenges in global health. Science. 2003;302:398–399. pmid:14563993
- 5. Hay SI, Guerra CA, Tatem AJ, Noor AM, Snow RW. The global distribution and population at risk of malaria: past, present, and future. The Lancet Infectious Diseases. 2004;4(6):327–336. pmid:15172341
- 6. Cohen JM, Smith DL, Cotter C, Ward A, Yamey G, Sabot OJ, et al. Malaria resurgence: a systematic review and assessment of its causes. Malaria Journal. 2012;11(1):122. pmid:22531245
- 7. Chiyaka C, Tatem A, Cohen J, Gething P, Johnston G, Gosling R, et al. The stability of malaria elimination. Science. 2013;339(6122):909–910. pmid:23430640
- 8. Smith DL, Cohen JM, Chiyaka C, Johnston G, Gething PW, Gosling R, et al. A sticky situation: the unexpected stability of malaria elimination. Philosophical Transactions of the Royal Society B: Biological Sciences. 2013;368(1623):20120145.
- 9. Tatem AJ, Gething PW, Smith DL, Hay SI. Urbanization and the global malaria recession. Malaria Journal. 2013;12(1):133. pmid:23594701
- 10. Maudea RJ, Lubella Y, Socheatd D, Yeunge S, Saralambaa S, Pongtavornpinyoa W, et al. The role of mathematical modelling in guiding the science and economics of malaria elimination. International Health. 2010;2(4):239–246.
- 11. McKenzie FE, Samba EM. The role of mathematical modeling in evidence-based malaria control. American Journal of Tropical Medicine and Hygiene. 2004;71(Suppl 2):94–96. pmid:15331824
- 12. McKenzie FE. Why model malaria? Parasitology Today. 2000;16(12):511–516.
- 13. Smith DL, Battle KE, Hay SI, Barker CM, Scott TW, McKenzie FE. Ross, Macdonald, and a theory for the dynamics and control of mosquito-transmitted pathogens. PLoS Pathogens. 2012;8(4):e1002588. pmid:22496640
- 14. Reiner RC, Perkins TA, Barker CM, Niu T, Chaves LF, Ellis AM, et al. A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970–2010. Journal of The Royal Society Interface. 2013;10(81):20120921.
- 15. Anderson RM, May RM. Infectious Diseases of Humans. Dynamics and Control. Oxford University Press; 1992.
- 16. Zanzibar Malaria Control Program. Malaria elimination in Zanzibar: a feasibility assessment; 2009. [Online; accessed October-2009]. http://www.malariaeliminationgroup.org/sites/default/files/MalariaEliminationZanzibar.pdf.
- 17. Alonso PL, Brown G, Arevalo-Herrera M, Binka F, Chitnis C, Collins F, et al. A research agenda to underpin malaria eradication. PLoS Medicine. 2011;8(1):e1000406. pmid:21311579
- 18. Hasibeder G, Dye C. Population dynamics of mosquito-borne disease: persistence in a completely heterogeneous environment. Theoretical Population Biology. 1988;33(1):31–53. pmid:2897726
- 19. Rodríguez DJ, Torres-Sorando L. Models of infectious diseases in spatially heterogeneous environments. Bulletin of Mathematical Biology. 2001;63:547–571. pmid:11374305
- 20. Cosner C, Beier JC, Cantrella RS, Impoinvilc D, Kapitanskia L, Pottsd MD, et al. The effects of human movement on the persistence of vector-borne diseases. Journal of Theoretical Biology. 2009;258(4):550–560. pmid:19265711
- 21. Perkins TA, Scott TW, Le Menach A, Smith DL. Heterogeneity, mixing, and the spatial scales of mosquito-borne pathogen transmission. PLoS Computational Biology. 2013;9(12):e1003327. pmid:24348223
- 22. Qiu Z, Kong Q, Li X, Martcheva M. The vector-host epidemic model with multiple strains in a patchy environment. Journal of Mathematical Analysis and Applications. 2013;405(1):12–36.
- 23. Grenfell BT, Dobson AP. Ecology of infectious disease in natural populations. Cambridge University Press; 1995.
- 24. Torres-Sorando L, Rodríguez DJ. Models of spatio-temporal dynamics in malaria. Ecological Modelling. 1997;104:231–240.
- 25. Smith DL, Dushoff J, McKenzie FE. The risk of mosquito-borne infection in a heterogeneous environment. PLoS Biology. 2004;2(11):e368. pmid:15510228
- 26. Woolhouse ME, Dye C, Etard JF, Smith T, Charlwood JD, Garnett GP, et al. Heterogeneities in the transmission of infectious agents: implications for the design of control programs. Proceedings of the National Academy of Sciences USA. 1997;94(1):338–342.
- 27. Shaw DJ, Grenfell BT, Dobson AP. Patterns of macroparasite aggregation in wildlife host populations. Parasitology. 1998;117:597–610. pmid:9881385
- 28. Smith DL, Drakeley CJ, Chiyaka C, Hay SI. A quantitative analysis of transmission efficiency versus intensity for malaria. Nature Communications. 2010;1:108. pmid:21045826
- 29. Mbogo CM, Mwangangi JM, Nzovu J, Gu W, Yan G, Gunter JT, et al. Spatial and temporal heterogeneity of Anopheles mosquitoes and Plasmodium falciparum transmission along the Kenyan coast. American Journal of Tropical Medicine and Hygien. 2003;68(6):734–742.
- 30. Kazembe LN, Kleinschmidt I, Holtz TH, Sharp BL. Spatial analysis and mapping of malaria risk in Malawi using point-referenced prevalence of infection data. International Journal of Health Geographics. 2006;5:41. pmid:16987415
- 31. Adams B, Kapan DD. Man bites mosquito: Understanding the contribution of human movement to vector-borne disease dynamics. PLoS ONE. 2009;4(8):e6763. pmid:19707544
- 32. Bousema T, Griffin JT, Sauerwein RW, Smith DL, Churcher TS, Takken W, et al. Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Medicine. 2012;9(1):e1001165. pmid:22303287
- 33. Hotez PJ, Molyneux DH, Fenwick A, Kumaresan J, Ehrlich S, Sachs JD, et al. Control of neglected tropical diseases. The New England Journal of Medicine. 2007;357(10):1018–1027. pmid:17804846
- 34. Lambrechts L, Knox TB, Wong J, Liebman KA, Albright RG, Stoddard ST. Shifting priorities in vector biology to improve control of vector-borne disease. Tropical Medicine and International Health. 2009;14(12):1505–1514. pmid:19807899
- 35. Miranda-Franco R, Casta-Vélez A. La erradicación de la malaria en Puerto Rico. Revista Panamericana de Salud Pública. 1997;2(2):146–150. pmid:9312420
- 36. Stratton L, O’Neill MS, Kruk ME, Bell ML. The persistent problem of malaria: Addressing the fundamental causes of a global killer. Social Science and Medicine. 2008;67(5):854–862. pmid:18583009
- 37. Barbu C, Dumonteil E, Gourbière S. Evaluation of spatially targeted strategies to control non-domiciliated Triatoma dimidiata vector of Chagas disease. PLoS Neglected Tropical Diseases. 2011;5(5):e1045. pmid:21610862
- 38. Baragatti M, Fournet F, Henry MC, Assi S, Ouedraogo H, Rogier C, et al. Social and environmental malaria risk factors in urban areas of Ouagadougou, Burkina Faso. Malaria Journal. 2009;8:30.
- 39. Hernández-Avila JE, Rodríguez MH, Betanzos-Reyes AF, Danis-Lozano R, Méndez-Galván J, Vel azquez-Monroy OJ, et al. Determinant factors for malaria transmission on the coast of Oaxaca State, the main residual transmission focus in México. Salud Pública de México. 2006;48(5):405–417. pmid:17063824
- 40. Grenfell B, Harwood J. (Meta)population dynamics of infectious diseases. Trends in Ecology and Evolution. 1997;12(10):395–399. pmid:21238122
- 41. Smith DL, McKenzie FE, Snow RW, Hay SI. Ross, Macdonald, and a theory for the dynamics and control of mosquito-transmitted pathogens. PLoS Pathogens. 2012;8:e1002588. pmid:22496640
- 42. Freedman H, Moson P. Persistence definitions and their connections. Proceedings of the american mathematical society. 1990;109(4):1025–1033.
- 43. Thieme HR. Persistence under relaxed point-dissipativity (with application to an endemic model). SIAM Journal on Mathematical Analysis. 1993;24(2):407–435.
- 44. Prosper O, Ruktanonchai N, Martcheva M. Assessing the role of spatial heterogeneity and human movement in malaria dynamics and control. Journal of Theoretical Biology. 2012;303:1–14. pmid:22525434
- 45. Diekmann O, Heesterbeek JAP, Metz JAJ. On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. Journal of Mathematical Biology. 1990;28:365–382. pmid:2117040
- 46. van den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Mathematical Biosciences. 2002;180(1–2):29–48. pmid:12387915
- 47. Auger P, Kouokam E, Sallet G, Tchuente M, Tsanou B. The Ross-Macdonald model in a patchy environment. Mathematical Biosciences. 2008;216(2):123–131. pmid:18805432
- 48. Li H, Reynolds J. On definition and quantification of heterogeneity. Oikos. 1995;73:280–284.
- 49. Loyola EG, González-Cerón L, Rodríguez MH, Arredondo-Jiménez JI, Bennett S, Bown DN. Anopheles albimanus (Diptera: Culicidae) Host Selection Patterns in Three Ecological Areas of the Coastal Plains of Chiapas, Southern Mexico. Journal of Medical Entomology. 1993;30:518–523. pmid:8510111
- 50. Beier JC, Davis JR, Vaughan JA, Noden BH, Beier MS. Quantitation of Plasmodium falciparum sporozoites transmitted in vitro by experimentally infected Anopheles gambiae and Anopheles stephensi. American Journal of Tropical Medicine and Hygiene. 1991;44:564–570. pmid:2063960
- 51. Bonnet S, Gouagna C, Safeukui I, Meunier JY, Boudin C. Comparison of artificial membrane feeding with direct skin feeding to estimate infectiousness of Plasmodium falciparum gametocyte carriers to mosquitoes. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2000;94:103–106. pmid:10748913
- 52. Collins WE, Warren M, Skinner JC, Richardson BB, Kearse TS. Infectivity of the Santa Lucia (El Salvador) strain of Plasmodium falciparum to different anophelines. Journal of Parasitology. 1977;63:57–61. pmid:403273
- 53. Graves PM, Burkot TR, Saul AJ, Hayes RJ, Carter R. Estimation of anopheline survival rate, vectorial capacity and mosquito infection probability from malaria vector infection rates in villages near Madang, Papua New Guinea. Journal of Applied Ecology. 1990;27:134–147.
- 54. Rodríguez MH, Bown DN, Arredondo-Jiménez JI, Villarreal C, Loyola EG, Frederickson CE. Gonotrophic cycle and survivorship of Anopheles albimanus (Diptera: Culicidae) in southern Mexico. Journal of Medical Entomology. 1992;29:395–399. pmid:1625288
- 55. Bekessy A, Molineaux L, Storey J. Estimation of incidence and recovery rates of Plasmodium falciparum parasitaemia from longitudinal data. Bulletin of the World Health Organization. 1976;54:685–693. pmid:800968
- 56. Molineaux L, Gramiccia G. The Garki project: research on the epidemiology and control of malaria in the Sudan savanna of West Africa. WHO Publications, Albany, NY, USA; 1980.
- 57. Collins WE, Jeffery GM. A retrospective examination of mosquito infection on humans infected with Plasmodium falciparum. American Journal of Tropical Medicine and Hygiene. 2003;68:366–371. pmid:12685646
- 58. Dye C, Hasibeder G. Population dynamics of mosquito-borne disease: effects of flies which bite some people more frequently than others. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1986;80:69–77. pmid:3727001
- 59. Jansen VAA, Lloyd AL. Local stability analysis of spatially homogeneous solutions of multi-patch systems. Journal of Mathematical Biology. 2000;41:232–252. pmid:11072757
- 60. Formont E, Pontier D, Langlais M. Disease propagation in connected host populations with density-dependent dynamics: the case of Feline Leukemia Virus. Journal of Theoretical Biology. 2003;223:465–475.
- 61. Smith DL, Battle KE, Hay SI, Barker CM, Scott TW. Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biology. 2007;5:e42. pmid:17311470
- 62. Gurarie D, Seto YW. Connectivity sustains disease transmission in environments with low potential for endemicity: modelling schistosomiasis with hydrologic and social connectivities. Journal of the Royal Society Interface. 2009;6:495–508.
- 63. Hess G. Disease in metapopulation models: implications for conservation. Ecology. 1996;77:1617–1632.
- 64. McCallum H, Dobson A. Disease and connectivity. In: Connectivity Conservation. Cambridge University Press; 2006. p. 479–501.
- 65. Smith D, McKenzie FE. Statics and dynamics of malaria infection in Anopheles mosquitoes. Malaria Journal. 2004;3(1):13. pmid:15180900
- 66. Dietz K. The estimation of the basic reproduction number for infectious diseases. Statistical Methods in Medical Research. 1993;2(1):23–41. pmid:8261248
- 67. Dietz K. Models for vector-borne parasitic diseases. Lecture Notes in Biomathematics. 1980;39:264–277.
- 68. Dietz K, Wernsdorfer W, McGregor I. Mathematical models for transmission and control of malaria. Malaria: principles and practice of malariology Volume 2. 1988;p. 1091–1133.
- 69. Smith DL, Dushoff J, Snow RW, Hay SI. The entomological inoculation rate and its relation to prevalence of Plasmodium falciparum infection in African children. Nature. 2005;438:492–495. pmid:16306991
- 70. Stoddard ST, Morrison AC, Vazquez-Prokopec GM, Soldan VP, Kochek TJ, Kitron U, et al. The role of human movement in the transmission of vector-borne pathogens. PLoS Neglected Tropical Diseases. 2009;3(7):e481. pmid:19621090
- 71. Bautista CT, Chan AS, Ryan JR, Calampa C, Roper MH, Hightower AW, et al. Epidemiology and spatial analysis of malaria in the Northern Peruvian Amazon. American Journal of Tropical Medicine and Hygiene. 2006;75:1216–1222. pmid:17172396
- 72. Ernst KC, Adoka SO, Kowuor DO, Wilson ML, John CC. Malaria hotspot areas in a highland Kenya site are consistent in epidemic and non-epidemic years and are associated with ecological factors. Malaria Journal. 2006;5:78. pmid:16970824
- 73. Eagle N, Pentland A, Lazer D. Inferring friendship network structure by using mobile phone data. Proceedings of the National Academy of Sciences. 2009;106:15274–15278.
- 74. Weslowski A, Eagle N, Noor AM, Snow RW, Buckee CO. Heterogeneous mobile phone ownership and usage patterns in Kenya. PLoS ONE. 2012;7(4):e35319.
- 75. Wesolowski A, Eagle N, Tatem AJ, Smith DL, Noor AM, Snow RW, et al. Quantifying the impact of human mobility on malaria. Science. 2012;338(6104):267–270. pmid:23066082