Abstract

In this paper, a novel age-structured delayed mathematical model to control Aedes aegypti mosquitoes via Wolbachia-infected mosquitoes is introduced. To eliminate the deadly mosquito-borne diseases such as dengue, chikungunya, yellow fever, and Zika virus, the Wolbachia infection is introduced into the wild mosquito population at every stage. This method is one of the promising biological control strategies. To predict the optimal amount of Wolbachia release, the time varying delay is considered. Firstly, the positiveness of the solution and existence of both Wolbachia present and Wolbachia free equilibrium were discussed. Through linearization, construction of suitable Lyapunov–Krasovskii functional, and linear matrix inequality theory (LMI), the exponential stability is also analyzed. Finally, the simulation results are presented for the real-world data collected from the existing literature to show the effectiveness of the proposed model.

1. Introduction

Mosquito-borne diseases represent the vertical transmission of bacteria and viruses from mosquitoes to human while female mosquito taking a blood meal. Mosquito-borne diseases such as dengue, chikungunya, yellow fever, Zika virus, and Japanese encephalities cause over one million deaths per annum [1, 2]. Gubler in [3, 4] explained that the dengue and dengue hemorrhagic fever are the most common issues for public health. The primary vector for most of the mosquito-borne diseases is Aedes aegypti, and recently Aedes albopictus also add as a secondary vector [5]. In the past sixty years, the spread of mosquito-borne diseases has increased dramatically [6]. More than that, per year, dengue causes nearly 20 thousand deaths all over the world [7]. Also, nearly 112 countries are attacked by mosquito-borne diseases [8].

In recent years, there are several articles are available to control vectors by genetic modifications [9]. Moreover, some biological control methods to replace the wild mosquitoes by releasing genetically modified mosquitoes are also tried by some researchers. Those biological control methods are sterilization of male mosquitoes [10, 11] and genetic modification to reduce the reproduction and increase the life-shortening bacteria Wolbachia [12]. Furthermore, via finding the reproduction number of a mathematical model which depicts the virus transmission via human sexual contact was analyzed in [13]. In [14], the author tried some other types of control agents such as bed nets, mosquito repellents, indoor residual spray, condoms during sex, medically treating infected human, and quarantine. However, in [15], the author tried to control mosquitoes by making modifications in feeding behaviours.

In this environment, our main aim is to control the vector population that transmits the virus to the uninfected human while taking a blood meal. There is a life-shortening bacterium called Wolbachia which will be very useful to reach our aim; in [16], Mcmeniman et al. analyzed the stable introduction of Wolbachia in Aedes mosquitoes. Wolbachia is a Gram-negative bacterium, and it is first reported in the tissues of the mosquito Culex pipients (Hertig and Wolbach, 1924) [17]. In recent results, they found that yellow fever virus can also be blocked by Wolbachia [18].

If a mosquito carries this bacterium, then the virus inside the mosquito does not get transmitted into the uninfected human. It blocks the virus inside the mosquito at salivary gland. This can be understood with Figure 1.

In Figure 1, the process of releasing Wolbachia bacteria into mosquito population is as follows:(1)In laboratory, the Wolbachia pipients are injected into eggs, larvae, and pupae of Aedes aegypti via microinjection.(2)Cytoplasmic incapability (CI): the adult Wolbachia-infected mosquitoes which are reared at the laboratory are released to the wild mosquito population of Aedes aegypti. Through this process, there exist three types of possibilities which are(i)If the Wolbachia-infected female mosquitoes mate with the Wolbachia-infected male, then the progeny should have the Wolbachia by birth which is compatible.(ii)If the Wolbachia-infected female crosses with Wolbachia-uninfected male, then the progeny face the same problems as in (i).(iii)If Wolbachia-uninfected female crosses with Wolbachia-infected male, then there is no viable progeny.

These two processes can be virtually understood by Figure 2.

In eggs, larvae, and pupae population, we can microinject the Wolbachia and release this in patches at dengue-suspected areas. This process is practically done by placing “Zancu kits” around the people living areas. And the adult mosquitoes which are reared at lab can also be released into wild mosquito population. This process is called “introgression.”

Various mathematical models have been studied to understand the interplay among Wolbachia and non-Wolbachia mosquitoes. In [19], the author considered the Wolbachia bacteria as a mechanism to control arbovirus, and his experimental studies show that the spread of Zika virus among mosquitoes and human was notably reduced. The same process for dengue virus spread was studied by Segoli et al. in [20]. In [21], the author has created a mathematical model considering only adult female mosquitoes and converted the model into endoepidemic model consisting of adult female mosquitoes and human population. In [22], the authors proposed a mathematical model depicting the life stages of mosquitoes with Wolbachia and proved that Wolbachia has excellent quality to control dengue virus spread. In that work, Koiler et al. discussed the virus as well as Wolbachia in both mosquitoes and human. Also, the objective is to predict the appropriate release of this bacterium, and the basic reproduction number was analyzed. Supriatna et al. in [23] developed a mathematical model to express the dynamics of dengue virus in both human and Aedes aegypti mosquitoes. In that, the human vaccination and Wolbachia introduction were used as optimal control methods. In [24], authors discussed the birth and death rate impulsive model to control mosquito-borne diseases using Wolbachia via the stroboscopic map method. The integer-order mathematical model which describes the interplay among the wild and Wolbachia-infected mosquitoes was analyzed in [25]. In that work, the author divided the mosquito population into two groups: one is aquatic and another one is adult. In [26], the author proposed a mathematical model to describe the persistence of Wolbachia via two-sex stage-structured model.

Hence, with a full understanding of the interplay among the Wolbachia and non-Wolbachia mosquitoes in our work, we have created a mathematical model consisting of 10 stages to ensure the success of the proposed strategy. With reference of the practical results in [18], we can release the Wolbachia in every stages in the forms of “Zancu kits” and “Introgression.” So to obtain an optimal control, it is necessary to consider each and every variable. Because we know that in mathematical modeling, each and every parameter plays an important role. Up to our knowledge, this is the first article considering the control inputs in 10 stages incorporated with time-varying delays.

By motivated by the above discussions, the main contribution of this paper is as follows:(i)The main aim is to establish a novel mathematical model to describe the interplay among the both non-Wolbachia (wild mosquitoes) and Wolbachia-infected mosquitoes with time-varying delay.(ii)We found that our method will increase Wolbachia-infected mosquitoes in terms of CI rescue team and non-Wolbachia mosquitoes go to annihilation.(iii)We have proved that the releasing of adult female Wolbachia-infected mosquitoes is more beneficial than the releasing of adult male Wolbachia-infected mosquitoes.(iv)There exists no literature on exponential stability results of delayed Wolbachia and non-Wolbachia age-structured model. This gap is filled by our work. Finally, by using real-world data, we checked the dynamics of the proposed model by using MATLAB LMI tool box.

The rest of the paper is arranged as follows: in Section 2, the novel mathematical model which describes the interplay among the Wolbachia free and Wolbachia present mosquito population is proposed. In Section 3, the analysis of the model such as positiveness of the solution and existence of equilibrium points are discussed. In Section 4, the exponential stability results of the linearized delayed system with time-varying delay is presented. In Section 5, numerical simulation results are presented. The work is concluded in Section 6.

2. Preliminaries and Model Formulation

In this section, some basic definitions and lemmas which are used to derive our results are presented. And the interaction between wild mosquitoes and Wolbachia-infected mosquitoes is modeled.

Definition 1 (see [27]). A model is said to be exponentially stable at its equilibrium point, if there exists such that

Lemma 1. (Schur Complement, see [28]). Let us denote three matrices as , where and . Then if and only if or .

Lemma 2. (see [29]). For any scalar , , and matrix , then

2.1. Modeling the Life Stages of Wild Mosquito Population

In a common environment, Aedes aegypti mosquito population has five important life stages such as eggs , larvae , pupae , matured female mosquitoes , and adult male mosquitoes (Figure 3). These life stages with respect to time can be modeled as follows:

2.2. Modeling the Interaction between Wild and Wolbachia-Infected Mosquitoes

In this subsection, we modeled the release of Wolbachia-infected mosquitoes into the wild mosquitoes in the mosquito-borne disease suspected areas. The Wolbachia infection is released in both aquatic (Zanku Kits) and ariel stages. Therefore, by using the language of mathematics, we can model the interaction between wild and Wolbachia mosquitoes as follows:

To understand the system of equations in (4), refer Figure 3.

In the above mathematical model, the production of eggs can be found by the term . That is, eggs without Wolbachia infection are produced by the mating between wild female and wild male mosquitoes, where is the total population which can be calculated by the following expression:

Along with this, the terms (natural mortality rate of non-Wolbachia eggs) and (maturation rate of non-Wolbachia eggs) denote the limitations in the growth of wild mosquito eggs. At the same time, after release of Wolbachia-infected mosquitoes (in both aquatic and ariel stages) in a common environment, the production of Wolbachia-infected mosquito eggs depends on mating between Wolbachia-infected female and non-Wolbachia male and from mating between Wolbachia-infected female and Wolbachia-infected male . This implies that the birth rate of Wolbachia-infected mosquito egg population with the reproduction rate is

Similarly, the increase in growth is limited by the natural mortality rate and the maturation rate (that is, the rate in which the corresponding compartment moved into the next stage).

Furthermore, is added to the wild mosquito larvae population. Because the terms and denote the probability of getting larvae with and without Wolbachia, respectively. Similarly, and denote the probability of getting pupae with and without Wolbachia, respectively, and denote the probability rate of having Wolbachia infection in adult mosquitoes by introgression. That is, be the probability of getting Wolbachia-infected adults (with  = probability of getting male and  = probability of getting female). Because of these reasons, the terms , , , and are added to the corresponding stages, and similarly, the terms , and are removed from the corresponding stages. And the other parameters used in this model are described in Table 1.

3. Analysis of the Model

In this section, we analyze the positivity, existence of equilibrium points, and stability of the system of equations in (4).

3.1. Positivity of Solutions

Lemma 3. For all , the solutions are all nonnegative if the initial values , , , , , , , , , .

Proof. Let us defineFrom (7), . Let us consider the first equation of model (4), that is,The integrating factor isMultiply (9) with (8) on both sidesThat is,Integrating on both sides with respect to ,This implies that .
Similarly, we can prove that for the positive initial values, the solution is positive.

3.2. Existence of Equilibrium Points

In this section, the possible cases of existing equilibrium points are discussed. Furthermore, the null mosquitoes equilibrium point is omitted because this case does not exist in nature.

3.2.1. Wolbachia Free Equilibrium

In this subsection, we can find Wolbachia free equilibrium by equating the system of equations in (2) to zero and putting , , , , and .

That is,

By solving equations (13), we get the equilibrium point aswhere

3.2.2. Wolbachia Present Equilibrium

In this subsection, we can find Wolbachia present equilibrium by equating the system of equations in (4) to zero. That is,

By solving the system of equations in (16), we get the following equilibrium point in terms of :with , both roots can be found from the quadratic equation:wherewhere

3.2.3. Existence of Wolbachia Free and Wolbachia Present Equilibrium

The equilibrium points exist if it satisfies the following conditions:(i) exists provided that and .(ii) and exist provided that and .

4. Stability Analysis

In this section, the stability results of the system of equations in (2) by linearization, Lyapunov–Karasovskii functional, and LMI approach. To optimize our differential model, we consider the delay terms into account as discussed in [30] i.e., where is the produced wild egg density at time from which the wild larvae density at time is produced, is the produced wild larvae density at time from which the wild pupae density at time is produced, is the produced wild pupae density at time from which the wild female or male density at time is produced, is the produced Wolbachia-infected egg density at time from which the Wolbachia-infected larvae density at time is produced, is the produced Wolbachia-infected larvae density at time from which the Wolbachia-infected pupae density at time is produced, and is the produced Wolbachia-infected pupae density at time from which the Wolbachia-infected female or male density at time is produced. These considerations make our system of differential equations to system of delay differential equations:

By linearizing, the following system describes the interaction between wild mosquitoes and Wolbachia mosquitoes aswhere

Assume that , and it is bounded, differentiable, and nonnegative. Now, consider an equilibrium point of system (22) asLet , and this implies that . And . Then the modified system iswith be the initial condition.

In general, stability is an important requirement to study dynamical systems [3133]. In following theorem, we derived the sufficient for exponential stability results for system (25).

Theorem 1. The system (25) is said to be exponentially stable if there exist positive definite matrices and and positive scalars and such thatwith . Moreover,

Proof. Let us consider the following Lyapunov–Krasovskii function:Now, the time derivative of the Lyapunov–Krasovskii functional along with the trajectories of system of equation (25) Put , the linear feedback control,By using Lemma 2, we can get the following inequalities:, wherewhere .Hence, by Lemma 1, . Moreover, . Here,where is the maximum eigen value of .Here,Now,This implies thatand by Definition 1, our system (25) is exponentially stable.

Remark 1. If the system (25) is without time-varying delay terms, then the results of Theorem 1 have changed as follows: the possible Lyapunov function is with the sufficient conditionsThen the system (25) is without time-varying delay is exponentially stable.

5. Numerical Simulations

The numerical simulation results of system of equations in (21), based on the real-world data, are mentioned in Table 2.

Consider the following system:with numerical valueswith the control gain matrixwith and .

Remark 2. The numerical simulations derived in our paper shows that the allowable upper bound of the time-varying delay . However, the existing literature used the classical method such as Pontryagin maximum principle and Jacobian matrix method, by finding the sign of eigen values, fractional-order mathematical models [38, 39], and through finding the reproduction number to analyze the stability of the model. In our work, we have used LMI approach. However, compared with the other algebraic methods LMI approach is less conservative.
Figures 4 and 5 depict the population dynamics of wild mosquitoes and Wolbachia-infected mosquitoes at the initial release of Wolbachia, respectively. From Figure 6, we can observe that the Wolbachia-infected population also tends to decrease in numbers. Our aim is to maintain the Wolbachia population in certain quantity to spread Wolbachia among wild mosquitoes naturally. Because of this reason, we have released another batch of Wolbachia-infected mosquitoes as an optimal control. After this release, we can observe from Figures 7 and 8 that the level of Wolbachia-infected mosquitoes is increased to certain quantity and wild mosquito population is decreased to zero.

Remark 3. For the optimal control to ensure the successful release of Wolbachia, the data from Table 2 are substituted in the system of equations in (21). From this we obtained the following results: Figure 4 denotes the dynamics of non-Wolbachia mosquitoes, and Figure 5 denotes the dynamics of Wolbachia-infected mosquitoes. We can observe that the Wolbachia release in two main stages is not enough to control wild mosquitoes. In this, we noticed that the Wolbachia-infected mosquito population also goes to annihilation. Refer Figure 6, when comparing wild mosquitoes with Wolbachia-infected mosquitoes, it is increased and Wolbachia is not carried out to next generation. By implementing our proposed strategy, if we released Wolbachia in every stage, the Wolbachia-infected population is increased to certain level and maintained at that level (Figure 7) for the remaining period; at the same time, the Wolbachia population is annihilated without any increase (Figure 8).

6. Conclusion

A novel integer-ordered age-structured delayed mathematical model to describe the interplay among the wild and Wolbachia-infected mosquitoes was proposed. The positiveness of the solution and existence of both Wolbachia present and Wolbachia free equilibriums were analyzed. By using some mathematical techniques such as linearization, Lyapunov–Krasovskii functional, and LMI theory, the exponential stability results of the proposed model were established. By using MATLAB, the numerical simulation results were presented to show the effectiveness of the created mathematical model. Our results show that the release of adult female Wolbachia-infected mosquitoes is more beneficial than the release of adult male Wolbachia-infected mosquitoes.

Data Availability

Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This article has been written with the joint partial financial support of SERB-EEQ/2019/000365, the National Science Centre in Poland (Grant DEC-2017/25/B/ST7/02888, RUSA Phase 2.0 (Grant No. F 24–51/2014-U), Policy (TN Multi-Gen), Dept.of Edn. Govt. of India, UGC-SAP (DRS-I) (Grant No. F.510/8/DRS-I/2016(SAP-I)), DST-PURSE 2nd Phase programme vide letter No. SR/PURSE Phase 2/38 (G), DST (FIST—level I) 657876570 (Grant No. SR/FIST/MS-I/2018/17), and Prince Sultan University for funding this work through research group Nonlinear Analysis Methods in Applied Mathematics (NAMAM) group number RG-DES-2017-01-17.