Computational and Mathematical Methods in Medicine
Volume 2018 (2018), Article ID 2575017, 5 pages
https://doi.org/10.1155/2018/2575017

Research Article

Clustered Breeding Sites: Shelters for Vector-Borne Diseases

J. C. A. Dias¹ and L. H. A. Monteiro1,2

¹Universidade Presbiteriana Mackenzie, PPGEEC, São Paulo, SP, Brazil
²Universidade de São Paulo, Escola Politécnica, São Paulo, SP, Brazil

Correspondence should be addressed to L. H. A. Monteiro; luizhenrique.monteiro@mackenzie.br

Received 4 March 2018; Accepted 3 June 2018; Published 9 July 2018

Academic Editor: Michele Migliore

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.

Abstract.

Here, the propagation of vector-borne 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 public-health 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 vector-borne pathogen transmission began to be mathematically analyzed in 1908 by Ross [6], who formulated the so-called “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 vector-borne 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 SIR-type 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 disease-prevention 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 three-dimensional torus is formed from this two-dimensional 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 th-column, if is even and if is odd.(iv)One-cluster distribution: in each cell of a region of the CA lattice, then ; outside this region, .(v)Two-clusters distribution: in each cell of two regions , then ; outside these regions, . These two regions are horizontally separated by cells.(vi)Four-clusters 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)Time-invariant 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.

Figure 1: Time evolution of (dotted line), (thick line), and (thin line) from , , and . Parameter values: , , , and , time-invariant uniform distribution with . In this case, the disease endemically persists at a constant level (around 0.1565, as shown in Table 1).

Table 1: Normalized infected population in steady state as a function of for six time-invariant distributions.

distribution

uniform
0.1565 0.0004
0.1667 0.0004
0.1712 0.0003

random
0.1511 0.0005
0.1626 0.0004
0.1679 0.0003

column
0.1480 0.0004
0.1607 0.0004
0.1667 0.0003

1-cluster
0.0455 0.0001
0.0458 0.0001
0.0459 0.0002

2-clusters
0.0462 0.0002
0.0466 0.0001
0.0466 0.0002

4-clusters
0.0476 0.0002
0.0480 0.0002
0.0482 0.0002

Table 1 shows the endemic steady-states 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 public-health interventions.

Table 2: Critical amount of vectors for six time-invariant distributions.

distribution

uniform
1200

random
1000

column
1000

1-cluster
300

2-clusters
300

4-clusters
300

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.

Table 3: Number of eradications in 100 simulations in function of for six seasonal distributions.

distribution

uniform
100
100

random
100
100

column
100
100

1-cluster
0
18

2-clusters
0
41

4-clusters
0
82

Figure 2: Parameter values: , , , , periodic uniform distribution with , , and (thus, oscillates between 4000 and 180 with period 30). In this case, the disease naturally disappears (that is, as the time passes).

Figure 3: Parameter values: , , , , periodic one-cluster distribution with , , and . In this case, the disease prevalence tends to a periodic oscillation.

4. Conclusions

The influence of spatiotemporal heterogeneities of vectors on the propagation of vector-borne 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 vector-borne diseases [7, 8, 11–13, 20]. Our simulations suggest that vector-control 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 vector-control programmes should be to eliminate the large clustered breeding sites of infected big cities, because these sites are shelters for vector-borne 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/2014-6).

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.

M. T. D. O. Mota, A. C. Terzian, M. L. C. R. Silva, C. Estofolete, and M. L. Nogueira, “Mosquito-transmitted viruses – the great Brazilian challenge,” Brazilian Journal of Microbiology, vol. 47, pp. 38–50, 2016.

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.

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.

C. Possas, R. M. Martins, R. L. de Oliveira, and A. Homma, “Urgent call for action: Avoiding spread and re-urbanisation of yellow fever in Brazil,” Memórias do Instituto Oswaldo Cruz, vol. 113, no. 1, pp. 1-2, 2018.

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 susceptible-infective-recovered model,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 83, no. 3, Article ID 037101, 4 pages, 2011.

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.

C. J. Dommar, R. Lowe, M. Robinson, and X. Rodó, “An agent-based model driven by tropical rainfall to understand the spatio-temporal heterogeneity of a chikungunya outbreak,” Acta Tropica, vol. 129, no. 1, pp. 61–73, 2014.

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.

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.

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.

L. B. Santos, M. C. Costa, S. T. Pinho et al., “Periodic forcing in a three-level cellular automata model for a vector-transmitted disease,” Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 80, no. 1, Article ID 016102, 9 pages, 2009.

M. K. Enduri and S. Jolad, “Dynamics of dengue disease with human and vector mobility,” Spatial and Spatio-Temporal Epidemiology, vol. 25, pp. 57–66, 2018.

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. 1-2, pp. 258–262, 2006.

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.

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.

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.

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.

P. E. Parham, J. Waldock, G. K. Christophides et al., “Climate, environmental and socio-economic change: weighing up the balance in vector-borne disease transmission,” Philosophical Transactions of the Royal Society B: Biological Sciences, vol. 370, no. 1665, Article ID 20130551, 17 pages, 2015.

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 mark-release-recapture in the village of Shauri Moyo in eastern Kenya.,” Journal of Medical Entomology, vol. 32, no. 1, pp. 27–33, 1995.

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.

L. Regis, A. M. Monteiro, M. A. V. De Melo-Santos 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.

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.

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.

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.