Difference between revisions of "Article:Current Visceral Leishmaniasis Research: A Research Review to Inspire Future Study (6076917)"
(Created page with "{{articleheader  Wikidata_code = Q57176530  title = Current Visceral Leishmaniasis Research: A Research Review to Inspire Future Study  publication_date = 20180710  auth...") 
m (Protected "Article:Current Visceral Leishmaniasis Research: A Research Review to Inspire Future Study (6076917)" ([Edit=Allow only administrators] (indefinite))) 
(No difference)

Latest revision as of 11:45, 15 April 2019
This page is the ScienceSource HTML version of the scholarly article described at https://www.wikidata.org/wiki/Q57176530. Its title is Current Visceral Leishmaniasis Research: A Research Review to Inspire Future Study and the publication date was 20180710. The initial author is Kaiming Bi.
Fuller metadata can be found in the Wikidata link, which lists all authors, and may have detailed items for some or all of them. There is further information on the article in the footer below. This page is a reference version, and is protected against editing.
Journal Information
Title: BioMed Research International
Current Visceral Leishmaniasis Research: A Research Review to Inspire Future Study
 Kaiming Bi
 Yuyang Chen
 Songnian Zhao
 Yan Kuang
 ChihHang John Wu
Department of Industrial and Manufacturing Systems Engineering, Kansas State University, 2061 Rathbone Hall., Manhattan, KS 66506, USA
Publication date (collection): /2018
Publication date (epub): 7/2018
Abstract
Visceral leishmaniasis (VL), one of the deadliest parasitic diseases in the world, causes more than 50,000 human deaths each year and afflicts millions of people throughout South America, East Africa, South Asia, and Mediterranean Region. In 2015 the World Health Organization classified VL as a neglected tropical disease (NTD), prompting concentrated study of the VL epidemic using mathematical and simulation models. This paper reviews literature related to prevalence and prevention control strategies. More than thirty current research works were reviewed and classified based on VL epidemic study methods, including modeling approaches, control strategies, and simulation techniques since 2013. A summarization of these technical methods, major findings, and contributions from existing works revealed that VL epidemic research efforts must improve in the areas of validating and verifying VL mathematical models with realworld epidemic data. In addition, more dynamic disease control strategies must be explored and advanced simulation techniques must be used to predict VL pandemics.
Paper
1. Introduction
Visceral leishmaniasis (VL), or kalaazar, is a protozoan disease that, second only to malaria in numbers of fatalities, afflicts millions of people worldwide [^{[1]}]. VL is primarily distributed in East Africa, South Asia, South America, and Mediterranean Region, with an estimated 50,000 to 90,000 new VL cases each year. Ninety percent of reported VL cases occur in Brazil, Ethiopia, India, Kenya, Somalia, South Sudan, and Sudan [^{[2]}]. In 2015 the World Health Organization (WHO) classified VL as a neglected tropical disease (NTD) due to relatively minimal granted attention from the public, resulting in high mortality rates (more than 20,000 in 2015), and endemic spreading in povertystricken regions around the world [^{[3]}, ^{[4]}].
VL is one of the most widespread human diseases, with more than 20 Leishmania species identified worldwide [^{[5]}]. Unlike other common forms of leishmaniasis, such as cutaneous leishmaniasis (CL), mucocutaneous leishmaniasis (ML), and postkalaazar dermal leishmaniasis (PKDL), VL symptoms usually occur internally [^{[1]}], meaning VL is more difficult to detect and cure than other leishmaniasis. Based on different kinds of susceptible species, VL can be classified as anthroponotic visceral leishmaniasis (AVL) or zoonotic visceral leishmaniasis (ZVL). AVL, which is transmitted between humans via vector carriers, is primarily caused by L. donovani throughout East Africa and the Middle East, especially Sudan, Somalia, Yemen, and Saudi Arabia [^{[3]}]. Since most VL vaccinations for humans are still in research or minimally effective [^{[6]}], AVL preventions include only early detection and treatment or use of longterm insecticide nets. ZVL, however, is transmitted between humans and other mammals, such as dogs. L. donovani, L. infantum, and L. archibaldi can result in ZVL, with specific concentrations in East Africa, South America, the Mediterranean basin, and South Asia [^{[6]}]. Because dogs are the most common mammal carriers of ZVL, ZVL control strategies such as culling dogs, dog vaccinations, and use of insecticide collars are more prevalent than any AVL control strategies [^{[7]}–^{[10]}].
The use of mathematical models to describe and predict epidemic transmissions has become a recent trend in disease research area. Mathematical models intuitively exhibit complex VL transmission processes, and they measure variables and system parameters to reveal VL spreading dynamics and related dominating factors. Rapid advancements in computer technologies have resulted in computeraided simulation that helps mathematical models directly predict future VL prevalence. Using results from model analysis, parametric estimation, and simulation experiments, researchers can study and anticipate disease transmission dynamics and identify disease control strategies to fight a VL pandemic. Consequently, an increasing number of studies have focused on mathematical modeling and corresponding analysis for VL disease dynamics. Approaches used in these studies can be generally categorized as system dynamic models, including ordinary differential equation (ODE) or partial differential equation (PDE) models, as well as statistic models, or machine learning models. The main contributions and results are concentrated in precise prediction tested by validation, determining the key parameters by sensitive analysis and analyzing the bifurcation point of the disease reproduction number.
A welldefined mathematical model can be used to develop disease control strategies that are ascertained by solving the mathematical model or using numerical experiments. Numerical control strategies are robust and reliable approaches because potential bias from empirical data is not included. Conversely, implementation of reallife control strategies can be cost prohibitive, irreversible, and difficult to apply in a largescale format, especially for developing countries that lack public health resources. However, computeraided simulations that compare possible control strategies derived from a mathematical model can be carried out, and they are relatively inexpensive and can be performed repeatedly to examine system sensitivity and determine optimal control parameters. Almost half of corresponding research used mathematical modeling approaches to study potential disease control strategies, including dog culling, use of insecticidal dog collars, vector controls, and insecticide spraying strategies. Using optimal control, parametric analysis, or stochastic control methods, research results provided welldeveloped guidelines for disease control centers to prevent or mitigate a VL pandemic.
The rest of this paper is organized into comprehensive sections. Section 2 introduces current worldwide VL pandemic situations, and Section 3 presents VL control strategies in existing literature and recent breakthroughs in this field. Section 4 reviews papers on VL mathematical modeling and summarizes new developments and significant contributions since 2013. Section 5 reviews papers on control strategies and the use of numerical simulation, and Section 6 concludes the paper and suggests possible areas for future VL pandemic research. The research tree of this paper is shown in Figure 1.
2. VL Pandemic Recent History
Pigott et al. collected and summarized prevalence data reports of CL and VL from 1960 to 2012 [^{[16]}]. According to their summarized database, the worldwide VL pandemic has affected at least 55 countries and territories (Figure 2), and another 45 countries have reported unspecified and borderline VL cases. Affected countries are located primarily in Latin America, the Mediterranean basin, Northeastern Africa, and Asia. Moreover, VL outbreaks historically occur most often in developing or agricultural countries because citizens are more likely to be in contact with disease transmission vectors such as dogs and mosquitoes [^{[17]}]. Worse still, VL traps many families in these developing countries into vicious circle. Families affected by VL have to spend more direct cost (like treatment) and indirect cost (like loss of household income) during the VL epidemic.
According to the WHO neglected tropical diseases (NTD) report in 2007, VL is identified as one of the NTDs [^{[18]}]. The primary reason WHO classified VL as an NTD is that confirmed VL cases have decreased worldwide from approximately 60,000 to around 20,000 since around 2010 [^{[11]}], as shown in Figure 3. However, thousands of VL cases may not be covered in the WHO VL estimation report [^{[19]}] since some countries without public health information systems do not submit their infection report data to WHO (e.g., Chad, the Central African Republic, and the Democratic Republic of the Congo). Even though for the countries with completed public health information system, the reported epidemiological data could be incomplete, and official figures are likely to underestimate grossly [^{[8]}]. Based on the estimation, there are 500,000 new cases of VL and more than 50,000 deaths from the VL each year [^{[17]}]. In 2012, another research group from the WHO Leishmaniasis Control Team corrected the number of VL estimation cases into 0.20.4 million, and the number of VL deaths into 20,00040,000 [^{[20]}]. Therefore, VL continues to be a serious and present threat to public health worldwide.
By observing the VL epidemiological situations for each country in Figure 4 based on the report from WHO [^{[11]}], the significant drop in reported VL cases can primarily be attributed to a significant decrease in reported VL cases in India and Bangladesh. Between 2006 and 2016 the number of reported VL cases in India decreased from 39,173 to 6,249 and reported cases in Bangladesh declined from 9,379 to 255. The major reason of this decreasing is the wide utilization of insecticidetreated nets [^{[43]}, ^{[44]}]. However, VL prevalence did not change significantly for other countries. For example, the number of annually reported VL cases in Brazil were approximately 2,700 throughout the 10year reporting period.
Conversely, reported VL cases in several African countries have shown significant increase since 2006, as shown in Figure 5. For example, Somalia reported less than 100 cases in 2006 and 781 cases in 2016. Although the population of Somalia (14.32 million) is much less than the population of India (13.24 billion), the proportions of VL cases reported in these two countries were similar in 2016. If no immediate measures are taken in these selected African countries, largescale VL outbreaks are imminent.
Therefore, VL is still a serious disease which threatens people lives and health especially in the developing countries. To strain the transmission of VL, WHO and health organizations in VL afflicted countries should apply effective prevention and control strategies. In the next section, this paper will review several existing strategies against VL.
3. VL Mitigation and Prevention Methods
Since 1995, researchers have focused on ZVL when investigating the intervention and prevention of VL transmission because many ZVL control strategies are related to animals. Tesh categorized former ZVL control strategies into three main classes: early detection of human cases, destruction or treatment of infected dogs, and vector control [^{[10]}]. Tesh's paper mentioned that treating infected people did not affect the parasite transmission and the drug resistance of L. infantum made the treatment even more difficult. For the dog control, many serologically positive dogs are harddetectible. And the VL infected dogs need expensive and continual treatments. Tesh's paper also proposed that preventing the disease in dogs population is the best control strategies, since it can interrupt the transmission cycle of ZVL. In 2002, Guerin et al. asserted that principal interventions of VL via early diagnosis and treatments, dog controls, and vector controls are effective control strategies [^{[8]}]. However, Guerin et al. also pointed out that vector controls such as residualinsecticide house spraying are cost prohibitive and rarely used. Especially in India, Phlebotomus argentipes (the dominated species of sand fly) is becoming resistant to the insecticide. The authors also mentioned some special challenges in VL endemic areas: lack of financial support (India and Bangladesh), remote places (Brail), and civil war (Sudan). In 2006, DantasTorres et al. introduced the Brazilian Leishmaniasis Control Program (BLCP), which includes diagnosis and early treatment of human cases, immunological screening of seropositive dogs, and insecticide spraying [^{[7]}]. DantasTorres et al. indicated that the culling of seropositive dogs has limited effects, as proven by experiments and mathematical models. Authors summarized several key points of unsuccessful culling dogs: the limitation of the immunological screenings to detect antiLeishmania antibodies, the opposition of dog owners to the culling of asymptomatic dogs, and the lack of evidence that it is an effective intervention strategy. A paper by Quinnell et al. suggested that tropical insecticidescollars and pourons can be used to reduce VL incidences for dogs by more than 83% [^{[9]}]. However, like the insecticide spraying, the cost of this strategy limits the feasibility of the tropical insecticidescollars. The authors offered an alternative method, insecticidal bath, which can protect dogs for at least 3.5 months against Phlebotomus chinensis.
In 2014, Werneck considered the effectiveness of control strategies on the basic reproduction number R0 [^{[22]}]. The author found that vector controls (e.g., controlling vector density, the ratio of female vector, and the extrinsic incubation period of L. infantum in sand flies) and dog controls (e.g., culling infected dogs, dog vaccinations, and insecticidereleasing dog collars) can cause the disease reproduction rate R0 to decrease to less than 1, meaning that the number of infected individuals will eventually decrease to zero. However, the author did not find the relevant data to support these control strategies. He pointed out some potential implementation difficulties for these strategies, such as costs issues associated with the continual uses of tropical insecticidescollars. Werneck also voiced concern that the effectiveness of using insecticide collars in the large community (like Brazil) may not work so well, since the insecticide collars strategy has the relatively shortterm effect and consequent need.
Due the high cost of indoor residual spraying, insecticidetreated nets (ITNs) were introduced as an alternative control strategy for ZVL [^{[21]}]. Experimental trials in subSaharan Africa, Latin American, Thailand, Pakistan, and Iran show that ITNs could reduce infections with malaria by 17%–62%. In 2015, Picado et al. summarized the results of the KALANET project to analyze ITN effectiveness. The KALANET project is a clusterrandomized controlled trial in India and Nepal [^{[23]}]. The KALANET project was conducted in 26 highincidence regions (in India and Nepal) with more than 20,000 inhabitants observed over 24 months. Results showed that use of longlasting insecticidal nets resulted in a 50% reduction of L. donovani infections. However, this study also showed that some of those nets were untreated, many were damaged, and most families did not have enough nets to protect all family members. The authors suggested standardizing the color and size of insecticide nets; they also want the government to replace the untreated and damaged nets by new treated nets for free for publics. Although the longlasting insecticidal nets can prevent VL transmissions while people are sleeping, recent entomological findings in India indicated that vectors are more exophilic than previously thought, meaning that people engaged in outdoor work have more opportunities to become infected [^{[45]}].
In 2016, Özbel et al. analyzed the geographical distribution, ecological aspects, and species habits of VL vectors (sand flies) in Bangladesh [^{[24]}]. In general, two genera of sand fly (Phlebotomus and Sergentomyia) and a total of 13 species were recorded in Bangladesh, among which Phlebotomus argentipes is the dominant vector species in VL endemic areas of Bangladesh (>94%). Researchers also found that the Phlebotomus. argentipes population peaks around monsoon season and reaches the lowest ebb during winter and summer in Bangladesh. They also determined that eight ecological parameters (soil temperature and moisture, rainfall, air temperature, relative humidity, soil pH, soil organic carbon, and wind speed) can influence the Phlebotomus. argentipes population. Future research must enact control strategies based on their ecological aspect. The potential application of this research can provide early warning of the incoming VL epidemic and narrow the range of the insecticide spraying.
As the most effective control strategy for infectious diseases, the successful vaccination on VL is longawaited for the VL afflicted countries. The experiment on VL vaccination was started in 1990s; researchers tried to utilize the proteins from L. donovani to develop vaccines [^{[46]}]. At the beginning of the 21st century, researchers investigated the possibility of vaccination based on DNA and genetically attenuated parasites [^{[47]}, ^{[48]}]. More recently, Duthie et al. studied several different vaccine antigens for VL transmissions using recombinant proteins from E. coli [^{[25]}]. Their research results have shown that several potential vaccine antigen candidates are verified in different platforms. The authors consolidated seven vaccine candidates as recombinant proteins in E. coli, and they validated the effectiveness of E. coli to L. donovani via experiments on mice. However, their research pertained to only nonhuman experiments. Significant research and advancements are still needed to obtain effective vaccination for humans against VL parasites.
For now, even though many contributions have been done for the VL controls and prevention, an effective, feasible, and economical control strategy is still an ongoing effort. The current VL control strategies and corresponding deficiencies are summarized in Table 1. Since the most of VL severely afflicted countries are developing countries, how to balance the effectiveness and costs involved in such VL control plan is delicate tradeoff. This paper will discuss more about this particular issue in Section 5.
4. Mathematical Epidemiology Models for VL
4.1. System Dynamic Model
In 1996 Dye first introduced a 4equation ODE susceptiblelatentinfectiousremoved (SLIR) model to describe the VL epidemic [^{[49]}]. SLIR, respectively, represents the susceptible, latent, infectious, and resistant populations of VL, and the model considered the transitions between these populations. However, the author did not consider the behaviors of dog and vector in his model. Courtenay et al. improved Dye's model by considering the dog population as a key parameter in the SLIR model [^{[50]}]. Variations of this parameter directly influence the human infection rate. Ribas et al. built an ODE model to describe VL transmissions among humans, dogs, and vectors [^{[14]}]. The model utilized 11 ODE equations, including the susceptible (dog, sand fly, and human), latent (dog, sand fly, and human), infectious (dog, sand fly, and human), and recovered (dog and human) populations. Although their model was presented without data, simulation, and analysis, it was the first model to describe behaviors for all species involved in VL.
Since WHO's designation of VL as an NTD in 2015 [^{[51]}], an abundance of research and studies have focused on developing mathematical models of VL. In 2016 Zhao et al. developed a 12equation ODE model to comprehensively describe a VL epidemic (Figure 6) [^{[12]}]. They improved the model by including a hospitalized population into the ODE system; this population has a higher probability of survival than infections without hospitalization due to the systematic treatment. The authors utilized a backward bifurcation method to analyze VL equilibrium behavior and the basic reproduction rate R0. Their research showed that VL equilibrium is highly related to a critical model parameter, Rc, as the epidemic threshold value for R0. Similarly, Subramanian et al. proposed a compartmentbased ODE model of VL transmission to explain disease transmissions in symptomatic VL, asymptomatic VL, and PKDLinfection classes [^{[26]}]. Sensitivity analysis of model parameters found that the biting and birthing rates of sandflies and the recovery rate of symptomatic humans are dominating factors for VL epidemic control.
Biswas simplified the 12equation ODE model from Zhao et al. to an 8equation ODE model by dividing the nonhuman populations (dog and vector) into susceptible and infected population groups [^{[27]}]. With a simplified ODE model, researchers reduced the complexity of system sensitivity analysis and reduced the numbers of assumed or estimated model parameters. This model also successfully reproduced the 2011 VL epidemic in South Sudan. Shimozako et al. transferred the ODE model in Ribas et al. by considering the dog population as the only source of infection since vectors could not transfer ZVL without dogs [^{[28]}]. Therefore, their mathematical model contained eight variables corresponding to the susceptible, latent, infectious, and recovery populations for dogs and humans. Le Rutte et al. compared three ODE models and corresponding simulation results while considering indoor residual spraying [^{[29]}]. The primary difference between the three models was how relationships between PKDL and the recovery population were modeled. Their research predicted that, using 60%–80% IRS coverage, VL epidemics could be eliminated within three years in Bihar, India. In addition, other researchers have made incremental contributions using various ODE VL epidemic models [^{[30]}–^{[33]}].
Although many ODE models describe VL epidemics and transmission, the development of a novel dynamic model is an active area of research in the investigation of complicated transmission behaviors of VL under various situations and the development of improved mitigation and control strategies. Bi et al. introduced a twodimensional PDE model based on an existing ODE model [^{[13]}]. The model presented in their research considered human age structures and time as two dimensions since historical data have strongly suggested that VL infection rates are highly correlated to human age groups. For example, children and teenagers (aged from 0 to 20) are more likely to become infected compared to other age groups. This research used computer simulation and mathematical analysis to explain this phenomenon and recover the historical VL endemic data published by WHO.
4.2. Models Based on RealWorld Data
VL attracted significant epidemiology research attention; abundant statistical data were collected and reported on the current VL pandemic worldwide by scholars and researchers. Many researchers realized the importance of data utilization in VL model development. Disease data is generally utilized in three ways: use of reported data to build statistical models, use of historical data to predict future prevalence, and use of existing data to calibrate model parameters in mathematical epidemic models.
The primary objective of building a VL statistical model is to statistically identify key parameters in the VL transmission process and determine relationships between the number of parameters and the number of infected population. Werneck et al. used consolidated census tracts to analyze VL disease prevalence data from different regions of Brazil [^{[52]}]. The authors developed a spherical covariance structure model based on census data from 1992 to 1996 in Teresina, Brazil. By exploring spatial correlation structures of the census data, they found a positive correlation between reported VL incidences and residence in areas of green vegetation, especially in favela. In 2007 Werneck et al. extended their previous work by analyzing and comparing the results from 21 statistical models [^{[53]}]. According to human and canine VL data in Teresina, Brazil, the study found significant correlations between residence in areas with green vegetation and infected dogs and between the human infection rate and urbanization index or socioeconomic status index.
Thompson et al. studied relationships between climate and VL epidemics by establishing the statistical regression model [^{[54]}]. Their research found that rainfall is the most significant parameter statistically correlated to VL incidences. The influence of geographical features of areas of residence (e.g., cities, plains, or foothills) on VL transmission was also considered. The foothill population statistically revealed a higher risk of VL infection than other populations. De Araújo et al. considered the statistical model with data and found that spatial data is more reliable and accurate for VL epidemic study and analysis [^{[34]}]. Therefore, they applied spatial statistical modeling and the Bayesian approach to model and estimate risks of VL using historical data from Belo Horizonte, Brail. Their research showed that the relative risk of VL is correlated to income, education, and the ratio of infected dogs to inhabitants in Belo Horizonte. However, as opposed to their previous research findings, residence in areas with green vegetation did not show significant correlations to the risk of VL.
Ecological niche modeling (ENM), stemming from the genetic algorithm [^{[55]}], has been widely used to predict VL prevalence since 2006. Nieto et al. first used ENM to analyze VL data from northeastern Brazil [^{[56]}]. Using the geographic information system (GIS), the ENM model can predict prevalence risks in three levels (high, moderate, and low). When being validated with historical data from Bahia, Brazil, the predictive model demonstrated high accuracy (more than 90%) on highlevel and moderatelevel data. Similar approaches were used to predict VL prevalence in North America and Middle Eastern regions [^{[57]}, ^{[58]}]. Several other methods were utilized to predict the trend of VL epidemics with various factors. Elnaiem et al. summarized data from eastern Sudan in GIS and then used that data to build predictive models of VL infections based on rainfall and corresponding distance to a river [^{[59]}]. Oshaghi et al. built a predictive degreeday model for VL using the single triangulation method [^{[60]}]. This model predicted temporal and spatial distribution of VL infection density and generations of sandflies. KaragiannisVoules et al. employed Bayesian geostatistical models to fit the VL incidence data from Brazil, and they identified environmental and socioeconomic predictors using Bayesian variables [^{[15]}]. Their research results predicted that regions with humid climates and dense vegetation distributions are more vulnerable to VL than other regions.
Parameter estimation is another essential application when validating VL mathematical models using realworld data since the use of assumed system parameters in the model may reduce model reliability. Bi et al. summarized age structures of VL infections in various regions [^{[13]}]. They revised the constant infection rate into agedependent distribution of infection rates by studying historical VL prevalence data from existing literature from diverse regions over time (Figure 7). Use of this agedependent distribution function allowed their PDE epidemiology model to reflect VL human prevalence in both age and time periods. Biswas et al. established the posterior distribution of different parameters and initial parameters based on observation data [^{[27]}]. Use of parameter distribution allowed the simulation experiments to reflect more than one result with unique possibilities. Even though most mathematical models still use assumptions or estimations with existing literature as their system parameters [^{[12]}, ^{[26]}, ^{[28]}, ^{[29]}], an increasing number of studies are utilizing realworld data to more accurately estimate system parameters and validate their models.
5. Optimal Control Strategies Using VL Models
5.1. Parameter Control Strategy
The most generalized control strategy in a VL mathematical model is the parameter control strategy, which assumes that the key parameters in the model are adjustable. When the parameters are adjusted, the model outputs become dependent variables; therefore, the parameter adjustment process can be considered a corresponding realworld control strategy. In 2002 Courtenay et al. introduced the numerical control strategy to the field of VL mathematical modeling [^{[50]}]. They assumed that the parameter of dog density could be controlled by culling dogs in specific areas, and their research findings showed that dog culling could effectively reduce the proportion of infectious population. Ribas et al. improved the parameter control method by considering additional control parameters [^{[14]}], including the dog treatment rate, the insecticide collar utilization rate, the dog culling rate (natural mortality rate of dogs), the dog vaccination rate, and the mortality rate of sandflies (vector control). Using the control parameter method, they used simulation experiments to compare the efficiencies of each control strategy, as shown in Figure 9, where vector control proved to be more effective than control strategies such as culling and vaccine. Shimozako et al. incorporated control strategies such as dog treatment, dog vaccination, and use of insecticidal dog collars into an ODE mathematical model of VL [^{[37]}]. However, they increased the dimension of the model by using control variables to replace the control parameters in their model. A costeffect analysis and simulation experiments showed that use of insecticidal dog collars should be the most utilized control strategy.
5.2. Optimal Control Strategy
Lev Pontryagin established the optimal control strategy in the 1950s [^{[61]}, ^{[62]}]. This strategy provides optimality criterion by maximizing or minimizing a given objective function subject to constraints defined in the differential equation mathematical model. Zhao et al. introduced optimal control into their 12equation ODE mathematical model of VL [^{[12]}], including the susceptible, latent, and infectious populations for sandflies, humans, and dogs; recovered human and dog populations; and hospitalized human population. Their study included three control strategies in the model: the control dog prevention level (vaccine protection or dog culling), the control insecticide usage level (insecticide sprayed around buildings), and the control personal protection level (longlasting insecticide). Consideration of exposed human, infected human, and sandfly populations as well as control cost as the control objective function caused the final control strategy to effectively lower the negative influence of VL by approximately 80% (Figure 9). Agusto et al. considered the use of fabrics and insecticidal animal collars as an additional control [^{[33]}]. Their research was innovative because they used various combinations of controls instead of utilizing control strategies identical to previous research. Their simulated experimental results showed that use of an ODE mathematical model allowed their control strategies to effectively reduce VL and PKDL infected human populations. Biswas et al. revised a previous optimal control by suggesting a strategies mechanism and incorporating various combinations of optimal strategies into the ODE mathematical model with a given objective function [^{[36]}]. This objective function was a linear combination of VL and PKDL infections with the cost of control strategies. They calculated the corresponding infection averted ratios (IAR) and incremental costeffectiveness ratios (ICER) for each control strategy, where IAR is the ratio of the number of infections averted to the number of recoveries and ICER is the additional cost per additional health outcome. The researchers then selected the strategy with the highest IAR and lowest ICER as the optimal strategy. An optimal control strategy specifically targeting the human VL vaccination was also analyzed, but no evidence revealed the effectiveness of a human VL vaccination [^{[27]}].
5.3. Control Strategy Selection Using Simulation
Simulation comparison is the most common method of VL mathematic control modeling in simulation. Figure 8 compares the human infected population with control (dog culling, sandfly control, and human protection) and without control [^{[12]}]. The solid line in the figure shows epidemic performance (which includes the exposed human Eh, the infected human Ih, and the total population of sand flies Nf) of the control strategy, which reduced disease prevalence to 80%. This simulation proved the effectiveness of the combined control strategies.
Efficacy comparison between control strategies is another type of simulation comparison. Ribas et al. compared the human prevalence influence using vector control, insecticidal collars, dog culling, dog vaccines, and dog treatment [^{[14]}], as shown in Figure 9. When the parameter was changed to manipulate control levels, the human prevalence by vector control decreased most rapidly, proving that human prevalence is the most sensitive to vector controls. Therefore, vector control is the most effective control strategy.
Spatial simulation is a simulation estimation method that provides spatial information throughout the model behaviors. Using GIS, spatial simulation can exhibit VL prevalence information from various regions, as shown in Figure 10. KaragiannisVoules et al. utilized historical data (2001–2009) from Brazil to build a statistical model [^{[15]}]. Their simulation used GIS to predict VL prevalence in Brazil in 2010, thereby directly reflecting high infection density in the eastern region of Brazil.
6. Discussion and Conclusion
Although globally reported, VLconfirmed cases have decreased since 2011; VL prevalence has not improved significantly worldwide except in South Asia (i.e., India and Bangladesh). However, VL outbreaks have increased in Ethiopia, Somalia, and Kenya since around 2008. Moreover, public health agencies in underdeveloped African countries such as Chad and the Central African Republic do not have resources and capabilities to collect and report VL incidences. However, because these countries are near regions severely afflicted with VL, such as Sudan and South Sudan, the number of global VL cases reported from WHO may be underestimated.
This paper reviewed current VL epidemiological research ranging from VL epidemic control strategies to VL mathematical models and related optimal control strategies. The research demonstrated how to use numerical methods such as modeling and sensitivity analysis, as well as equilibrium/stability studies and simulation experiments, to assist mitigation and prevention strategies for a worldwide VL pandemic. Governments and health organizations can utilize the modeling and simulation results to predict or estimate impacts of various control strategies.
Despite significant research efforts using mathematical models for the VL epidemic, research gaps still exist and many areas of study remain unexplored. Table 2 summarizes 21 literature works related to mathematical modeling and VL disease control strategies since 2013. Most of the reviewed research used or proposed system dynamic mathematical models or statistical models; approximately half of these research works considered realworld data and studied possible control strategies. Only one paper included realworld data, control strategies, mathematical models, and numerical simulation experiments.
Future work, thorough VL epidemic research using mathematical or statistical models, ought to consider the four following aspects:
building more sophisticated mathematical models to explain underlying infectious disease transmission dynamics,
including realworld data to aid model validation and verifications,
exploring possible disease control/mitigation strategies to increase understanding of model maneuverability and robustness,
using numerical simulation experiments as a predictive tool to verify the feasibility of model and control strategies.
Moreover, future work in these four aspects of VL mathematical modeling must utilize modern analytical tools. The disadvantage of current modeling is the limited diversity of model types. A majority of existing VL mathematical models are ODE models, which are widely used but produce limited predicted results without details. Therefore, more statistical, machine learning, and PDE models are needed to build sophisticated, comprehensive mathematical models of VL. Statistical and machine learning models can more advantageously utilize realworld data to ensure model prediction accuracy, while use of a PDE model can enrich predicted results with age, gender, socioeconomic group, ethics, and spatial information. For the second aspect, the inclusion of realworld data, most test data currently used to validate and verify underlying mathematical models are estimated or assumed, consequently limiting the mathematical model to reflect only data from previous VL epidemic episodes. Future research efforts should utilize recent epidemic data with temporal and spatial data during the modeling phase, making the modeling process increasingly dynamic and reflecting realtime data while predicting possible trends of an ongoing epidemic. The current primary disadvantage of the third aspect, exploring possible control strategies, is that the control strategies lack of applicability in the real world. In fact, the most effective control strategies suggested by the mathematical models may not be operable or they may be too cost prohibitive to be implemented. Operable control strategies should be carefully quantized, such as specific consideration of the optimal level of canine culling in a particular time frame or the level of insecticide spraying in each area affected by VL. For the fourth aspect, current studies using numerical simulation experiments frequently provide insufficient information from simulation results. Most simulations of VL models can only predict the trend of VL infections. Future research should focus on spatial simulation and agentbased simulation as well as the study of the interactions between multiple regions or environments.
In conclusion, the use of mathematical models to study, analyze, and predict VL epidemics and to explore effective and implementable control strategies remains an active and studyworthy area of future research. However, research results from more comprehensive studies that use modern analytical tools will help public health organizations understand and prevent the VL disease.
References
 ↑ W. AlSalemJ. R. HerricksP. J. HotezA review of visceral leishmaniasis during the conflict in South Sudan and the consequences for East African countriesParasites & Vectors201691, article no. 4602s2.08498335656610.1186/s130710161743727549162
 ↑ J. AlvarS. YactayoC. BernLeishmaniasis and povertyTrends in Parasitology200622125525572s2.03375084532210.1016/j.pt.2006.09.00417023215
 ↑ J. A. R. PostigoLeishmaniasis in the world health organization eastern mediterranean regionInternational Journal of Antimicrobial Agents2010361S62S6510.1016/j.ijantimicag.2010.06.0232s2.07795793724220728317
 ↑ W. H. OrganizationWorking to overcome the global impact of neglected tropical diseases: first WHO report on neglected tropical diseases, 2010
 ↑ A. OryanM. AkbariWorldwide risk factors in leishmaniasisAsian Pacific Journal of Tropical Medicine20169109259322s2.08499273146710.1016/j.apjtm.2016.06.02127794384
 ↑ K. JainN. K. JainVaccines for visceral leishmaniasis: A reviewJournal of Immunological Methods20154221122s2.08493043502510.1016/j.jim.2015.03.01725858230
 ↑ F. DantasTorresS. P. BrandãoFilhoVisceral leishmaniasis in Brazil: Revisiting paradigms of epidemiology and controlRevista do Instituto de Medicina Tropical de São Paulo20064831511562s2.03374597549210.1590/S00364665200600030000716847505
 ↑ P. J. GuerinP. OlliaroS. SundarVisceral leishmaniasis: current status of control, diagnosis, and treatment, and a proposed research and development agendaThe Lancet Infectious Diseases20022849450110.1016/s14733099(02)00347x2s2.0003631971212150849
 ↑ R. J. QuinnellO. CourtenayTransmission, reservoir hosts and control of zoonotic visceral leishmaniasisParasitology200913614191519342s2.07744910481810.1017/S003118200999115619835643
 ↑ R. B. TeshControl of zoonotic visceral leishmaniasis: is it time to change strategies?The American Journal of Tropical Medicine and Hygiene19955232872922s2.0002891455210.4269/ajtmh.1995.52.2877694973
 ↑ WHO.World Health Organization. Leishmansishttp://www.who.int/leishmaniasis/burden/en/
 ↑ S. ZhaoY. KuangC.H. WuD. BenAriehM. RamalhoOrtigaoK. BiZoonotic visceral leishmaniasis transmission: modeling, backward bifurcation, and optimal controlJournal of Mathematical Biology201673671525156010.1007/s002850160999zMR3556842Zbl1350.920592s2.08496374714727084184
 ↑ K. BiY. ChenS. ZhaoC.H. WuA New Zoonotic Viseral Leishmaniasis Dynamic Transmission Model with AgeStructurePLOS Neglected Tropical Diseases2018
 ↑ L. M. RibasV. L. ZaherH. J. ShimozakoE. MassadEstimating the optimal control of zoonotic visceral leishmaniasis by the use of a mathematical modelThe Scientific World Journal2013
 ↑ D.A. KaragiannisVoulesR. G. C. ScholteL. H. GuimarãesJ. UtzingerP. VounatsouBayesian Geostatistical Modeling of Leishmaniasis Incidence in BrazilPLOS Neglected Tropical Diseases2013752s2.084878498943
 ↑ D. M. PigottN. GoldingJ. P. MessinaGlobal database of leishmaniasis occurrence locations, 19602012Scientific Data201412s2.08495953283310.1038/sdata.2014.3614003625984344
 ↑ F. ChappuisS. SundarA. HailuVisceral leishmaniasis: what are the needs for diagnosis, treatment and control?Nature Reviews Microbiology200751187388210.1038/nrmicro17482s2.03534887102817938629
 ↑ B. ThimphuReview of The Decisions And Resolutions of The Sixtieth World Health Assembly And The 120Th And 121St Sessions of The Executive Board
 ↑ W. H. OrganizationReport of a WHO informal consultation on liposomal amphotericin B in the treatment of visceral leishmaniasis, Rome, Italy, 2007
 ↑ J. AlvarI. D. VélezC. BernLeishmaniasis worldwide and global estimates of its incidencePLoS ONE20127510.1371/journal.pone.0035671e356712s2.084861665791
 ↑ C. LengelerInsecticidetreated bed nets and curtains for preventing malariaCochrane Database of Systematic Reviews20042CD0003632s2.03843070027
 ↑ G. L. WerneckVisceral leishmaniasis in Brazil: Rationale and concerns related to reservoir controlRevista de Saúde Pública20144858518562s2.08491212252110.1590/S00348910.201404800561525372177
 ↑ A. PicadoB. OstynS. RijalLonglasting Insecticidal Nets to Prevent Visceral Leishmaniasis in the Indian Subcontinent; Methodological Lessons Learned from a Cluster Randomised Controlled TrialPLOS Neglected Tropical Diseases2015942s2.084929462820
 ↑ Y. ÖzbelC. SanjobaY. MatsumotoGeographical distribution and ecological aspect of sand fly species in BangladeshKala Azar in South Asia2016199209
 ↑ M. S. DuthieM. FavilaK. A. HofmeyerStrategic evaluation of vaccine candidate antigens for the prevention of Visceral LeishmaniasisVaccine20163425277927862s2.08496482837810.1016/j.vaccine.2016.04.06727142329
 ↑ Y. WangH. WeiJ. WangElectropolymerized polyaniline/manganese iron oxide hybrids with an enhanced color switching response and electrochemical energy storageJournal of Materials Chemistry A2015341207782079010.1039/c5ta04439a
 ↑ S. BiswasMathematical modeling of visceral leishmaniasis and control strategiesChaos, Solitons & Fractals201710454655610.1016/j.chaos.2017.09.005MR3715120
 ↑ H. J. ShimozakoJ. WuE. MassadMathematical modelling for Zoonotic Visceral Leishmaniasis dynamics: a new analysis considering updated parameters and notified human Brazilian dataInfectious Disease Modelling20172214316010.1016/j.idm.2017.03.00229928734
 ↑ E. A. Le RutteL. A. C. ChapmanL. E. CoffengElimination of visceral leishmaniasis in the Indian subcontinent: a comparison of predictions from three transmission modelsEpidemics20171867802s2.08501448504810.1016/j.epidem.2017.01.00228279458
 ↑ D. N. C. C. CostaC. T. CodeçoM. A. SilvaG. L. WerneckCulling Dogs in Scenarios of Imperfect Control: Realistic Impact on the Prevalence of Canine Visceral LeishmaniasisPLOS Neglected Tropical Diseases2013782s2.08488344322810.1371/journal.pntd.0002355e235523951375
 ↑ A. P. SeváF. G. OvallosM. AmakuCanineBased Strategies for Prevention and Control of Visceral Leishmaniasis in BrazilPLoS ONE2016117p. e01600582s2.08501456069010.1371/journal.pone.016005827471852
 ↑ L. ZouJ. ChenS. RuanModeling and analyzing the transmission dynamics of visceral leishmaniasisMathematical Biosciences and Engineering20171456158516042s2.08502053521210.3934/mbe.2017082Zbl1365.9214029161877
 ↑ F. B. AgustoI. M. ELmojtabaOptimal control and costeffective analysis of malaria/visceral leishmaniasis coinfectionPLoS ONE20171222s2.08501171136910.1371/journal.pone.0171102e0171102
 ↑ V. E. M. de AraújoL. C. PinheiroM. C. M. de AlmeidaRelative risk of visceral leishmaniasis in Brazil: a spatial analysis in urban areaPLOS Neglected Tropical Diseases201371110.1371/journal.pntd.0002540e25402s2.084888222929
 ↑ E. MillerA. WarburgI. NovikovQuantifying the Contribution of Hosts with Different Parasite Concentrations to the Transmission of Visceral Leishmaniasis in EthiopiaPLOS Neglected Tropical Diseases20148102s2.08492050704410.1371/journal.pntd.000328825356795
 ↑ S. BiswasA. SubramanianI. M. ELMojtabaJ. ChattopadhyayR. R. SarkarOptimal combinations of control strategies and costeffective analysis for visceral leishmaniasis disease transmissionPLoS ONE20171222s2.08501346059810.1371/journal.pone.0172465e017246528222162
 ↑ H. J. ShimozakoJ. WuE. MassadThe Preventive Control of Zoonotic Visceral Leishmaniasis: Efficacy and Economic EvaluationComputational and mathematical methods in medicine2017
 ↑ A. StauchH.P. DuerrA. PicadoModelBased Investigations of Different VectorRelated Intervention Strategies to Eliminate Visceral Leishmaniasis on the Indian SubcontinentPLOS Neglected Tropical Diseases2014842s2.08490128419110.1371/journal.pntd.0002810e281024762676
 ↑ M. ZamirG. ZamanA. S. AlshomraniControl strategies and sensitivity analysis of anthroponotic visceral leishmaniasis modelJournal of Biological Dynamics201711132333810.1080/17513758.2017.133983528648137
 ↑ F. BoukhalfaM. HelalA. LakmecheMathematical Analysis of Visceral Leishmaniasis Model. Research in Applied Mathematics20171
 ↑ K. K. GorahavaJ. M. RosenbergerA. MubayiOptimizing insecticide allocation strategies based on houses and livestock shelters for visceral leishmaniasis control in Bihar, IndiaThe American Journal of Tropical Medicine and Hygiene20159311141222s2.08493699323410.4269/ajtmh.14061225940194
 ↑ K. S. RockR. J. QuinnellG. F. MedleyO. CourtenayProgress in the Mathematical Modelling of Visceral LeishmaniasisAdvances in Parasitology201694491312s2.08499236598410.1016/bs.apar.2016.08.00127756459
 ↑ B. OstynV. VanlerbergheA. PicadoVector control by insecticidetreated nets in the fight against visceral leishmaniasis in the Indian subcontinent, what is the evidence?Tropical Medicine & International Health2008138107310852s2.04834908942310.1111/j.13653156.2008.02110.x18564350
 ↑ A. PicadoS. P. SinghS. RijalLonglasting insecticidal nets for prevention of Leishmania donovani infection in India and Nepal: paired cluster randomised trialBMJ (Clinical research ed.)2010341p. c676010.1136/bmj.c67602s2.078651394675
 ↑ R. M. PochéR. GarlapatiD.E. A. ElnaiemD. PerryD. PochéThe role of Palmyra palm trees (Borassus flabellifer) and sand fly distribution in northeastern IndiaJournal of Vector Ecology20123711481532s2.08486045265310.1111/j.19487134.2012.00211.x22548548
 ↑ C. L. JaffeN. RachamimR. SarfsteinCharacterization of two proteins from leishmania donovani and their use for vaccination against visceral leishmaniasisThe Journal of Immunology199014426997062s2.000250088362295807
 ↑ V. F. AmaralA. TevaM. P. OliveiraNetoStudy of the safety, immunogenicity and efficacy of attenuated and killed Leishmania (Leishmania) major vaccines in a rhesus monkey (Macaca mulatta) model of the human diseaseMemórias do Instituto Oswaldo Cruz2002977104110482s2.0003678322810.1590/S00740276200200070001912471434
 ↑ L. KedzierskiY. ZhuE. HandmanLeishmania vaccines: progress and problemsParasitology2006133S2S87S11210.1017/s00311820060018312s2.03384684777717274851
 ↑ C. DyeThe logic of visceral leishmaniasis controlThe American Journal of Tropical Medicine and Hygiene19965521251302s2.0002978218810.4269/ajtmh.1996.55.1258780448
 ↑ O. CourtenayR. J. QuinnellL. M. GarcezJ. J. ShawC. DyeInfectiousness in a cohort of Brazilian dogs: why culling fails to control visceral leishmaniasis in areas of high transmissionThe Journal of Infectious Diseases200218691314132010.1086/3443122s2.0003683828012402201
 ↑ W. H. OrganizationInvesting to Overcome the Global Impact of Neglected Tropical Diseases: Third Who Report on negLected Tropical Diseases 2015World Health Organization20153
 ↑ G. L. WerneckJ. H. MaguireSpatial modeling using mixed models: an ecologic study of visceral leishmaniasis in Teresina, Piauí State, BrazilCadernos de Saúde Pública200218363363710.1590/S0102311X20020003000072s2.0003658341912048589
 ↑ G. L. WerneckC. H. N. CostaA. M. WalkerJ. R. DavidM. WandJ. H. MaguireMultilevel modelling of the incidence of visceral leishmaniasis in Teresina, BrazilEpidemiology & Infection2007135219520110.1017/s09502688060068812s2.03384724794216824254
 ↑ R. A. ThompsonJ. W. De Oliveira LimaJ. H. MaguireD. H. BraudD. T. SchollClimatic and demographic determinants of American visceral leishmaniasis in northeastern Brazil using remote sensing technology for environmental categorization of rain and region influences on leishmaniasisThe American Journal of Tropical Medicine and Hygiene20026766486552s2.0003691257610.4269/ajtmh.2002.67.64812518857
 ↑ K. M. McNysetUse of ecological niche modelling to predict distributions of freshwater fish species in KansasEcology of Freshwater Fish20051432432552s2.02814446395810.1111/j.16000633.2005.00101.x
 ↑ P. NietoJ. B. MaloneM. E. BaviaEcological niche modeling for visceral leishmaniasis in the state of Bahia, Brazil, using genetic algorithm for ruleset prediction and growing degree daywater budget analysisGeospatial Health20061111512610.4081/gh.2006.2862s2.05444909513918686237
 ↑ M. G. ColaciccoMayhughP. M. MasuokaJ. P. GriecoEcological niche model of Phlebotomus alexandri and P. papatasi (Diptera: Psychodidae) in the Middle eastInternational Journal of Health Geographics20109, article 210.1186/1476072x922s2.076949105829
 ↑ C. GonzálezO. WangS. E. StrutzC. GonzálezSalazarV. SánchezCorderoS. SarkarClimate change and risk of leishmaniasis in North America: predictions from ecological niche models of vector and reservoir speciesPLOS Neglected Tropical Diseases201041, article e58510.1371/journal.pntd.00005852s2.076249128178
 ↑ D.E. A. ElnaiemJ. SchorscherA. BendallRisk mapping of visceral leishmaniasis: the role of local variation in rainfall and altitude on the presence and incidence of kalaazar in Eastern SudanThe American Journal of Tropical Medicine and Hygiene200368110172s2.00037263934
 ↑ M. A. OshaghiN. M. RavasanE. JavadianApplication of predictive degree day model for field development of sandfly vectors of visceral leishmaniasis in northwest of IranJournal of Vector Borne Diseases20094642472542s2.07344911671319959849
 ↑ L. S. PontryaginOn the zeros of some elementary transcendental functionsAmerican Mathematical Society Translations1955195110MR0073686
 ↑ L. S. PontryaginThe mathematical theory of optimal processes and differential gamesTrudy Matematicheskogo Instituta imeni VA Steklova1985169119158MR836572