A mathematical model for assessing the effectiveness of controlling relapse in Plasmodium vivax malaria endemic in the Republic of Korea
Authors:
Sungchan Kim aff001; Jong Hyuk Byun aff001; Anna Park aff001; Il Hyo Jung aff001
Authors place of work:
Department of Mathematics, Pusan National University, Geumjeong-Gu, Busan 46241, Republic of Korea
aff001; Finance · Fishery · Manufacture Industrial Mathematics Center on Big Data, Pusan National University, Geumjeong-Gu, Busan 46241, Republic of Korea
aff002
Published in the journal:
PLoS ONE 15(1)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0227919
Summary
Malaria has persisted as an endemic near the Demilitarized Zone in the Republic of Korea since the re-emergence of Plasmodium vivax malaria in 1993. The number of patients affected by malaria has increased recently despite many controls tools, one of the reasons behind which is the relapse of malaria via liver hypnozoites. Tafenoquine, a new drug approved by the United States Food and Drug Administration in 2018, is expected to reduce the rate of relapse of malaria hypnozoites and thereby decrease the prevalence of malaria among the population. In this work, we have developed a new transmission model for Plasmodium vivax that takes into account a more realistic intrinsic distribution from existing literature to quantify the current values of relapse parameters and to evaluate the effectiveness of the anti-relapse therapy. The model is especially suitable for estimating parameters near the Demilitarized Zone in Korea, in which the disease follows a distinguishable seasonality. Results were shown that radical cure could significantly reduce the prevalence level of malaria. However, eradication would still take a long time (over 10 years) even if the high-level treatment were to persist. In addition, considering that the vector’s behavior is manipulated by the malaria parasite, relapse repression through vector control at the current level may result in a negative effect in containing the disease. We conclude that the use of effective drugs should be considered together with the increased level of the vector control to reduce malaria prevalence.
Keywords:
Infectious disease control – Parasitic diseases – Malaria – Plasmodium – Malarial parasites – Korea – Mosquitoes – South Korea
Introduction
Malaria has persisted in the Republic of Korea since the 1993 Plasmodium vivax (P. vivax) malaria re-emergence [1]. According to KCDC (Korea Centers for Disease Control & Prevention), the Korean Government Health Authorities’ continued efforts had managed to reduce the number of malaria patients to a few hundreds, but recently the number has started to increase again.
Globally, P. vivax malaria was one of the most neglected diseases in the world due to its low rate of fatality, but recently there has been much interest ever since it has been found that severe symptoms and deaths due to vivax malaria are not uncommon [2], even in Korea [3]. In particular, P. vivax is known to produce hypnozoites, which are dormant forms of the parasite residing in the liver of the victim [4]. Activation of these hypnozoites makes future relapse possible even if blood-stage malaria is removed from primary infection. Hypnozoites are produced when the conditions for transmission of the parasite are less favorable, for instance, during short seasons of appropriate temperatures and vectors available at high latitudes [5]. The relapse rate caused by these hypnozoites is around 68% across the world, without any popular treatment for such relapses [6]. Therefore, treatment to prevent malaria relapse by removing the remaining hypnozoites in the liver is imperative alongside blood-stage malaria treatment.
Treatment of P. vivax malaria in Korea generally uses Chloroquine for schizonticide in the blood and Primaquine for detoxification in the liver, termed ‘radical cure’ [7]. Until now, Primaquine has been the only drug used for radical cure. Primaquine has a short half-life, so patients faced considerable difficulty in taking their daily medications during the recommended period of 14 days during which dormant hypnozoites in the liver were treated [8]. Thus, the misuse of drugs in patients treated with blood-stage malaria often resulted in poor treatment rate for hypnozoites [8]. However, with Tafenoquine being approved as a new drug in relapse treatment by the U.S. FDA in 2018. Tafenoquine only requires a single dose to prevent relapse of malaria due to its long half-life of about 16 days [9]. The development of this therapy reduces the effort required from an individual during his treatment. Therefore, proper use of Tafenoquine in addition to the use of Primaquine may lower the relapse rate even at the current level [10].
Meanwhile, a model for evaluating malaria relapse patterns has been introduced by Roy et al. [11]. In this model, the authors developed a mathematical model that uses gamma-distributed lag-time to cover the mosquito data and the absence of relapse dynamics, demonstrating a new technique to estimate the rate of relapse of infections in northwest India. In addition, based on the estimated values of the parameters, the effect of anti-relapse on the area was evaluated. However, the model presented above is not suitable for use in our research because, unlike northwest India, the northernmost part of Korea has a yearly variation in temperature of more than 30°C and is not a mosquito-friendly environment between November and February. Additionally, due to the higher latitudes of the region, the distribution of the period between the initial infection and the first relapse has to be expressed in a multimodal form. Furthermore, the distribution of latent period of malaria in humans is expressed in short-term and long-term bimodal forms. The Gamma distributed model—a generalized exponential model—has been widely used for modeling realistic distributions that the exponentially distributed model cannot be used for [12–16]. However, it was revealed that the exponential or gamma distributions, which are commonly used to model the incubation and relapse periods, are inadequate in capturing the complexity of the event-time distribution of P. vivax malaria infections in the temperate regions being considered [5]. Therefore, for this study, a transmission model that considers multimodal distributions is essential. In this paper, we develop a transmission model to evaluate the rate of malaria relapse infections in the northern part of Korea and to examine its effect at the population-level on radical cure. This model takes into account data regarding the mosquito population of the region and the intrinsic distribution of relapses based on existing literature on malaria models in Korea [17–19]. To model intrinsic distributions of incubation periods and the time between the initial infection and the first relapse, we introduce a Coxian distribution, which can express multimodal phenomena. Our model measures the effectiveness of the new mass treatment for relapse at the population-level.
Materials and methods
Data
Incidence in targeted area
In this study, we focus on the outbreak of P. vivax malaria in Yeoncheon-gun, Gyeonggi-do in Korea in 2007. Yeoncheon-gun is one of the places in Korea that suffers from the most severe malarial endemic. Geographically, it is located near the Demilitarized Zone (DMZ). Anopheles sinensis is the most popular malaria-bearing mosquito in the region. We used weekly malaria incidence data from the Infectious Disease Portal of Korea Centers for Disease Control & Prevention (KCDC). This region has four distinct seasons with average temperatures in summer and in winter. According to Korea Meteorological Administration (KMA), the lowest average monthly temperature is -7.0°C in January, and the highest average monthly temperature is 23°C in August. However, the rainfall is heavy enough to reach 1,300 mm, which is concentrated mostly in summer. Therefore, mosquitoes thrive during the summer, but their population is significantly reduced in winter. Thus, malaria in humans is most prevalent during summer, especially in July and August, and conversely, is least prevalent in winter, especially in December and February.
Incubation periods and the time to first relapse
First, a study containing incubation period (IP) data for a total of 229 individuals was identified [20]. This study provides analysis of the data of the incubation periods from cases of civilians infected with P. vivax malaria in Korea using data from the malaria epidemiological survey conducted between 2001 to 2010. As malaria in Korea is mainly concentrated in the northwestern region, near Yeoncheon-gun, it is reasonable to use malaria incubation data of the entire country for this study. Data points were recorded from the histogram presented in the study. Next, the time to first relapse (TTFR) data was obtained from [21]. This paper presents raw TTFR data for 46 relapsed samples from South and North Korea.
Fig 1 depicts the histograms indicating the relative frequencies of IPs and TTFRs for P. vivax malaria in human. The distribution of the histogram clearly explains why the malaria parasite survives, although mosquitoes cannot survive in the region during winter. Both incubation and relapse occur in two forms—short-term and long-term. Malaria parasites undergo short-term incubations of about 1-2 weeks and long-term incubations of about a year [22], and relapse due to dormant hypnozoites in the liver occurs about 8-10 weeks after the primary attack in the short-term and about a year afterwards in the long-term, depending on the specific strain of virus [23, 24]. The spread of the disease during the succeeding year is also affected by this.
Mathematical model
We formulate a mathematical model to capture the trend of relapse of P. vivax malaria at the population-level, which has been schematically depicted in Fig 2:
In the model (1), the component related to humans is formulated by dividing them into five classes (Table 1): S, E, Ip, L and Ir, which describe the groups of susceptible, exposed, primary infected, humans who carry liver-stage dormant hypnozoites, and infectious individuals who have relapsed from the re-emergence of hypnozoites. The population M of female Anopheles mosquitoes is divided into three classes: SM, standing for susceptible vectors, EM, standing for exposed vectors, and IM, standing for infectious vectors. The total rate in each case remains constant, i.e., S(t) + E(t) + Ip(t) + L(t) + Ir(t) = 1 and SM(t) + EM(t) + IM(t) = 1. The ratio of the number of mosquitoes to that of humans, ρ, represents the seasonality of transmission, and therefore ρ ≡ ρ(t). As demonstrated in the previous section, one of the most important considerations in modeling P. vivax malaria in temperate regions is expressing the multimodal nature of IP and TTFR. Therefore, we consider P(t) and Q(t), which are the survival functions of the incubation period and the duration of the overall liver stage, TTFR, respectively, to reflect realistic distributions.
When a malaria-infected mosquito bites a susceptible human, he/she becomes infected at a rate of λ(t) and undergoes an incubation period, which follows the survival probability, P(t) in the body. Following that, he/she becomes infectious and, with subsequent treatment, leaves the Ip-class at the rate of γ. One section of this population enters the L-class at the rate of ξ, as the treatment of hypnozoites in the liver is still not complete, and the other section enter the S-class again at the rate of 1 − ξ after successful completion of treatment. The period that each individual spends in the L-class until the re-emergence of malaria follows the survival probability Q(t). The relapsed malaria patient in the Ir-class then begins to infect susceptible mosquitoes with a force of infection of λM(t). Relapses can occur for three main reasons: re-infection from a mosquito, re-emergence from untreated hypnozoites in the liver, and from untreated blood-stage malaria. However, it is known that blood-stage malaria has a short treatment period and a high rate of treatment. Therefore, the last reason can be deemed to be statistically insignificant in our case, and we only consider re-infections caused by mosquitoes and untreated hypnozoites. Detailed descriptions of relevant parameters are recorded in Table 2.
We consider the control parameter κ to observe the effect of group relapse prevention control, which prevents relapse due to hypnozoites. We assume that control at a rate κ reduces Ip-to-L transmission rate ξ by (1 − κ)ξ. Therefore, control at a rate 100 × κ% enables individuals from the Ip-class to enter L-class at a rate (1 − κ)ξγ, and to enter the S-class at a rate (1 − (1 − κ)ξ)γ. Note that κ = 0 represents the current level of radical cure, and κ = 1 represents the successful control of all primary infected patients occurs.
In addition, we introduce our chosen functions P and Q as the survival functions of the Coxian distribution. The Coxian distribution, one of the phase-type distributions, is chosen because of its property that any type of non-negative distributions can be approximated to a Coxian distribution [25]. The detailed derivation of the Coxian distributed model is as follows: Let X and Y be random variables representing the dwelling-times in the exposed stage(E) and stage of carrying dormant hypnozoites in the liver (L) in humans, respectively. Then the probability density functions of X and Y, fX and fY, respectively, are defined to be
We successfully derived the Ordinary Differential Equation model from (1) in the case in which P and Q are assumed to correspond to Coxian distributions (4). The techniques for derivation are as follows: we put the 1 × n vector φ as
In many cases, the threshold principle of in a constant environment is dealt with by using Diekmann-Heesterbeek-Metz(DHM) definition [27] of the basic reproduction number, R0. However, it is not applicable to the formulation of threshold principles for population growth in a periodic environment [28]. In particular, the definition of R0 in a periodic environment was given by Bacaër and Guernaoui [29], which can be interpreted as the asymptotic per generation. Motivated by the work, H. Inaba presented a new definition for the basic reproduction number as the spectral radius of the generation evolution operator (GEO) [28]. Moreover, it was shown that the supercritical condition R0>1 implies existence of positive periodic solution [30]. These facts together yield the existence of a periodic solution to our model that is positive (1′) at which the reproduction number, according to the GEO definition, is bigger than unity.
Estimating model parameters
As the Coxian chains in the model (1′) cannot be biologically interpreted because of phenomena that are going out of the absorption state without going through the all the chains [31], we estimate model parameters using two procedures as follows:
- (A)
Fit survival functions P and Q to the IP and TTFR data, respectively.
- (B)
Fit other model parameters pH, pM and ξ in the model (1′) to the parameter values related to P and Q which are obtained in (A).
More detailed descriptions are as follows: First, we construct the IP and TTFR data as depicted in Fig 1 as an empirical cumulative distribution function (ECDF) and fit these to Coxian CDF (1 − P) and (1 − Q) in (4). That is, we want to find P and Q to minimize
To estimating pH, pM and ξ, we performed a fit of averaged weekly cases for each of the 4 weeks to the incidence data for 4 weeks using our model. Using lsqnonlin in Matlab, we found the values of the parameters that minimize the following objective functional:
The seasonal reproduction number
The seasonal (or effective) reproduction number Rs—an alternative form of the DHM definition of the basic reproduction number—is used when when parameters in the model are seasonally periodic [34]. The method of calculation of Rs is as follows: Rs is determined by the spectral radius of the next-generation operator. Moreover, the next generation operator is given by the following 2 × 2 block matrix in the case of a single vector-borne diseases:
Motivated by [31], we derive the Rs as follows:
Motivated by [35], we consider the proportion of the contribution of relapse in the seasonal reproduction number, Rs. Using (7), we can reformulate Rs as the square root of the sum of the primary-based part, Rp, and the relapse-based part, Rr:
Results
Parameters estimation
Parameters consisting of P and Q and fitted to ECDFs of IP and TTFR to minimize (5), are given in Table 3. Fitted curves agree well with the data, as shown in Fig 3, and satisfactorily express multimodal phenomena well.
Next, the model parameters are set as follows: we set the force of mortality of a human as μ, with the assumption that the death age of of death an individual is distributed exponentially, so that 1/μ is the average life expectancy. We followed concepts of mosquito physiological parameters, such as μM and ρ(t), as discussed in [17]. These parameters were fitted to the observed 2007 mosquito population data in Paju-si and Cheolwon-gun, near Yeoncheon-gun, with exponential and Gaussian assumptions. Since the period of recovery for humans is about 2 weeks [17], we assume 1/γ = 2 = 1/γr. We also set 1/ϵM = 7/9 because the mean incubation period for mosquitoes is 9 days [36]. The biting rate of a mosquito with respect to a human, b, is assumed to be the average number of bites for that mosquito per week multiplied by the associated Human Blood Index (HBI), which is defined to be by the proportion of the blood in a mosquito population obtained from human [37]. We set b = 0.138 with the HBI is assumed to be in the range [0.01, 0.1] and the gonotrophic cycle of mosquito is assumed as about 2.5 days following [17, 38]. The optimal values of the parameters that minimize (6) are pH = 0.643, pM = 0.492, ξ = 0.133 and Ur = 2.21 and the curve fits with the observed data are as shown in Fig 4.
Detailed values of parameters are listed in Table 4.
Effectiveness of relapse control
Contribution of relapse in the seasonal reproduction number
In (9), R% does not depend on seasonality, and is constant. The left panel in Fig 5 depicts the values of the seasonal reproduction number with changing time. The maximum value occurs in summer between July and August, when its value is 2.58 in the absence of relapse control, and 2.42 in the presence of relapse control with complete coverage, respectively. The right panel shows the percentage of the contribution of relapse in the seasonal reproduction number, Rs, in the presence of a relapse treatment with coverage rate κ. It is shown that relapse treatment reduces the contribution proportion of relapse. With the current coverage of relapse while Primaquine is being used, the model shows a maximum relapse contribution proportion of 11.7% in Korea. However, if the coverage of the relapse treatment were 90%, it is possible to reduce R% to 1.3%.
Malaria prevalence
Fig 6 illustrates the results of the malaria prevalence model using the parameters estimated in the previous section. During a year, 11.8% of cases resulted from relapse in total prevalence. Overall, primary infection is more prevalent than relapsed infection, but in winter, infected humans who have relapsed are more common. While total prevalence of malaria in humans tends to be similar to the prevalence of primary infections, relapsed prevalence remains more or less constant throughout the year, besides having a peak during summer, just like the case of primary infections. This is because relapsed infections are less affected by changes in the mosquito population than primary infections.
Fig 7 depicts the effect of constant control over 5 years. The gradual weakening of the scale of outbreaks over several consecutive years suggests the possibility of complete elimination of the parasites. 90% coverage of relapse also reduces relapsed prevalence to zero within 5 years. A further finding is that control not only reduces the relapse prevalence but also reduces the rate of primary infections. The upper right panel of the figure shows that the peak of prevalence of primary infections reduces by more than 2 times after 5 years of relapse treatment with 90% coverage.
Eradication of malaria
Fig 8 summarizes the results related to the percentages of reduction in yearly cases after application of control that could explain the eradication of malaria. Theoretically, control at an application rate of 20% for a very long time can reduce cases by more than 90%, and can exterminate malaria at an application rate of approximately 40%. However, in the case of low coverage control over a period of 5 years or 10 years, control can reduce a great proportion of malaria incidences, but extinction is not possible. This means that malaria extinction is impossible via short-term intensive control.
Potential effect of manipulated behavior of mosquito caused by malaria parasites
Recent studies have shown that manipulated behaviors caused by malaria parasites in mosquitoes can increase the force of infection of malaria by more than twofold [40, 41]. These behavioral changes are caused by a survival instinct in malaria parasites, and therefore they tend to intensify as malaria is controlled [42]. Hence, the absence of these manipulation phenomena in the usual scenario can be an obstacle to identify appropriate intervention schedules for vector-borne disease control efforts [43].
We describe the potential effect of this feeding behavioral change over long periods of control in Fig 9. The solid line in the left panel shows that the bifurcation between the percentages of reduction in yearly cases with a 5-year control is positive and negative. Although a high rate of control has been implemented, it shows that a slight increase in the force of infection caused by control methods can have a negative manipulative effect without being significantly affected by the delay at the start of the manipulation.
Discussion
We used a transmission model to characterize the relapse of malaria considering realistic intrinsic distributions of incubation periods and the time to first relapse. We used Coxian distributed IP and TTFR distribution to model bimodality of incubation and relapse, which are distinct from other countries. Thus, the performance of our suggested method for statistical inferring of relapse related parameters from prevalence data is lower when considering with other countries. To the best of our knowledge, there has been no prior modeling study considering heuristic time to first relapse data for malaria dynamics in Korea. Moreover, incubation time and realistic relapse time distributions are considered more than others in the existing literature (refer to the S1 Appendix).
The seasonal reproduction number was used to calculate the proportion of contribution of relapse, indicating the effect of relapse on the epidemiology of malaria in the Republic of Korea. By dividing the number by two parts, primary infection based part and relapse based part, we could get the proportion explicitly. Next, we performed the study to demonstrate the effectiveness of anti-relapse treatment on the assumption that relapse treatments with greater coverage compared to the present could be used. The treatment was assumed to have a higher coverage rate than in the present, and we conducted the study under the assumption that the treatment rate was κ at the current rate. The results show that the control of relapse reduces total prevalence of malaria by reducing the numbers of both relapsed and primary infected individuals even though relapsed cases are significantly fewer in number than primary infected cases. We also studied whether malaria eradication is possible through relapse treatment. Theoretically, eradication can be achieved with a relapse coverage rate of 40% after a very long time, but eradication of malaria within 10 years is not possible even if the relapse is controlled at a very high rate. Therefore, further developments in vector control or prevention strategies are necessary. Moreover, considering the manipulated behavior of mosquitoes, our results demonstrate that negative effects on control might be yielded even if a high rate of treatment is implemented. We did not consider drug resistance in this study. If this were considered, a strategy for eradication of malaria in the near future might be necessary, which is difficult to control by merely controlling relapse. Therefore, it is proposed that mass treatment using new drugs should be performed with gradual increase in vector control considering the mosquitoes’ feeding behavioral changes. The study will also help disease prevention authorities in implementing timely and effective malaria control measures when new malaria drugs become commercially available.
Our work has several limitations. We used a data fitting method to quantify the current level of relapse treatment with absence of time series relapsed data and time to first relapse data for each patient. In addition, population of malaria mosquitoes in the Republic of Korea have a bimodal form that rapidly increases around May, temporarily decreases during the rainy season and then increases rapidly after the rainy season ends. As such, the bimodality of the number of patients shown in Fig 4 is affected by the number of mosquitoes. However, in this model, the influence of the peak around May is designed only by relapse and incubation. This reason may cause overestimation of total relapse prevalence in this model although about 5% of estimated cases of official relapse in 2018 were reported in the Republic of Korea [7]. If further data are investigated, suggested method may be applied to the current management of malaria in Korea by investigating the specific contents that can be applied centering on policy application.
Conclusion
In this work, we constructed a mathematical model to quantify the relapse of malaria in Korea and to evaluate the effectiveness of control methods in reducing relapse. Our model was able to express the distribution of multimodal incubation periods and the time to first relapse in continental climate zones where seasonality is evident. In this study, it was found that control of malaria relapse in Korea and the progression of control under the current modes of vector control can reduce malaria in the short-term but can have potentially negative effects in the long-term. The results suggest that control at higher vector control levels may result in reduction of malaria.
A number of recommendations for further research may be provided. For example, considering the monthly differences in the intrinsic distribution, the model generated by the age structured PDE may be reconstructed to exhibit the temporal effect on malaria prevalence. In addition, combinations with higher vector control will allow the reconstruction of the model to determine the appropriate timings and strategies for malaria control.
Supporting information
S1 Appendix [pdf]
Comparison of distributions fit with ones in existing literatures.
Zdroje
1. Feighner BH, Pak SI, Novakoski WL, Kelsey LL, Strickman D. Reemergence of Plasmodium vivax malaria in the republic of Korea. Emerging Infectious Diseases. 1998;4(2):295. 9621202
2. Naing C, Whittaker MA, Wai VN, Mak JW. Is Plasmodium vivax malaria a severe malaria?: a systematic review and meta-analysis. PLoS Neglected Tropical Diseases. 2014;8(8):e3071. doi: 10.1371/journal.pntd.0003071 25121491
3. Im JH, Kwon HY, Baek J, Park SW, Durey A, Lee KH, et al. Severe Plasmodium vivax infection in Korea. Malaria Journal. 2017;16(1):51. doi: 10.1186/s12936-017-1684-4 28129766
4. Markus MB. Malaria: origin of the term “hypnozoite”. Journal of the History of Biology. 2011;44(4):781–786. doi: 10.1007/s10739-010-9239-3 20665090
5. Lover AA, Zhao X, Gao Z, Coker RJ, Cook AR. The distribution of incubation and relapse times in experimental human infections with the malaria parasite Plasmodium vivax. BMC Infectious Diseases. 2014;14(1):539. doi: 10.1186/1471-2334-14-539 25280926
6. Cappellini MD, Fiorelli G. Glucose-6-phosphate dehydrogenase deficiency. The Lancet. 2008;371(9606):64–74. doi: 10.1016/S0140-6736(08)60073-2
7. 2019 malaria control manual. Korea Centers for Disease Control and Prevention (KCDC); 2019. Available from: http://www.cdc.go.kr/CDC/cms/content/mobile/06/143706_view.html.
8. Yeom JS, Park YK. Treatment of Korean vivax malaria in Korea. Journal of the Korean Medical Association. 2007;50(1):88–92. doi: 10.5124/jkma.2007.50.1.88
9. Lacerda MV, Llanos-Cuentas A, Krudsood S, Lon C, Saunders DL, Mohammed R, et al. Single-dose tafenoquine to prevent relapse of Plasmodium vivax malaria. New England Journal of Medicine. 2019;380(3):215–228. doi: 10.1056/NEJMoa1710775 30650322
10. Llanos-Cuentas A, Lacerda MV, Hien TT, Vélez ID, Namaik-larp C, Chu CS, et al. Tafenoquine versus primaquine to prevent relapse of Plasmodium vivax malaria. New England Journal of Medicine. 2019;380(3):229–241. doi: 10.1056/NEJMoa1802537 30650326
11. Roy M, Bouma MJ, Ionides EL, Dhiman RC, Pascual M. The potential elimination of Plasmodium vivax malaria by relapse treatment: insights from a transmission model and surveillance data from NW India. PLoS Neglected Tropical Diseases. 2013;7(1):e1979. doi: 10.1371/journal.pntd.0001979 23326611
12. Lloyd AL. Realistic distributions of infectious periods in epidemic models: changing patterns of persistence and dynamics. Theoretical Population Biology. 2001;60(1):59–71. doi: 10.1006/tpbi.2001.1525 11589638
13. Wearing HJ, Rohani P, Keeling MJ. Appropriate models for the management of infectious diseases. PLoS Medicine. 2005;2(7):e174. doi: 10.1371/journal.pmed.0020174 16013892
14. Brooks-Pollock E, Cohen T, Murray M. The impact of realistic age structure in simple models of tuberculosis transmission. PLoS One. 2010;5(1):e8479. doi: 10.1371/journal.pone.0008479 20062531
15. Vergu E, Busson H, Ezanno P. Impact of the infection period distribution on the epidemic spread in a metapopulation model. PLoS One. 2010;5(2):e9371. doi: 10.1371/journal.pone.0009371 20195473
16. Krylova O, Earn DJ. Effects of the infectious period distribution on predicted transitions in childhood disease dynamics. Journal of The Royal Society Interface. 2013;10(84):20130098. doi: 10.1098/rsif.2013.0098
17. Nah K, Kim Y, Lee JM. The dilution effect of the domestic animal population on the transmission of P. vivax malaria. Journal of Theoretical Biology. 2010;266(2):299–306. doi: 10.1016/j.jtbi.2010.06.032 20619273
18. Nah K, Nakata Y, Röst G. Malaria dynamics with long incubation period in hosts. Computers & Mathematics with Applications. 2014;68(9):915–930. doi: 10.1016/j.camwa.2014.05.001
19. Kim JE, Choi Y, Lee CH. Effects of climate change on Plasmodium vivax malaria transmission dynamics: A mathematical modeling approach. Applied Mathematics and Computation. 2019;347:616–630. doi: 10.1016/j.amc.2018.11.001
20. Kim SJ, Kim SH, Jo SN, Gwack J, Youn SK, Jang JY. The long and short incubation periods of Plasmodium vivax malaria in Korea: the characteristics and relating factors. Infection & Chemotherapy. 2013;45(2):184–193. doi: 10.3947/ic.2013.45.2.184
21. Battle KE, Karhunen MS, Bhatt S, Gething PW, Howes RE, Golding N, et al. Geographical variation in Plasmodium vivax relapse. Malaria Journal. 2014;13(1):144. doi: 10.1186/1475-2875-13-144 24731298
22. Chu CS, White NJ. Management of relapsing Plasmodium vivax malaria. Expert Review of Anti-Infective Therapy. 2016;14(10):885–900. doi: 10.1080/14787210.2016.1220304 27530139
23. Yeom JS, Park JW. Status of vivax malaria after re-emergence in South Korea. Infection and Chemotherapy. 2008;40(4):191–198. doi: 10.3947/ic.2008.40.4.191
24. Warrell DA, Gilles HM. Essential malariology. CRC Press; 2017.
25. Koole G. A formula for tail probabilities of Cox distributions. Journal of Applied Probability. 2004;41(3):935–938. doi: 10.1239/jap/1091543436
26. Buchholz P, Kriege J, Felko I. Input modeling with phase-type distributions and Markov models: theory and applications. Springer; 2014.
27. Diekmann O, Heesterbeek JAP, Metz JA. On the definition and the computation of the basic reproduction ratio R 0 in models for infectious diseases in heterogeneous populations. Journal of Mathematical Biology. 1990;28(4):365–382. doi: 10.1007/bf00178324 2117040
28. Inaba H. On a new perspective of the basic reproduction number in heterogeneous environments. Journal of Mathematical Biology. 2012;65(2):309–348. doi: 10.1007/s00285-011-0463-z 21842424
29. Bacaër N, Guernaoui S. The epidemic threshold of vector-borne diseases with seasonality. Journal of Mathematical Biology. 2006;53(3):421–436. doi: 10.1007/s00285-006-0015-0 16823580
30. Inaba H. The basic reproduction number R0 in time-heterogeneous environments. Journal of Mathematical Biology. 2019; p. 1–34.
31. Kim S, Byun JH, Jung IH. Global stability of an SEIR epidemic model where empirical distribution of incubation period is approximated by Coxian distribution. Advances in Difference Equations. 2019;2019(1):469. doi: 10.1186/s13662-019-2405-9
32. Kakizoe Y, Nakaoka S, Beauchemin CA, Morita S, Mori H, Igarashi T, et al. A method to determine the duration of the eclipse phase for in vitro infection with a highly pathogenic SHIV strain. Scientific reports. 2015;5:10371. doi: 10.1038/srep10371 25996439
33. Akaike H. Information theory and an extension of the maximum likelihood principle. In: Selected papers of hirotugu akaike. Springer; 1998. p. 199–213.
34. Massad E, Coutinho FAB, Lopez LF, Da Silva DR. Modeling the impact of global warming on vector-borne infections. Physics of Life Reviews. 2011;8(2):169–199. doi: 10.1016/j.plrev.2011.01.001 21257353
35. Gao D, Lou Y, He D, Porco TC, Kuang Y, Chowell G, et al. Prevention and control of Zika as a mosquito-borne and sexually transmitted disease: a mathematical modeling analysis. Scientific Reports. 2016;6:28070. doi: 10.1038/srep28070 27312324
36. Ishikawa H, Ishii A, Nagai N, Ohmae H, Harada M, Suguri S, et al. A mathematical model for the transmission of Plasmodium vivax malaria. Parasitology International. 2003;52(1):81–93. doi: 10.1016/s1383-5769(02)00084-3 12543150
37. Garrett-Jones C. The human blood index of malaria vectors in relation to epidemiological assessment. Bulletin of the World Health Organization. 1964;30(2):241. 14153413
38. Feng X, Zhang S, Huang F, Zhang L, Feng J, Xia Z, et al. Biology, bionomics and molecular biology of Anopheles sinensis Wiedemann 1828 (Diptera: Culicidae), main malaria vector in China. Frontiers in Microbiology. 2017;8:1473. doi: 10.3389/fmicb.2017.01473 28848504
39. Agency CI. The world factbook. Available from: http://www.cia.gov/library/publications/the-world-factbook/.
40. Cator LJ, Lynch PA, Read AF, Thomas MB. Do malaria parasites manipulate mosquitoes? Trends in Parasitology. 2012;28(11):466–470. doi: 10.1016/j.pt.2012.08.004 23044288
41. Cator LJ, Lynch PA, Thomas MB, Read AF. Alterations in mosquito behaviour by malaria parasites: potential impact on force of infection. Malaria Journal. 2014;13(1):164. doi: 10.1186/1475-2875-13-164 24885783
42. Churcher TS, Trape JF, Cohuet A. Human-to-mosquito transmission efficiency increases as malaria is controlled. Nature Communications. 2015;6:6054. doi: 10.1038/ncomms7054 25597498
43. Spence Beaulieu MR. The role of parasite manipulation in vector-borne diseases. Evolution, medicine, and public health. 2019;2019(1):106–107. doi: 10.1093/emph/eoz019 31289686
Článek vyšel v časopise
PLOS One
2020 Číslo 1
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Severity of misophonia symptoms is associated with worse cognitive control when exposed to misophonia trigger sounds
- Chemical analysis of snus products from the United States and northern Europe
- Calcium dobesilate reduces VEGF signaling by interfering with heparan sulfate binding site and protects from vascular complications in diabetic mice
- Effect of Lactobacillus acidophilus D2/CSL (CECT 4529) supplementation in drinking water on chicken crop and caeca microbiome