Full text

        
10.1186/s40409-017-0123-x

Research article - Vol. 23, 2017

 

Mathematical modelling of vector-borne diseases and insecticide resistance evolution

 

Maria Laura Gabriel Kuniyoshi1, Fernando Luiz Pio dos Santos1 

 

1 Department of Biostatistics, Institute of Biosciences of Botucatu, São Paulo State University (UNESP - Universidade Estadual Paulista), Street Prof. Dr. Irina Delanova Gemtchujnicov, no number, Rubião Júnior, zip code 18618-693, PO box 510, Botucatu, SP, Brazil

 

ABSTRACT

Background:

Vector-borne diseases are important public health issues and, consequently, in silico models that simulate them can be useful. The susceptible-infected-recovered (SIR) model simulates the population dynamics of an epidemic and can be easily adapted to vector-borne diseases, whereas the Hardy-Weinberg model simulates allele frequencies and can be used to study insecticide resistance evolution. The aim of the present study is to develop a coupled system that unifies both models, therefore enabling the analysis of the effects of vector population genetics on the population dynamics of an epidemic.

Methods:

Our model consists of an ordinary differential equation system. We considered the populations of susceptible, infected and recovered humans, as well as susceptible and infected vectors. Concerning these vectors, we considered a pair of alleles, with complete dominance interaction that determined the rate of mortality induced by insecticides. Thus, we were able to separate the vectors according to the genotype. We performed three numerical simulations of the model. In simulation one, both alleles conferred the same mortality rate values, therefore there was no resistant strain. In simulations two and three, the recessive and dominant alleles, respectively, conferred a lower mortality.

Results:

Our numerical results show that the genetic composition of the vector population affects the dynamics of human diseases. We found that the absolute number of vectors and the proportion of infected vectors are smaller when there is no resistant strain, whilst the ratio of infected people is larger in the presence of insecticide-resistant vectors. The dynamics observed for infected humans in all simulations has a very similar shape to real epidemiological data.

Conclusion:

The population genetics of vectors can affect epidemiological dynamics, and the presence of insecticide-resistant strains can increase the number of infected people. Based on the present results, the model is a basis for development of other models and for investigating population dynamics.

 

Keywords: Epidemiology; Population genetics; Tropical diseases; Insecticides; Theoretical modelling; Numerical simulation; ODE system

 

Received: December 21, 2016.

Accepted: June 13, 2017.

 

Correspondence: marialgk@aluno.ibb.unesp.br

 

Authors’ contributions

MLGK participated in the development of the mathematical model, execution and design of the numerical simulations, analysis of the results and preparation of the manuscript. FLPS participated in the design of the numerical simulations, analysis of the results and preparation of the manuscript. All authors read and approved the final manuscript.

 

Competing interests

The authors declare that they have no competing interests.