Research Article  Open Access
Clustered Breeding Sites: Shelters for VectorBorne Diseases
Abstract
Here, the propagation of vectorborne diseases is modeled by using a probabilistic cellular automaton. Numerical simulations considering distinct spatial distributions and time variations of the vector abundance are performed, in order to investigate their impacts on the number of infected individuals of the host population. The main conclusion is as follows: in the clustered distributions, the prevalence is lower, but the eradication is more difficult to be achieved, as compared to homogeneous distributions. This result can be relevant in the implementation of preventive surveillance measures.
1. Introduction
Contagious diseases transmitted by insects remain a serious publichealth problem in many countries. For instance, in the past few years, the Brazilian people have been infected by Chagas disease [1], chikungunya, dengue, zika [2], leishmaniasis [3], malaria [4], and yellow fever [5].
The vectorborne pathogen transmission began to be mathematically analyzed in 1908 by Ross [6], who formulated the socalled “mosquito theorem.” According to this theorem, malaria would be naturally eradicated from a region if the mosquito abundance was reduced below a critical value in such a region. In this seminal work, however, the spatial dimension of this region was not explicitly taken into account.
Theoretical investigations on the spread of vectorborne diseases based on cellular automaton (CA) have been carried out [7–14]. In CA models, the spatial features of the host and vector populations, such as mobility patterns and heterogeneities, can be naturally taken into consideration in the CA lattice and in the rules of state transition.
Here, in our CA model, each cell of the CA lattice is occupied by one individual of the host population and in each cell there may be an amount of vectors. At each time step , each individual is in one of three states: susceptible (), infected (), or recovered (). The state transitions between the time steps and of this SIRtype epidemic model are driven by probabilistic rules. The goal is to examine the influence on the infected host group of distinct spatial distributions and time variations of the vector abundance, by running computer simulations.
This paper is organized as follows. In Section 2, the CA model is described. In Section 3, the results obtained from numerical simulations are presented. In Section 4, the relevance of the results for diseaseprevention campaigns is discussed.
2. The CA Model
In our CA model, the host population lives in a square lattice with cells. To eliminate edge effects, the top edge of the lattice contacts the bottom edge and the left edge contacts the right edge (thus, a threedimensional torus is formed from this twodimensional lattice). Each cell is occupied by one individual and each individual is in contact with its eight surrounding neighbors, which is usually known as Moore neighborhood of unit radius [15]. Note that, due to the boundary conditions chosen for the CA lattice, all individuals have the same number of neighbors. In addition, suppose that the amount of vectors associated with the th cell is . Therefore, the total amount of vectors in the CA lattice is .
The time evolution of the proposed SIR model is ruled by the following set of probabilities of state transition (for similar models, see, for instance, [16, 17]). At each time step , the probability of a individual being infected is given by , in which is the sum of considering the eight surrounding neighbors and is number of neighbors. Note that if and/or ; hence, the transition between the time steps and cannot occur if there are not vectors and/or infected individuals in the neighborhood of such a individual. Note also that monotonously increases with and/or . For “high” values of , then .
An individual has probability per time step of becoming cured and permanently protected against the infection; if not cured, then this individual has probability per time step of dying (due to the disease). A individual has probability per time step of dying (for other causes). Assume that when and individuals perish, individuals replace them. Consequently, the total number of individuals remains constant, which is a convenient assumption for modeling host populations in which the deaths are balanced by the births. Notice that the probabilities , , and correspond to the state transitions , , and , respectively. The states of all individuals are simultaneously updated throughout a simulation [18].
In many models, the vector population is divided into infected and noninfected subgroups [7–9, 11–14, 19, 20]. In our model, this division is not made; only the vector abundance matters. The rationale for this simplifying assumption is that the higher the amount of transmitters, the higher the probability of an neighbor being bitten and indirectly infecting a individual.
In short, the parameters of the model are the probabilities , , and , the total number of individuals , the total amount of vectors (for instance, mosquitoes) , the spatial distribution, and the time variation of .
The following spatial distributions are considered in the simulations:(i)Uniform distribution: in each cell of the CA lattice, .(ii)Random distribution: in each cell of the CA lattice, there is a 50% chance of and a 50% chance of .(iii)Column distribution: in each cell of the thcolumn, if is even and if is odd.(iv)Onecluster distribution: in each cell of a region of the CA lattice, then ; outside this region, .(v)Twoclusters distribution: in each cell of two regions , then ; outside these regions, . These two regions are horizontally separated by cells.(vi)Fourclusters distribution: in each cell of four regions , then ; outside, . These clusters are horizontally separated by cells.
Note that, in all distributions, ; thus, the total amount of vectors is the same. Therefore, only the impact of the geographical spread of vectors is evaluated.
We also consider the following time dependencies of :(i)Timeinvariant function: the value of in each cell is not altered during the simulation; therefore, is kept fixed.(ii)Periodic function: oscillates between two numbers with period . In this case, is kept fixed for time steps, then it is reduced to (with ) for the next time steps, then it returns to the original value for the next time steps, and so on. Thus, the total amount of vectors periodically varies between and . This variation can be a consequence of seasonal oscillations of climatic variables, such as temperature and humidity [21].
Also, there is no migration of vectors on the lattice; that is, it is supposed that the vectors can move only a short distance from their breeding sites [22, 23].
3. Simulation Results
Computer simulations were performed by taking: , , , (hence, ), and , , and (therefore, , 6000, and 8000, respectively). When time variation in is considered, , , and or 30. Simulations with other parameter values were performed, but the results were qualitatively the same as reported below. In all simulations, the initial condition (at ) is , , and . The asymptotical solution, however, does not depend on this starting point (that is, the attractor is globally asymptotically stable). Figure 1 illustrates a simulation with 200 time steps in which an endemic solution is reached. In this figure, is represented by a dotted line, by a thick line, and by a thin line.

Table 1 shows the endemic steadystates reached for = constant. For each , the average value of obtained in the last 100 time steps is given. Observe that the percentage of infected individuals grows with , as intuitively expected. For the clustered distributions, this percentage increases with the number of clusters. The values found for uniform, random, and column distributions are similar to each other. The values found for the three clustered distributions are also similar to each other, but they are smaller than those found for the three homogeneous distributions. Obviously, the number of infected individuals can be increased, for instance, by reducing and/or [18].
Inspired by the mosquito theorem, another set of simulations was run to determine the critical number of vectors below which the disease disappears. For values of below the ones shown in Table 2, the infection was naturally eradicated in 10 consecutive simulations. These critical values were numerically found by varying with a step size of 100. Observe that these two tables lead to a surprising conclusion: the prevalence in clustered distributions is lower than in homogeneous distributions; however, clustered distributions impair disease eradication. That is, clustered distributions require a greater effort to eliminate the disease than homogeneous distributions, in which the prevalence is higher! This result should be taken into account in the planning of publichealth interventions.

Table 3 presents the number of eradications in 100 simulations for seasonal variation of the amount of vectors. In these simulations, periodically oscillates between 4000 and 180 with period or 30. Figures 2 and 3 illustrate eradication and persistence of the disease, respectively. These simulations show that the more the vectors are clustered, the lower the number of cases in which the disease disappears. Observe that, for , the disease endemically persists only in the clustered distributions (it is always eradicated in the three homogeneous distributions); for , the disease is always eliminated in the homogeneous distributions and, for the clustered distributions, the number of eradications grows with the number of clusters.

4. Conclusions
The influence of spatiotemporal heterogeneities of vectors on the propagation of vectorborne diseases has been experimentally [24, 25] and theoretically [8, 26] analyzed. Here, our simplistic epidemic model reveals that the more the vectors are clustered in a given region, the lower is the prevalence; however, the greater must be the effort to eradicate the disease from such a region.
Reducing the vector population is one of the key methods employed to control the transmission of vectorborne diseases [7, 8, 11–13, 20]. Our simulations suggest that vectorcontrol programmes should take into consideration the spatial distribution of vectors and not only the disease prevalence and/or the vector population size. A decrease in the number of infections can also be achieved by restricting the mobility of the host population in the infected region [9, 11, 14, 19], but this surveillance measure may be difficult to implement in practice.
Big cities are shelters for directly transmitted infections [16, 27]. In agreement with other studies [11, 24], this work proposes that the main focus of vectorcontrol programmes should be to eliminate the large clustered breeding sites of infected big cities, because these sites are shelters for vectorborne diseases.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
L. H. A. Monteiro is partially supported by CNPq (under Grant no. 305827/20146).
References
 E. C. Freitas, M. de Fátima Oliveira, A. S. O. de Barros Vasconcelos et al., “Analysis of the seroprevalence of and factors associated with chagas disease in an endemic area in Northeastern Brazil,” Journal of the Brazilian Society of Tropical Medicine, vol. 50, no. 1, pp. 44–51, 2017. View at: Publisher Site  Google Scholar
 M. T. D. O. Mota, A. C. Terzian, M. L. C. R. Silva, C. Estofolete, and M. L. Nogueira, “Mosquitotransmitted viruses – the great Brazilian challenge,” Brazilian Journal of Microbiology, vol. 47, pp. 38–50, 2016. View at: Publisher Site  Google Scholar
 I. P. S. F. de Carvalho, H. M. Peixoto, G. A. S. Romero, and M. R. F. de Oliveira, “Cost of visceral leishmaniasis care in Brazil,” Tropical Medicine & International Health, vol. 22, no. 2, pp. 1579–1589, 2017. View at: Google Scholar
 J. Recht, A. M. Siqueira, W. M. Monteiro, S. M. Herrera, S. Herrera, and M. V. G. Lacerda, “Malaria in Brazil, Colombia, Peru and Venezuela: Current challenges in malaria control and elimination,” Malaria Journal, vol. 16, no. 1, Article ID 273, 18 pages, 2017. View at: Google Scholar
 C. Possas, R. M. Martins, R. L. de Oliveira, and A. Homma, “Urgent call for action: Avoiding spread and reurbanisation of yellow fever in Brazil,” Memórias do Instituto Oswaldo Cruz, vol. 113, no. 1, pp. 12, 2018. View at: Google Scholar
 R. Ross, Report on the Prevention of Malaria in Mauritius, Waterlow and Sons, London, UK, 1908.
 T. Botari, S. G. Alves, and E. D. Leonel, “Explaining the high number of infected people by dengue in Rio de Janeiro in 2008 using a susceptibleinfectiverecovered model,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 83, no. 3, Article ID 037101, 4 pages, 2011. View at: Google Scholar
 B. Cissé, S. El Yacoubi, and S. Gourbiére, “A cellular automaton model for the transmission of Chagas disease in heterogeneous landscape and host community,” Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 40, no. 2, pp. 782–794, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 C. J. Dommar, R. Lowe, M. Robinson, and X. Rodó, “An agentbased model driven by tropical rainfall to understand the spatiotemporal heterogeneity of a chikungunya outbreak,” Acta Tropica, vol. 129, no. 1, pp. 61–73, 2014. View at: Publisher Site  Google Scholar
 D. O. Gerardi and L. H. A. Monteiro, “System identification and prediction of dengue fever incidence in Rio de Janeiro,” Mathematical Problems in Engineering, vol. 2011, Article ID 720304, 13 pages, 2011. View at: Publisher Site  Google Scholar
 L. C. de Castro Medeiros, C. A. R. Castilho, C. Braga, W. V. de Souza, L. Regis, and A. M. V. Monteiro, “Modeling the dynamic transmission of dengue fever: Investigating disease persistence,” PLOS Neglected Tropical Diseases, vol. 5, no. 1, Article ID e942, 14 pages, 2011. View at: Publisher Site  Google Scholar
 F. M. Pereira and P. H. T. Schimit, “Dengue fever spreading based on probabilistic cellular automata with two lattices,” Physica A: Statistical Mechanics and its Applications, vol. 499, pp. 75–87, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 L. B. Santos, M. C. Costa, S. T. Pinho et al., “Periodic forcing in a threelevel cellular automata model for a vectortransmitted disease,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 80, no. 1, Article ID 016102, 9 pages, 2009. View at: Publisher Site  Google Scholar
 M. K. Enduri and S. Jolad, “Dynamics of dengue disease with human and vector mobility,” Spatial and SpatioTemporal Epidemiology, vol. 25, pp. 57–66, 2018. View at: Publisher Site  Google Scholar
 S. Wolfram, Cellular Automata and Complexity: Collected Papers, Westview Press, Boulder, USA, 1994.
 L. H. A. Monteiro, H. D. B. Chimara, and J. G. C. Berlinck, “Big cities: Shelters for contagious diseases,” Ecological Modelling, vol. 197, no. 12, pp. 258–262, 2006. View at: Publisher Site  Google Scholar
 P. H. T. Schimit and L. H. A. Monteiro, “On the basic reproduction number and the topological properties of the contact network: An epidemiological study in mainly locally connected cellular automata,” Ecological Modelling, vol. 220, no. 7, pp. 1034–1042, 2009. View at: Publisher Site  Google Scholar
 J. C. A. Dias and L. H. A. Monteiro, “Investigando a influência da distribuição espacial de mosquitos na propagação de doenças via autômato celular,” in Proceedings of the XVII Workshop de Informática Médica, vol. 1, pp. 1–7, São Paulo, Brazil, 2017. View at: Google Scholar
 G. R. Phaijoo and D. B. Gurung, “Modeling Impact of Temperature and Human Movement on the Persistence of Dengue Disease,” Computational and Mathematical Methods in Medicine, vol. 2017, Article ID 1747134, 9 pages, 2017. View at: Publisher Site  Google Scholar
 Y. Cheng, X. Wang, Q. Pan, and M. He, “Modeling the parasitic filariasis spread by mosquito in periodic environment,” Computational and Mathematical Methods in Medicine, vol. 2017, Article ID 4567452, 10 pages, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 P. E. Parham, J. Waldock, G. K. Christophides et al., “Climate, environmental and socioeconomic change: weighing up the balance in vectorborne disease transmission,” Philosophical Transactions of the Royal Society B: Biological Sciences, vol. 370, no. 1665, Article ID 20130551, 17 pages, 2015. View at: Publisher Site  Google Scholar
 M. Trpis, W. Häusermann, and G. B. Craig Jr., “Estimates of population size, dispersal, and longevity of domestic Aedes aegypti aegypti (Diptera: Culicidae) by markreleaserecapture in the village of Shauri Moyo in eastern Kenya.,” Journal of Medical Entomology, vol. 32, no. 1, pp. 27–33, 1995. View at: Publisher Site  Google Scholar
 M. W. Service, “Mosquito (Diptera: Culicidae) Dispersal  The Long and Short of It,” Journal of Medical Entomology, vol. 34, no. 6, pp. 579–588, 1997. View at: Publisher Site  Google Scholar
 L. Regis, A. M. Monteiro, M. A. V. De MeloSantos et al., “Developing new approaches for detecting and preventing Aedes aegypti population outbreaks: Basis for surveillance, alert and control system,” Memórias do Instituto Oswaldo Cruz, vol. 103, no. 1, pp. 50–59, 2008. View at: Publisher Site  Google Scholar
 S. Wanji, F. F. Mafo, N. Tendongfor et al., “Spatial distribution, environmental and physicochemical characterization of Anopheles breeding sites in the Mount Cameroon region,” Journal of Vector Borne Diseases, vol. 46, no. 1, pp. 75–80, 2009. View at: Google Scholar
 A. M. Lutambi, M. A. Penny, T. Smith, and N. Chitnis, “Mathematical modelling of mosquito dispersal in a heterogeneous environment,” Mathematical Biosciences, vol. 241, no. 2, pp. 198–216, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 M. S. Bartlett, “Measles Periodicity and Community Size,” Journal of the Royal Statistical Society. Series A (General), vol. 120, no. 1, pp. 48–70, 1957. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 J. C. A. Dias and L. H. A. Monteiro. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.