Application of Mathematical Methods in Nature-Inspired ComputationView this Special Issue
Combined Trojan Y Chromosome Strategy and Sterile Insect Technique to Eliminate Mosquitoes: Modelling and Analysis
Sterile insect technique has been successfully applied in the control of agricultural pests; however, it has a limited ability to control mosquitoes. A promising alternative approach is the Trojan Y Chromosome strategy, which works by manipulating the sex ratio of a population through the introduction of YY supermales that guarantee male offspring. To take the advantages of both approaches, a combined Trojan Y chromosome strategy and sterile insect technique (TYC-SIT) strategy considering intraspecies competition is modeled. The pure TYC method is compared with the pure SIT method by cancelling one-state variable. The dynamical analysis leads to results on both local and global stabilities of this combined TYC-SIT model. Optimal control analysis is also implemented to investigate the optimal mechanisms for extinction of mosquitoes. In particular, the numerical results affirm that the combined TYC-SIT enables near elimination of mosquitoes and works better than the pure TYC or pure SIT method. These conclusions have great significance for species controls with an XX-XY sex determinism system or ZZ-ZW system.
Mosquito-borne diseases are transmitted by mosquitoes infected with viruses, such as Zika virus, yellow fever virus, West Nile fever virus, and dengue fever virus [1, 2]. The spread of mosquito-borne virus in humans is mainly through the bite of mosquitoes infected with the virus such as Anopheles sinensis, Anopheles anthropophagus, Aedes aegypti, Aedes albopictus, and Culex mosquitoes [3–5]. Frequent outbreaks of mosquito-borne diseases in tropical and subtropical regions have caused large numbers of deaths and economic losses. Chemical treatments such as pesticides have been implemented for many years to eradicate mosquitoes. Heavy spraying of insecticides has effectively suppressed mosquito density. However, environmental problems caused by excessive use of pesticides , insecticide resistance [7, 8], and combined with the lack of vaccines [9, 10], have called for alternative environment-friendly and sustainable approaches , such as radiation-based sterile insect technique (SIT) and Trojan Y chromosome strategy (TYC).
SIT works by releasing radiation-sterilized males to an existing population, to mate with wild females so that they have no viable offspring [11, 12]. SIT has been successfully applied in the control of several agricultural pests such as invasive fruit flies, lepidopteran, and Hemiptera [13–17]; however, it has a limited ability to control mosquitoes due to the three issues. (1) It requires introducing large quantities of sterile mosquitoes. (2) The health of sterile male mosquitoes is significantly impaired, thus leading to the decline in mating competitiveness. (3) Mosquitoes are very prolific; the population may be supported even with a low mating ratio of normal females to normal males [18–21].
A promising approach that has been put forward for eliminating mosquitoes is TYC, in which supermales (containing two Y chromosomes) are released into the field to mate with wild females, resulting in a sharp sex imbalance in the offspring [22–25]. Figure 1 illustrates that an equal proportion of females and males are produced in the wild; however, the offspring is guaranteed to be male when a natural female mates with a YY supermale. The gradual reduction in females and increase in males may lead to eventual extinction of the targeted population after several generations . The advantages of TYC strategy are as follows. (1) It only requires to recombine sex chromosomes of mosquitoes instead of changing their genes. (2) Only a fixed number of modified mosquitoes are required. (3) The intensity of this strategy can be controlled.
However, it is expensive to design and produce YY supermales, that is why TYC is not widely used in controlling mosquitoes or agricultural pests [27, 28]. In this study, a combined TYC and SIT model is proposed to take advantages of both approaches. Intraspecies competition is especially considered, which is omitted in many literature. By setting different values to the release rate of radiation-based sterile males and the release rate of supermales with YY chromosome structure, we can switch to a pure SIT model or a pure TYC model. To the best of our knowledge, it is the first mathematical investigation to compare the combined SIT-TYC to TYC or SIT.
2. Materials and Methods
2.1. Mathematical Modelling
The population of mosquitoes in the wild tends to be numerous because their reproduction rates are high and matings occur constantly; therefore, continuous models of ordinary differential equations (ODEs) can be established to describe the population dynamics of mosquitoes . Parameters will be used in this work and are firstly explained in Table 1.
The combined TYC-SIT model, in which both YY supermales and radiation-based sterile males are introduced, is described by a system of four ODEs for state variables: a wild-type XX female , a wild-type XY male , a sterile male , and a YY supermale . Half female and half male offspring are produced in the wild; however, only males can be produced if a wild female mates with a YY supermale, and no viable offspring can be produced if mates with a sterile male. Thus, the set of equations that describes the system iswhere define the number of individuals in each associated class, and
The intraspecies competition for female mates caused by the introduction of YY supermales and radiation-based sterile males is considered and modeled as
If , that is, , it implies that there is no mating pressure on wild males. Obviously,
The larger the value of is, the more competitive wild males are. Similarly, the larger the value of is, the more competitive YY supermales are.
2.2. Equilibrium and Stability Analysis
The stability of the TYC-SIT model (1) is now investigated. There is one equilibrium, , on the boundary. It is clear is global stable. To get the positive interior equilibrium explicitly, i.e., , it is equivalent to solve the following equations:
By solving equations (5), , and can be obtained as follows:
And is the solution of
Theorem 1. Let . The interior equilibrium is locally asymptotically stable if
Proof. The Jacobian matrix of model (1) about is given aswhereThe corresponding characteristic equation iswhereAccording to Routh–Hurwitz criterion, is locally asymptotically stable if
Theorem 2. In the case of , the trivial equilibrium of TYC-SIT model (1) is globally asymptotically stable if .
Proof. Consider the Lyapunov function . Note that due to the positivity of its solutions and is radially unbounded. It is left to show that for all . Taking the derivative of about yieldsTherefore, we only need to show to get global stability, and this requires , which completes the proof.
2.3. Optimal Control Analysis
The goal in this section is to investigate the mechanisms in TYC-SIT system to lead to an optimal level of female mosquitoes density. We assume that the influx of sterile males and YY supermales and are not known a priori and enter them into the system as time-dependent controls and . In fact, male mosquitoes do not bite human and only female mosquitoes bite and spread diseases [30–32], which implies that we do not have to kill all mosquitoes, eliminating females is enough instead. Furthermore, we also hope the production of sterile males and YY supermales is minimized. Herein, the following objective function is chosen:where subject to the governing equation (1).
Optimal strategies are derived from the objective function, where the female population is minimized and also the introduction of both radiation-based sterile males and supermales with YY chromosome structure are minimized. We search for the optimal controls within the set , which is given by
The goal is to find the optimal controls such that
Theorem 3. There exists the optimal controls and for the minization problem (17).
Proof. It is obvious that is concave in the argument and ; all state variables are bounded by the carring capacity , and the control variables and are also assumed bounded due to the cost of releasing sterile males and YY supermales; therefore, the existence of for minimization problem (16) is guaranteed based on [33–35].
Theorem 4. An optimal control of system (1) that minimizes the objective function is characterized by
Proof. Here, Pontryagin’s minimum principle is used to derive the necessary conditions on this problem. The Hamiltonian in this problem isThe Hamiltonian is used to find the adjoint functions :To find the optimal , minimize pointwise, that is,Note that cancels with 2 which comes from the square of the controls and . Furthermore, the problem is indeed minimization asHence, the optimal solutions areA compact way of writing the optimal control isNow, the proof is completed.
3. Numerical Simulations
The numerical simulations are investigated by MATLAB R2019b with the values of initial conditions and parameters shown in Table 2. The ode15s solver was used to get numerical solutions of the combined TYC-SIT system. The TOMLAB Base Module and TOMLAB/SNOPT are also used to solve the optimal control problems of our dynamic systems.
The combined TYC-SIT system is modeled, and the initial conditions and parameters in Table 2 are utilized to observe the relative population decline of females in response to the addition of the radiation-based sterile males and YY supermales.
Under the condition of , that is, sterile males or YY supermales are added to the population only once at time and no additional will be introduced, it is observed that no matter how large the initial introduction of sterile males or YY supermales are, the system of the combined TYC-SIT cannot achieve extinction and instead leads to an equilibrium state. As shown in Figure 2, the population did decline for some time with large enough (i.e., purple/green star line) influx of sterile males or YY supermales; however, the population recovers soon and reaches the equilibrium state at approximately 172, which can be also calculated from equations (5) and (6). Extinction can occur with continuous introducing modified males; two examples are provided in Figure 3.
To further validate the effectiveness of the proposed TYC-SIT model, comparative analysis among the pure TYC, the pure SIT, and the combined TYC-SIT is carried out.
If we solely continue adding sterile males (i.e., ) or solely continue adding YY supermales (i.e., ) after the initial introduction of the both, how many modified males would be introduced has a great impact on the population decline of mosquitoes. With a relatively low influx of modified males, purely introducing sterile males or YY supermales cannot drive the population to extinction; however, the effectiveness of the continuous introduction of YY supermales is better because the population drops much more, see part (a) in Figure 4. As or is increased to 10, the initial decline rate with is faster than , but then the situation is reversed after the threshold point in Figure 4 part (b). As we can see in part (c) and (d) in Figure 4, the population decline rate with either or is similar if the influx of modified males is large enough.
The left panel of Figure 5 shows the female density change varies with time under different influx of sterile males, and the right panel shows the female density change varies with time under different influx of YY supermales. Identical parameters are used to compare the effectiveness of SIT and TYC. TYC is much more effective with low influx of modified males. Furthermore, combined TYC and SIT requires less input of modified males, as shown in Figure 6.
Considering the cost of the production of sterile males and YY supermales, optimal strategies and its corresponding optimal states have also been numerically simulated to compare the effects of the combined TYC-SIT strategy with the pure SIT and the pure TYC strategy, see Figures 7–9. The population can be driven to extinction with an optimal time controls or by using either SIT , or TYC , or the combined TYC-SIT approach.
Figure 9 shows the optimal eradication strategy of the combined TYC-STT model requires to introduce YY supermales at 20% of the initial female population (or 5% of the carrying capacity) for a short time to bring the population of females below some threshold, then gradually drops to a relatively low level after introducing sterile males at time , and finally turns off until the entire population vanishes. However, both pure TYC and pure SIT require to release at 30% of the initial female population, as shown in Figures 7 and 8. In overall, the combined TYC-SIT approach works better in eradicating mosquitoes, which reduces the whole releasing rate by 33%.
4. Conclusions and Discussion
In this study, we established a mathematical model of the combined TYC and SIT approaches. The population was divided by the following four state variables: wild female mosquitoes , wild male mosquitoes , radiation-based sterile males , and YY supermales . Six parameters, the birth coefficient , the death coefficient , the logistic term L, the carrying capacity K, the influx of sterile males , and the influx of sterile males , are included. The intraspecies competition for female mates caused by the introduction of modified male mosquitoes is especially considered, which is omitted in many studies. By setting or , pure SIT is compared with pure TYC strategy. If only a single method is considered, the TYC strategy is better than SIT with low influx of modified males, and the situation is reversed with high influx of modified males. The dynamical analysis and optimal control analysis of the combined TYC-SIT model in this study are important to understand the efficiency of the combined strategy to eliminate mosquitoes. Optimal control analysis also provides the optimal release rates of YY supermales and sterile males. In particular, the comparative studies confirm that the combined TYC-SIT approach can indeed eliminate mosquitoes and perform better than pure TYC or pure SIT.
The combined TYC-SIT approach is a safe biological control method because it has the advantage of being species-specific and does not harm other nontarget species. In addition, older generations would not transfer genetically engineered genes to their offspring. The intensity of this method is also controllable because the number of radiation-based sterile males or supermales with YY chromosome structure to be released to the targeted population can be determined by the executor. Unlike other strategies, the combined TYC-SIT does not rely on eliminating all wild matings to affect the total population. Instead, it depends on the number of wild females dwindling over a period of several generation cycles. The results of this research are of great significance for the biological control of mosquitoes.
Besides mosquitoes, the combined TYC-SIT strategy can also be applied to eradicate species with an XX-XY sex determinism system or ZZ-ZW system. Our results will help experts devise suitable management strategies and measures to control invasive species or agricultural pests.
|SIT:||Sterile insect technique|
|TYC:||Trojan Y chromosome strategy|
|ODEs:||Ordinary differential equations.|
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare no conflicts of interest.
J.L., C.K., W.S., and M.G conceptualized the study; J.L. carried out formal analysis; J.L. and C.K. investigated the study; J.L. prepared and wrote the original draft; C.K., M.G., and W.S. reviewed and edited; J.L., C.K., and W.S. administrated the project; J.L. and M.G. carried out funding acquisition All authors have read and agreed to the published version of the manuscript.
The authors thank Dr. Chuan He for his support throughout this work. This research was funded by the Key Laboratory of Pattern Recognition and Intelligent Information Processing, Institutions of Higher Education of Sichuan Province (Grant no. MSSB-2021-09), and Talent Initiation Program of Chengdu University (Grant no. 2081921002).
J. Lyu, M. Gu, and S. Wang, “Combined trojan Y chromosome strategy and sterile insect technique to eliminate mosquitoes: modelling and analysis,” 2021, https://arxiv.org/abs/2108.08285.View at: Google Scholar
X. Yu, Y. Zhu, X. Xiao, P. Wang, and G. Cheng, “Progress towards understanding the mosquito-borne virus life cycle,” Trends in Parasitology, vol. 35, no. 12, pp. 1009–1017, 2019.View at: Publisher Site | Google Scholar
R. Ross, “Researches on malaria,” Journal of the Royal Army Medical Corps, vol. 4, no. 6, pp. 705–740, 1905.View at: Google Scholar
L. Esteva and C. Vargas, “Analysis of a Dengue disease transmission model,” Mathematical Biosciences, vol. 150, no. 2, pp. 131–151, 1998.View at: Publisher Site | Google Scholar
M. Q. Benedict, R. S. Levine, W. A. Hawley, and L. P. Lounibos, “Spread of the tiger: global risk of invasion by the mosquito Aedes albopictus,” Vector Borne and Zoonotic Diseases, vol. 7, no. 1, pp. 76–85, 2007.View at: Publisher Site | Google Scholar
C. D. Milam, J. L. Farris, and J. D. Wilhide, “Evaluating mosquito control pesticides for effect on target and nontarget organisms,” Archives of Environmental Contamination and Toxicology, vol. 39, no. 3, pp. 324–328, 2000.View at: Publisher Site | Google Scholar
Y. Yu, “Studies on insecticde-resistance in mosquitoes II. Laboratory experiments on the selection of resistance to DDT and BHC in the adults of coquillett,” Acta Entomologica Sinica, vol. 13, no. 3, pp. 339–343, 1964.View at: Google Scholar
R. W. Fay, J. W. Kilpatrik, R. L. Crowell, and K. D. Quarterman, “A method for field detection of adult-mosquito resistance to ddt residues,” Bulletin of the World Health Organization, vol. 9, no. 3, pp. 345–351, 1953.View at: Google Scholar
R. S. Dhillon, D. Srikrishna, and A. K. Jha, “Containing zika while we wait for a vaccine,” BMJ, vol. 356, Article ID j379, 2017.View at: Publisher Site | Google Scholar
D. K. Rocha, O. Matosc, M. T. Novoa, and C. Moiteiro, “Larvicidal activity against Aedes aegypti of Foeniculum vulgare essential oils from Portugal and Cape Verde,” Natural product communications, vol. 10, 2015.View at: Google Scholar
V. A. Dyck, J. Hendrichs, and A. S. Robinson, Sterile Insect Technique: Principles and Practice in Area-wide Integrated Pest Management, Springer, Netherlands, pp. 3–36, 2005.
D. W. Dunn and P. A. Follett, “The sterile insect technique (SIT) - an introduction,” Entomologia Experimentalis et Applicata, vol. 164, no. 3, pp. 151–154, 2017.View at: Google Scholar
K. Bourtzis and J. Hendrichs, “Preface: development and evaluation of improved strains of insect pests for sterile insect technique (SIT) applications,” BMC Genetics, vol. 15, no. 2, 2014.View at: Publisher Site | Google Scholar
K. Bourtzis and M. J. B. Vreysen, “Sterile insect technique (SIT) and its applications,” Insects, vol. 12, no. 7, p. 638, 2021.View at: Publisher Site | Google Scholar
C. Pascacio-Villafán, L. Quintero-Fong, L. Guillén, J. P. Rivera-Ciprian, R. Aguilar, and M. Aluja, “Pupation substrate type and volume affect pupation, quality parameters and production costs of a reproductive colony of ceratitis capitata (Diptera: tephritidae) VIENNA 8 genetic sexing strain,” Insects, vol. 12, no. 4, p. 337, 2021.View at: Publisher Site | Google Scholar
A. B. Dhahbi, Y. Chargui, S. M. Boulaaras, and S. B. Khalifa, “A one-sided competition mathematical model for the sterile insect technique,” Complexity, vol. 2020, Article ID 6246808, 12 pages, 2020.View at: Publisher Site | Google Scholar
C. F. Marina, P. Liedo, J. G. Bond et al., “Comparison of ground release and drone-mediated aerial release of Aedes aegypti sterile males in southern Mexico: efficacy and challenges,” Insects, vol. 13, no. 4, p. 347, 2022.View at: Publisher Site | Google Scholar
D. A. Dame, C. F. Curtis, M. Q. Benedict, A. S. Robinson, and B. G. J. Knols, “Historical applications of induced sterilisation in field populations of mosquitoes,” Malaria Journal, vol. 8, no. 2, 2009.View at: Publisher Site | Google Scholar
R. S. Lees, J. R. Gilles, J. Hendrichs, M. J. Vreysen, and K. Bourtzis, “Back to the future: the sterile insect technique against mosquito disease vectors,” Current Opinion in Insect Science, vol. 10, pp. 156–162, 2015.View at: Publisher Site | Google Scholar
X. Zheng, D. Zhang, Y. Li et al., “Incompatible and sterile insect techniques combined eliminate mosquitoes,” Nature, vol. 572, pp. 56–61, 2019.View at: Publisher Site | Google Scholar
M. E. H. Helinski, A. G. Parker, and B. G. J. Knols, “Radiation biology of mosquitoes,” Malaria Journal, vol. 8, no. 2, p. S6, 2009.View at: Publisher Site | Google Scholar
J. B. Gutierrez and J. L. Teem, “A model describing the effect of sex-reversed YY fish in an established wild population: the use of a Trojan Y-Chromosome to cause extinction of an introduced exotic species,” Journal of Theoretical Biology, vol. 241, no. 2, pp. 333–341, 2006.View at: Publisher Site | Google Scholar
S. Cotton and C. Wedekind, “Control of introduced species using Trojan sex chromosomes,” Trends in Ecology & Evolution, vol. 22, no. 9, pp. 441–443, 2007.View at: Publisher Site | Google Scholar
J. B. Gutierrez, M. K. Hurdal, R. D. Parshad, and J. L. Teem, “Analysis of the Trojan Y chromosome model for eradication of invasive species in a dendritic riverine system,” Journal of Mathematical Biology, vol. 64, pp. 319–340, 2012.View at: Publisher Site | Google Scholar
X. Wang, J. R. Walton, R. D. Parshad, K. Storey, and M. Boggess, “Analysis of the Trojan Y-Chromosome eradication strategy for an invasive species,” Journal of Mathematical Biology, vol. 68, no. 7, pp. 1731–1756, 2014.View at: Publisher Site | Google Scholar
J. Lyu, P. Schofield, M. Beauregard, and R. D. Parshad, “A comparison of the trojan Y chromosome strategy to harvesting models for eradication of non-native species,” Natural Resource Modeling, vol. 33, no. 2, Article ID e12252, 2020.View at: Google Scholar
M. Jiang, X. Wu, K. Chen et al., “Production of YY supermale and XY physiological female common carp for potential eradication of this invasive species,” Journal of the World Aquaculture Society, vol. 49, no. 2, pp. 315–327, 2018.View at: Google Scholar
D. J. Schill, J. A. Heindel, M. R. Campbell, K. A. Meyer, and E. R. Mamer, “Production of a YY male brook trout broodstock for potential eradication of undesired brook trout populations,” North American Journal of Aquaculture, vol. 78, no. 1, pp. 72–83, 2016.View at: Google Scholar
J. L. Teem, J. B. Gutierrez, and R. D. Parshad, “A comparison of the Trojan Y Chromosome and daughterless carp eradication strategies,” Biological Invasions, vol. 16, pp. 1217–1230, 2014.View at: Publisher Site | Google Scholar
CDC, “Centers for disease control and prevention,” 2021, https://www.cdc.gov/mosquitoes/about/what-is-a-mosquito.html.View at: Google Scholar
M. A. Saleem and I. LobanovaA. I. Qureshi and O. Saeed, “Chapter 5-mosquito-borne diseases,” Dengue Virus Disease, Academic Press, Massachusetts, MA, USA, 2020.View at: Google Scholar
J. K. Biedler, B. A. Hall, X. Jiang, and J. ZhijianZ. Adelman, “Chapter 10 - exploring the sex-determination pathway for control of mosquito-borne infectious diseases,” Genetic Control of Malaria and Dengue, Academic Press, Massachusetts, MA, USA, 2016.View at: Google Scholar
W. H. Fleming and R. W. Rishel, Deterministic and Stochastic Optimal Control, Springer-Verlag, New York, NY, USA, 1975.
J. Macki and A. Strauss, Intoduction to Optimal Control Theory, Springer-Verlag, New York, NY, USA, 1982.
S. Lenhart and J. T. Workman, Optimal Control Applied to Biological Models, CRC Press, Boca Raton, FL, USA, 2007.