 Research
 Open access
 Published:
The influence of gender and temephos exposure on community participation in dengue prevention: a compartmental mathematical model
BMC Infectious Diseases volume 24, Article number: 463 (2024)
Abstract
Background
The use of temephos, the most common intervention for the chemical control of Aedes aegypti over the last half century, has disappointing results in control of the infection. The footprint of Aedes and the diseases it carries have spread relentlessly despite massive volumes of temephos. Recent advances in community participation show this might be more effective and sustainable for the control of the dengue vector.
Methods
Using data from the Camino Verde cluster randomized controlled trial, a compartmental mathematical model examines the dynamics of dengue infection with different levels of community participation, taking account of gender of respondent and exposure to temephos.
Results
Simulation of dengue endemicity showed community participation affected the basic reproductive number of infected people. The greatest shortterm effect, in terms of people infected with the virus, was the combination of temephos intervention and community participation. There was no evidence of a protective effect of temephos 220 days after the onset of the spread of dengue.
Conclusions
Male responses about community participation did not significantly affect modelled numbers of infected people and infectious mosquitoes. Our model suggests that, in the long term, community participation alone may have the best results. Adding temephos to community participation does not improve the effect of community participation alone.
Background
Control of the Aedes aegypti is the cornerstone of public health programs to prevent transmission [1, 2] of at least four systemic viral diseases carried by the mosquito, dengue fever [3], zika [4], chikungunya [5] and yellow fever [6]. For most of the last century, temephos has been at the centre of chemical control of the vector [7, 8]. Limited effectiveness may be related to the small window of chemical activity and increasing resistance to temephos [8, 9].
Several studies show low effectiveness of the pesticide on the dengue transmission, for reasons including larval resistance to temephos [8, 10, 11]. Aside from its inconclusive impact on dengue transmission, widespread use of temephos may produce a sense of “being protected” against dengue that could discourage participation in prevention activities [12]. Recent approaches to community participation seem to be more sustainable and effective [13,14,15]. A randomised controlled trial in Mexico and Nicaragua demonstrated that community participation adds effectiveness to government dengue control programs [12]. The same effect was not achieved across all intervention sites in that trial, however, and gender issues and social capital might influence the impact of community participation [16].
Many mathematical models of infection transmission rely on the 1927 proposal of Kermack and McKendrick. Summarised by the acronym SIR, this divides the population in three compartments: susceptible, without the infection but on risk of being infected; infected, with the disease and be able to spread it; recovered, those who recovers from the disease and have permanent immunity against the disease. In the case of dengue, several SIR models attempted to explain the dynamics of different dengue virus serotypes [17,18,19,20] and to analyse the effect of vector control strategies [21,22,23]. Compartmental models usually precede optimum control models, [24,25,26] used to analyse combination of control mechanisms for Aedes aegypti, to quantify the impact of different strategies while minimize or maximize the effect of each intervention in individual or in combination [27].
A previous study showed gender of respondents affected the results chain for dengue prevention behaviour [16] but we could not find published of mathematical models of gender and temephos on community participation in dengue vector control. We present a compartmental mathematical model that helps to explain the dynamics of the spread of virus dengue infection, considering different levels of community participation by gender and temephos exposure.
Methods
Formulation of the model
We used data from impact measurement survey of the Mexico arm of the Camino Verde [12] trial to explore community participation according to male and female respondents, and to explore the linkages between participation and temephos exposure. For this proposal, we understand community participation as specific groups with shared needs living in a defined geographic area actively pursue identification of their needs, take decisions and establish mechanisms to meet these needs [28]. The effect of community participation in our approach is the difference in vector indices between the intervention and control groups in the Camino Verde trial. We adapted the modelling approach of several authors [17, 18, 22, 25, 29, 30]. In formulation of the model, we assumed the infection is produced by only one serotype of dengue virus and that the human and Aedes aegypti populations are divided into compartments.
The compartments in which the populations are divided are susceptible, infectious and recovered. The dynamics for movement between these compartments are the following:

i)
The human population is assumed to be constant with size equal to N_{h}, with birth and death rate constant equal to μ_{h}. The population is divided in susceptible, infected and recovered classes, denoted by \({\overline{S} }_{h}\), \({\overline{I} }_{h}\), and \({\overline{R} }_{h}\), respectively, with \({\overline{S} }_{h}+{\overline{I} }_{h}+{\overline{R} }_{h}={N}_{h}\). Dengue transmission is sustained by the flows between humans and mosquito compartments. The human susceptible population is decreased following the infection force rate, which can be acquired via effective contact with an exposed or infective vector at a rate
$${\tau }_{h}=\alpha {\beta }_{vh}\left(1C\left(t\right)\right)\frac{{\overline{I} }_{v}}{{N}_{h}},$$(1)
where α is the biting rate per susceptible vector, β_{vh} is the transmission probability from an infective vector (\({\overline{I} }_{v}\)) to a susceptible human (\({\overline{S} }_{h}\)), and C(t) is the effort made by men and women in the control of Aedes aegypti through community participation. Infectious humans recover at a constant rate ρ.

ii)
Following Moulay [31], a stage structured model describes the vector population dynamics, which consists of three stages: egg (E), larva (L) and pupa (P). We differentiate between eggs (E) and aquatic stages (L and P), because they respond differently to the control measures. In order to study the influence of the use of temephos and community participation, the breeding sites of Aedes aegypti were divided in targeted by themephos household containers (thc)  tanks, drums, water tanks and cisterns  and nontarget household containers (nthc) that periodically renew the water or containers that can be considered as garbage to eliminate them. We assume that the average number of eggs is proportional to the number of female mosquitoes that will result, and that the availability of nutrients and space is different in each type of containers (thc and nthc); similarly, the number of eggs is proportional to the number of larvae and pupae.
The adult female mosquitoes M is divided in two compartments: susceptible and infective, whose are denoted by \({\overline{S} }_{v}\) and \({\overline{I} }_{v}\). The vector susceptible population is decreased following infection, which can be acquired via effective contact with an infected human at a rate
where β_{hv} is the transmission probability from an infectious human (\({\overline{I} }_{h}\)) to a susceptible vector (\({\overline{S} }_{v}\)).
Our model represents dengue transmission dynamics under the same assumptions as the Camino Verde trial: human and mosquito populations are mixed homogeneously, meaning each mosquito has the same probability of biting any individual in the population; since the outbreaks are relatively short, the population can be considered constant; reducing the mosquito population will reduce dengue cases; a single serotype is implicated in the outbreaks (allowing modelling of reduction of the susceptible population).
The entomological parameters of vector dengue, for each type of container i = {thc,nthc}, are the per capita oviposition rate θ_{i}, which is the average number of eggs by female mosquito; the transition rate from egg to larva ε_{i}; the transition rates from larva to pupa λ_{i} and pupa to adult γ_{i}; mortality rates of eggs, larvas, pupas and mosquitoes, μ_{E}, μ_{L}, μ_{P}, μ_{v}, respectively.
From the above assumptions, the model is given by the following system of differential equations:
where τ_{h} and τ_{v} are given by (1) and (2), respectively; A is the effect produced by temephos in the target household containers. A schematic of the model is shown in Fig. 1. The parameters are all strict positive constants described in Table 1.
Since \(M={\overline{S} }_{v}+{\overline{I} }_{v}\), we see from the last eight equations of (3) that mosquito dynamics from egg stage to adult form is given by
We consider that the maturation rates from pupae to adult in the different type of container, \({\gamma }_{thc}\left(t\right)\) and \({\gamma }_{nthc}\left(t\right)\), are modelled for a sinusoidal function with 1year period. This assumption allows us to simulate the peaks dengue seasons observed in [16] with the periodic function proposed in [33]
where \({\gamma }_{0}\) average reproduction maturation rate from pupae to adult in the different type of container, \(\zeta\) is the amplitude of the seasonal variation in the reproduction rate of vectors, with \(0\le \zeta <1\), and \(\varphi\) is the phase angle to adjust the peak season for mosquitoes.
Community participation effect
Recognizing that community participation intervenes in vector control differently than temephos [7, 36, 37], the model assumes that community participation affects: i) all stages of the mosquito life cycle, ii) both types of containers, with and without temephos, iii) reduces the load of infection by reducing the probability of humanmosquito transmission; and iv) that temephos only acts in the larval stage.
Based on the results obtained by Andersson et al. [12], we believe community participation affects infection risk slowly at the beginning of the intervention, as a critical mass knowledge and interest accumulates. Once people gain experience with their selected control strategies, participation has an increasing effect until it peaks, after which smaller changes maintain the impact. For that reason, we used a nonstandard logistic function to analyze this pattern on reduction of pupae per person index.
One way to model different collaborative scenarios by sex for the effectiveness of community participation would be given by C(t), the sum of female and male contribution in the intervention, modelled by nonstandard logistic functions for each group, where C_{0} is the preexisting effectiveness of community participation. Community participation is not a Boolean on/off switch. In its initial stages, it might be quite clumsy and inefficient. As participants experience the impact of their efforts, this gains momentum and effectiveness. r is the incremental rate of community participation effectiveness, modelling the idea that experiencing impact reinforces selfconfidence and motivation [12]. The mobilisation strategy might show little effectiveness at the beginning of community participation, then increase in effectiveness is a function of r and finally stability of effectiveness after some time. If k_{max} is the total maximum capacity of community participation effectiveness for the control of the dengue vector, then the women’s contribution capacity, k_{F}, is a proportion to this quantity, where p_{F} and p_{M}are the proportions of participation of women and men in control mosquito activities While participation of both women and men is important, several studies have reported greater involvement of women than men in community dengue vector control [16, 38, 39]. The general case is therefore p_{F} > p_{M} (see S1 Appendix).
C(t) represents the maximum effectiveness of community participation at time t ∈ [0,T] of 13%, according to superior value of 95%CI for risk difference of pupae per person index calculated in Nicaragua and Mexico [12]. To achieve this risk difference reduction, about 25% of women reported participation in control vector activities, while only 2% of men participated in such activities (see S2 Appendix).
Temephos efficacy depends of different factors related to temephos doses, water turnover rate, type of water, and environmental factors around water storage such as organic debris presence, temperature and exposure to sunlight [8]. In studies reporting risk difference in the pupa per person index between 1 and 5%, the residual effect of temephos lasted from two and three months [34, 35]. Our model assume a constant temephos effect of 3%, as an average of the values found in the literature, which we believe balances out the effectiveness of the intervention during different seasons of the year.
Our model considers the effects of different levels of male and female participation. In this way the nonstandard logistic function, C(t), provides evidence of the effectiveness of community participation, as a strategy for the control of Aedes aegypti, considering the separate contributions of man and women.
Analysis of the model
We define the net reproductive number of the mosquito population and then focus on the dynamics of dengue transmission applying a formula for the basic reproductive number and an estimate of the endemic equilibrium point.
Mosquito dynamics
We analyse the mosquito population dynamics given by system (4). Considering constant values of the parameters, system (4) has two equilibrium points, the mosquitoesfree state E_{0} = (0,0,0,0,0,0,0) which correspond to the trivial equilibrium, and the state characterized by the presence of mosquitoes denoted by
With the elements of E_{1} we obtained the basic offspring number, denote by R_{M} (details in S3 Appendix).
The basic reproductive number
By analysing the dengue dynamics, we will assume that R_{M} > 1, since otherwise the mosquito population is zero, and therefore there is no disease. For R_{M} > 1, solutions approach asymptotically the equilibrium E_{1}, and therefore we can assume that the mosquito population has already reached its equilibrium, and the total population of adult mosquitos is constant and equal to \(M^{\ast}\). Considering that \({S}_{v}=\frac{{\overline{S} }_{v}}{{M}^{*}}\), \({I}_{v}=\frac{{\overline{I} }_{v}}{{M}^{*}}\) and S_{v} = 1I_{v}, model (3) is equivalent to the system of differential equations for the proportions
Our model hinges on the effective reproduction number of the disease, R_{eff}, the number of people in a population who can be infected by an individual at any specific time (14). R_{eff} is a fraction of the basic reproduction number R_{0}, the average number of secondary cases that one case can produce if introduced into a susceptible population of humans and mosquitoes. If R_{0} < 1, less than one secondary case will arise from a primary case and the disease will fade out. If R_{0} > 1 an outbreak will start.
The model has two equilibria, the diseasefree state P_{0} = (1,0,0), and the denguepresent state denoted by
where
\({I}_{v}^{*}=\frac{{\mu }_{h}{N}_{h}\left(C+{\mu }_{v}\right)\left(\rho +{\mu }_{h}\right)\left({R}_{eff}1\right)}{\alpha {\beta }_{vh}\left(1C\right){M}^{*}\left[\left(C+{\mu }_{v}\right)\left(\rho +{\mu }_{h}\right)+{\mu }_{h}\alpha {\beta }_{hv}\left(1C\right)\right]}\) and R_{eff} is obtained from the follow lemma
Lemma 3.2.1 Considered that the maturation rates from pupae to adult in the different type of container, γ_{thc}(t) and γ_{nthc}(t), and the effectiveness of community participation, C(t), are continuos, bounded, positive, periodic (in the case of γ_{thc}(t)) and not identically zero functions of time. The average values \({\widehat{\gamma }}_{thc}\left(t\right)\) and \(\widehat{C}(t)\) taken over a cycle are \({\widehat{\gamma }}_{thc}\left(t\right)=\frac{1}{t}{\int }_{0}^{t}{\gamma }_{thc}\left(x\right)dx\) and \(\widehat{C}\left(t\right)=\frac{1}{t}{\int }_{0}^{t}C\left(x\right)dx\); if γ_{thc}(t) and C(t) are replaced by that average values \({\widehat{\gamma }}_{thc}\left(t\right)\) and \(\widehat{C}(t)\) then effective reproduction number of the disease is
where \({N}_{v}\) is the number of mosquitoes in the absence of intervention, \({R}_{0}\) is the basic reproduction number and \({\widehat{M}}^{*}\) is obtained according [40,41,42] to replaced \({\widehat{\gamma }}_{thc}\) and Ĉ(t). The average values \({\widehat{\gamma }}_{thc}\left(t\right)\) and Ĉ(t) are
Proof From (5) is obtained that \({\gamma }_{thc}\left(t\right)={\gamma }_{0}\left(1+\zeta {\text{cos}}\left(\frac{2\pi }{365}t+\varphi \right)\right)\), then by definition of \({\widehat{\gamma }}_{thc}\left(t\right)\) is obtained that \({\widehat{\gamma }}_{thc}\left(t\right)=\frac{1}{t}{\int }_{0}^{t}{\gamma }_{0}\left(1+\zeta {\text{cos}}\left(\frac{2\pi }{365}t+\varphi \right)\right)dx\). Resolving the integral is obtained that \({\int }_{0}^{t}{\gamma }_{0}\left(1+\zeta {\text{cos}}\left(\frac{2\pi }{365}t+\varphi \right)\right)dx=\left(t{\gamma }_{0}+\frac{365}{2\pi }\zeta {\gamma }_{0}{\text{sin}}\left(\varphi +\frac{2\pi }{365}t\right)\frac{365}{2\pi }\zeta {\gamma }_{0}{\text{sin}}\left(\varphi \right)\right)\), so
If C(t) = C_{F}(t) + C_{M}(t), then by definition Ĉ(t), and considering (6) and (7), getting \(\widehat{C}\left(t\right)=\frac{1}{t}{\int }_{0}^{t}\left[\frac{{k}_{F}{C}_{0}{e}^{rt}}{{k}_{F}+{C}_{0}\left({e}^{rt}1\right)}+\frac{\left({{k}_{{\text{max}}}k}_{F}\right){C}_{0}{e}^{rt}}{\left({{k}_{{\text{max}}}k}_{F}\right)+{C}_{0}\left({e}^{rt}1\right)}\right]dx\). Resolving and simplifying \({\int }_{0}^{t}\left(\frac{{k}_{F}{C}_{0}{e}^{rt}}{{k}_{F}+{C}_{0}\left({e}^{rt}1\right)}+\frac{\left({{k}_{{\text{max}}}k}_{F}\right){C}_{0}{e}^{rt}}{\left({{k}_{{\text{max}}}k}_{F}\right)+{C}_{0}\left({e}^{rt}1\right)}\right)dx=\left(\frac{{k}_{F}}{r}{\text{ln}}\left(1+\frac{{C}_{0}\left({e}^{rt}1\right)}{{k}_{F}}\right)+\frac{\left({k}_{{\text{max}}}{k}_{F}\right)}{r}{\text{ln}}\left(1+\frac{{C}_{0}\left({e}^{rt}1\right)}{\left({k}_{{\text{max}}}{k}_{F}\right)}\right)\right)\), from where
Estimation on the human population for Acapulco city
In the last century, the Acapulco’s population size had a nonlinear increase [43] though, from the year 2000, the increase was more stable. We therefore adjusted a nonlinear regression model to the data of the Acapulco population between the years 1910–2020. The nonlinear function considered was as follows:
where y is the Acapulco population, x is time in years, u_{0} is the initial population, K the saturation capacity and R the population growth rate. A value of 834 634 was obtained as a estimator of K with a 95% CI [803 423.7, 865 844.3], used as the population size of Acapulco N_{h}, and is assumed approximately constant in the next five years.
Results
Simulation of effectiveness of community participation
Figure 2 shows the modelling of total effectiveness of community participation. This considers a strategy aimed to trigger community participation achieves the greatest effectiveness in one and a half years (dotted line in panel ‘a’ of Fig. 2), and an estimator for incremental rate of community participation effectiveness is r = 0.015. Panel ‘b’ shows the community participation effectiveness divided by gender, dotted curves is the contribution of men for different r values.
The basic reproductive number for dengue cases is significantly affected by community participation.
Taking the maximum values of E_{thc}, E_{nthc}, L_{thc} and L_{nthc} when the populations of eggs and larvae stabilize over time, then when no intervention is implemented on Aedes aegyti control, the obtained value is R_{eff} = R_{0} = 3.3. If an intervention with temephos is implemented, considering a constant effectiveness of 3%, the obtained value is R_{eff} = 3.1, both scenarios are auspicious for the occurrence of a dengue outbreak. When social mobilization achieves an effectiveness of 8.5%, as result of community participation added to temephos intervention, the value obtained is R_{eff} < 1 (Fig. 3). Under these circumstances, control of dengue cases is achieved.
Simulation of dynamic of dengue
The simulation was carried out using the values of Table 1 and the initial values \({\overline{S} }_{h}\left(0\right)=833799\), \({\overline{I} }_{h}\left(0\right)=835\), R_{h}(0) = 0 in addition E_{1} = (974666, 866743, 866743, 317083, 317083, 487821, 487821) give by (9). We implemented the simulation using the “deSolve” package of the R [44].
Figure 4 shows the proportions of infected humans and infectious vectors from the simulation model (3). Up to forty days after the start of the dengue outbreak, the different interventions for Aedes aegypti control produce the same results. There is no evidence of protective effect of temephos intervention 220 days after outbreak initiation: the greater effect in short term  in terms of people infected with dengue  was the combination of temephos intervention and community participation (Fig. 4).
Simulation considering different levels of community participation by gender
Figure 5 shows different scenarios of the proportions of infected people and infectious mosquitoes. The simulation showed how increasing the contribution of men in community participation, to its initial participation of 2% and gradually increasing up to 25%, decreases more rapidly the proportion of infected people and infectious mosquitoes after day 120, it reaches its maximum difference 250 days after the start of the outbreak; in this period, there was no significant difference in the proportions of infected people and infectious mosquitoes, even when the percentage of participation of men increased by more than 10% (Fig. 5).
Discussion
Our model is well defined biologically and mathematically. It suggests that, in the long term, community participation on its own produces the strongest protection. Adding temephos to community participation does not improve its effect, possibly because temephos might produce a false feeling protection [15]. Our compartmental model is compatible with the association reported by Andersson et al. [12] between use of temephos and increased risk of dengue virus infection. The Camino Verde trial showed a 13% reduction in the pupal index per person with community engagement, with higher levels of protection in clusters where women engaged in vector control activities.
Our simulation of dengue endemicity reflects a basic reproductive number significantly affected by community participation. An initial reduction associated with temephos use can be explained by the combination of temephos intervention and community participation (Fig. 4). There was no protective effect of temephos intervention 220 days after the outbreak began, this agrees with the Nicaraguan report, where exposure to temephos was associated with higher entomological indices over time [40]. This was explained by several factors including the ecological adaptability of the vector, resistance of Aedes to the pesticide, operational deficiencies of the vector control program (lower quantities used at greater intervals in only a sample of reservoirs), a decrease in community motivation as a consequence of a false sense of protection when temephos is present in their water. Our simulation showed increasing the contribution of men in community participation by more than 10% might decrease slightly the proportions of infected people and infectious mosquitoes (Fig. 5).
Community participation affects the risk of infection in a variable way. Initially there may be little change. As a critical mass of interest, knowledge and confidence accumulates, the quality of participation changes and its effect is more noticeable. After people gain proficiency with a particular control strategy, implementing it fully, the biological potential of the specific activity to control reproduction reaches a peak, before adaptation of the vector mosquito.
We did simulations that calculate R_{0} for the different strategies, taking the maximum values of E_{thc}, E_{nthc}, L_{thc} and L_{nthc} when the populations of eggs and larvae stabilize over time. This generates an effectiveness of community participation of 8.5%, translating as a reduction of R_{0} < 1. The model assumes a constant biting rate but reduced humanvector contact, with a transmission rate \(\alpha {\beta }_{vh}\left(1C\left(t\right)\right){\overline{S} }_{h}\frac{{\overline{I} }_{v}}{{N}_{h}}\), where \(1C\left(t\right)\) describes the reduction in contacts between infected mosquitos and susceptible humans, reflecting the reduced mosquito population. This approach has been used in other models of vectorborne diseases [25, 26, 45].
Mathematical models have proved useful to understand the transmission of dengue and to help plan infection control strategies [17,18,19,20,21,22,23,24,25,26,27,28,29]. Optimizing the intervention mix of control strategies is a relatively recent trend in mathematical modelling. Our compartmental model, without identifying an optimum control approach, helps to explain transmission dynamics according to the reduction of the relative risk for the proportion of infected people and the rate of pupae per person.
Limitations
Our model does not consider the variations in the environment (temperature, humidity and height above sea level), and community location (urban, rural or suburbs) which are important if partially studied determinants of the dengue vector density. [29, 46] Most of Acapulco is less than 300 m above sea level, the environmental temperature is warm most of the year, and humidity differs depending with the rainy season. Our model took the variation in humidity into account by a sinusoidal function for estimating the pupal population. Other limitation is the lack of another compartment, in this specific case, the estimated population into latent period of the disease, which is not included in the SIR model.
Dengue mortality occurs mainly in children and the elderly and fortunately there is a low mortality rate due to severe dengue. We did not include the effect of mortality on transmission. We based simulations on the Acapulco general population, so not including mortality indicators will have little effect on the model.
Our lack of longitudinal data on infected individuals in intervention and control sites limits our ability to inform our model with actual dengue cases in the trial setting. We used historical dengue case reports from Acapulco. Both 2013 and 2015 reported two outbreaks in the year, the second with a higher incidence than the first. Shepard et al. [47] suggested an expansion factor of 1.4–3.3 to allow for the wellrecognized underreporting. The results of our simulation without the intervention thus coincide with 2–3 times the average number of cases reported during the peaks of the outbreaks in Acapulco in the years preceding the trial, 7 307 in 2013 and 5 586 in 2015.
We considered the temephos effect without any variation in effectiveness or use between households. Temephos effectiveness varies dramatically with the dose present and amount of water in the container, mosquito breeding density, time temephos is left in containers, rain season and other factors. Our model used average temephos effectiveness, possibly overestimating its impact without affecting the results of community participation. It was possible to consider the time dependence of the effect of temephos, as the Fig. 4 shows, after 220 days effectiveness of temephos dropped drastically.
We did not model different engagement dynamics for men and women [16]. We are aware of different roles in participatory vector control – men in Mexico would be more concerned with outdoor or community actions, with women more involved in household activities. We modelled the effect as additive, although other interactions (synergism or proportional functions) may have played a role. This aspect of our model could hide the real effect of gender in community participation not least, as pointed out by Andersson and colleagues [12], the possibility that the way Camino Verde engaged women acted as a disincentive to male engagement. There could be unidentified ripple effects of more equitable engagement, including better family communication, reduction of the women workload.
Community participation adds value to conventional control strategies for Aedes aegypti control for dengue fever. Our model suggests that, in the long term, use of temephos alone could fail to achieve Aedes aegypti control, whereas community participation might have a more sustained effect.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Abbreviations
 S1:

Supplement 1
 S2:

Supplement 2
 S3 :

Supplement 3
 SIR:

Susceptibleinfectedrecovered model
References
World Health Organization. Dengue and severe dengue. http://www.who.int/mediacentre/factsheets/fs117/en/. Accessed 15 June 2022.
Pan American Health Organization, World Health Organization. Scientists studying intensified vector control measures to combat Zika, dengue and chikungunya in the Americas. http://www.paho.org/hq/index.php?option=com_content&view=article&id=11780. Accessed 15 June 2022.
Bhatt S, et al. The global distribution and burden of dengue. Nature. 2013;496:504–17.
Hayes EB. Zika virus outside Africa. Emerg Infect Dis. 2009;15:1347–50.
Montero A. Chikungunya fever – a new global threat. Med Clin. 2015;145:118–23.
Gill GV, Beeching N. Tropical medicine: lecture notes. 7th ed. Oxford: Wiley; 2014. p. 283–9.
World Health Organization. Dengue: guidelines for diagnosis, treatment, prevention and control. http://whqlibdoc.who.int/publications/2009/9789241547871_eng.pdf. Accessed 17 June 2022.
George L, et al. Communityeffectiveness of temephos for dengue vector control: a systematic literature review. PLoS Negl Trop Dis. 2015;9:e0004006.
GarzaRobledo AA, et al. Effectiveness of Spinosad and Temephos for the control of mosquito larvae at a tire dump In Allende, Nuevo Leon, Mexico. J Am Mosq Control Assoc. 2011;27:404–7.
Grisales N, et al. Temephos resistance in Aedes aegypti in Colombia compromises dengue vector control. PLoS Negl Trop Dis. 2013;7:e2438.
Ocampo CB, et al. Insecticide resistance status of Aedes aegypti in 10 localities in Colombia. Acta Trop. 2011;118:37–44.
Andersson N, et al. Evidence based community mobilization for dengue prevention in Nicaragua and Mexico (Camino Verde, the Green Way): cluster randomized controlled trial. BMJ. 2015;351:h3267.
Heintze C, Velazco Garrido M, Kroeger A. What do communitybased dengue control programmes achieve? A systematic review of published evaluations. Trans R Soc Trop Med Hyg. 2007;101:317–25.
Erlanger TE, Keiser J, Utzinger J. Effect of dengue vector control interventions on entomological parameters in developing countries: a systematic review and metaanalysis. Med Vet Entomol. 2008;22:203–21.
AlvaradoCastro V, et al. Assessing the effects of interventions for Aedes aegypti control: systematic review and metaanalysis of cluster randomised controlled trials. BMC Public Health. 2017;17(Suppl 1):384.
Andersson N, et al. The women made it work: fuzzy transitive closure of the results chain in a dengue prevention trial in Mexico. BMC Public Health. 2017;17(Suppl 1):408.
Esteva L, Vargas C. Analysis of a dengue fever disease transmission model. Math Biosci. 1998;150:131–51.
Esteva L, Vargas C. A model for dengue disease with variable human population. J Math Biol. 1999;38:220–40.
Garba SM, Gumel AB, Bakar MR. Backward bifurcations in dengue transmission dynamics. Math Biosci. 2008;215:11–25.
LiMartín A, ReyesCarreto R, VargasDeLeón C. Dynamics of a dengue disease transmission model with twostage structure in the human population. Math Biosci Eng. 2023;20(1):955–74.
Predescu M, et al. On the dynamics of a deterministic and stochastic model for mosquito control. Appl Math Lett. 2007;20:919–25.
Burattini MN, et al. Modelling the control strategies against dengue in Singapore. Epidemiol Infect. 2008;136:309–19.
Abidemi A, AbdAziz M, Ahmad R. Vaccination and vector control effect on dengue virus transmission dynamics: modelling and simulation. Chaos Solit Fractals. 2020;133:109648.
Abboubakar H, Kamgang JC, Tieudjo D. Backward bifurcation and control in transmission dynamics of arboviral diseases. Math Biosci. 2016;278:100–29.
Buonomo B, Della MR. Optimal bed net use for a dengue disease model with mosquito seasonal pattern. Math Methods Appl Sci. 2018;41:573–92.
Asamoah JKK, et al. Optimal control and costeffectiveness analysis for dengue fever model with asymptomatic and partial immune individuals. Results Phys. 2021;31:104919.
Naidu DS. Optimal control systems. 1st ed. Florida: CRC Press; 2002. p. 432.
Ndekha A, et al. Community participation as an interactive learning process: experiences from a schistosomiasis control project in Zimbabwe. Acta Trop. 2003;85:325–38.
Esteva L, Yang HM. Assessing the effects of temperature and dengue virus load on dengue transmission. J Biol Syst. 2015;23:527–54.
Ai S, Li J, Lu J. Mosquitostagestructured malaria models and their global dynamics. SIAM J Appl Math. 2012;72:1213–37.
Moulay D, AzizAlaoui MA, Cadivel M. The chikungunya disease: modeling, vector and transmission global dynamics. Math Biosci. 2011;229:50–63.
Consejo Nacional de Población, México. Documento metodológico: proyecciones de la Población 2010–2050. 2012. http://www.conapo.gob.mx/work/models/CONAPO/Resource/1529/2/images/DocumentoMetodologicoProyecciones2010_2050.pdf. Accessed 2 July 2022.
Moulay D, AzizzAlaoui MA, Kwon H. Optimal control of chikungunya disease: larvae reduction, treatment and prevention. Math Biosci Eng. 2012;9:369–92.
Chadee DD. Impact of preseasonal focal treatment on population densities of the mosquito Aedes aegypti in Trinidad, West Indies: a preliminary study. Acta Trop. 2009;109:236–40.
MoralesPérez A, et al. Aedes aegypti breeding ecology in Guerrero: crosssectional study of mosquito breeding sites from the baseline for the Camino Verde trial in Mexico. BMC Public Health. 2017;17:61–70.
SztankayGulyás M. Mosquito control with integrated method. Wiad Parazytol. 1972;18:629–33.
Saide PM, et al. Technical document for the implementation of interventions based on generic operational scenarios for Aedes aegypti control. Washington, D.C.: PAHO; 2019. p. 49.
Gunn, et al. Current strategies and successes in engaging women in vector control: a systematic review. BMJ Glob Health. 2018;3(1):e000366.
MungallBaldwin C. Women’s participation in the prevention and control of dengue using environmental methods in the global south: a qualitative metasynthesis. Int J Equity Health. 2022;21(1):1–23.
Arosteguí J, et al. Beyond efficacy in water containers: temephos and household entomological indices in six studies between 2005 and 2013 in Managua, Nicaragua. BMC Public Health. 2017;17(1):85–92.
Chowell G, et al. Spatial and temporal dynamics of dengue fever in Peru: 1994–2006. Epidemiol Infect. 2008;136(12):1667–77.
Chowell G, et al. Estimation of the reproduction number of dengue fever from spatial epidemic data. Math Biosci. 2007;208(2):571–89.
National institute of statistics and geography, INEGI. Population and housing census and counting. http://en.www.inegi.org.mx/programas/ccpv/1900/. Accessed 24 June 2022.
Soetaert K, Petzoldt T, Setzer RW. Solving differential equations in R: package deSolve. J Stat Softw. 2010;33:1–25.
Abidemi A, et al. Lyapunov stability analysis and optimization measures for a dengue disease transmission model. Physica A. 2022;602:127646.
Mincham G, et al. Development of a mechanistic dengue simulation model for Guangzhou. Epidemiol Infect. 2019;147:e125.
Shepard DS, et al. Economic impact of dengue illness in the Americas. Am J Trop Med Hyg. 2011;84:200–7.
Acknowledgements
We thank all the households who responded to the questionnaire and agreed to the entomological survey in their premises under the Camino Verde project, and the members of the field teams who conducted the survey and carried out the entomological inspections.
Funding
This research was supported by the UBS Optimus Foundation. FomixCONACYTGUERRERO also supported the research. (Grant Number 200802108541).
Author information
Authors and Affiliations
Contributions
Conceptualization: VMAC CVDL SPS NA. Formal analysis: VMAC CVDL ALM. Funding acquisition: NA ENA. Investigation: VMAC SPS AMP. Methodology: VMAC CVDL. Project administration: NA. Resources: ENA AMP. Software: VMAC CVDL ALM. Supervision: SPS NA AC. Validation: VMAC SPS JLS. Visualization: VMAC SPS BMSG JLS. Writing – original draft: VMAC SPS BMSG JLS. Writing – review & editing: VMAC SPS NA AC.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
The study protocol of Camino Verde was approved by the Research Ethics Board of CIETcanada (which had received the resources for the study) on 16 November 2009 and the Ethics Committee of the Centro de Investigación de Enfermedades Tropicales (CIET) at the Autonomous University of Guerrero on 27 November 2009. Local authorities and community leaders gave their consent for the study. We also asked every household for verbal consent to apply the survey questionnaire, collect saliva samples and perform entomological inspections.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
AlvaradoCastro, V.M., VargasDeLeón, C., ParedesSolis, S. et al. The influence of gender and temephos exposure on community participation in dengue prevention: a compartmental mathematical model. BMC Infect Dis 24, 463 (2024). https://doi.org/10.1186/s1287902409341w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1287902409341w