Transmission network dynamics of Plasmodium Vivax Malaria P. Pongsumpun, and I.M.Tang Abstract One of the top ten killer diseases in the world is Malaria. In each year, there are between 3 to 5 million clinical episodes of malaria and 1.5 to 2.7 million deaths worldwide. The malaria disease is caused by the multiplication of protozoa parasite of the genus Plasmodium. Malaria in humans is due to 4 types of malaria parasites such that Plasmodium falciparum, Plasmodium vivax, Plasmodium malariae and Plasmodium ovale. P.vivax malaria differs from P. falciparum malaria in that a person suffering from P. vivax malaria can experience relapses of the disease. Between the relapses, the malaria parasite will remain dormant in the liver of the patient, leading to the patient being classified as being in the dormant class. In this paper, the dynamical model of P. vivax malaria is formulated to see the network distribution of this disease. Keywords Dynamical model, household, Plasmodium Vivax Malaria, relapse. M I. INTRODUCTION ALARIA is transmitted by the female anopheles mosquito. The Plasmodium genus of protozoa parasites has a life cycle which is split between a vertebrate host and an insect vector. The Plasmodium species, with the exception of Plasmodium malariae (which may affect the higher primates) are exclusively parasites of man. The sporozoites from the mosquito salivary gland are injected into the human as the mosquito must infect anticoagulant saliva to ensure an even flowing meal. Once in the human bloodstream, the sporozoites arrive in the liver and penetrate hepatocytes, where they remain for 9-16 days, multiplying within the cells. On release, they return to the blood and penetrate red blood cells in which they produce either merozoites or micro and macrogametocytes, which have no further activity within the human host. Another mosquito arriving to feed on the blood may suck up these gametocytes into its gut, where exflagellation of microgametocytes occurs, and the macrogametocytes are fertilized. The resulting ookinete penetrates the wall of a cell in the midgut, where it develops into an oocyst. Sporogeny within the oocyst produces many Manuscript received October 24, 28: Revised version received December 1, 28. P. Pongsumpun is with the Department of Mathematics and Computer Science, Faculty of Science, King Mongkut s Institute of Technology Ladkrabang, Chalongkrung road, Ladkrabang, Bangkok 152, Thailand. (corresponding author phone: 662-737-3 ext. 6196; fax: 662-326-4344 ext.284; e-mail: kppuntan@kmitl.ac.th). I.M. Tang is with the Department of Physics, Faculty of Science, Mahidol University, Rama 6 road, Bangkok 14, Thailand. sporozoites and, when the oocyst ruptures, the sporozoites migrate to the salivary gland, for injection into another host [1]. Not only are human beings the host (vertebrate) of human Plasmodium but also the Anopheles mosquitoes are also a host (invertebrate). The life cycle of Plasmodium vivax is complex and differs from those of the other malarias. Plasmodium vivax malaria is induced into bloodstream by the bite of an infected mosquito in the sporozoite form. The sporozoites travel to the liver and cells divide by asexual reproduction into what is known as merozoite within 7-8 days [2]. The merozoite can then invade the red blood cell and cause the illness. Some sporozoite will however develop into hypnozoite (there is no equivalent state in the other malarias). These hypnozoites will lay dormant in the liver for a long time. When the patient is exhausted, frail or more susceptible some of them will transform themselves into merozoites and produce the illness again. These relapses can occur up to several months or years (average 5 years) after the primary infection. The relapses peter out when there are no more hypnozoite in the liver [2]. The infection can also be transmitted accidentally through blood transfusion when the donating individual has the malaria parasite. This is one of the reasons why people who have been infected with disease can never donate blood. Congenital infection of a newborn from an infected mother also happens, but it is comparatively rare [2]. Natural transmission of malaria depends on the presence of the relationship between the three basic epidemiological factors: the host, the agent, and the environment. The most important environmental factors are temperature and humidity. Malaria parasites stop developing in the mosquitoes when the temperature is below 16 C. The best conditions for the development of diseases and the transmission of the infection are when the average temperature is within the range of 2 3 C, and the average relative humidity is at least 6% [3]. A high relative humidity lengthens the life of the mosquito and it enables them to live long enough to transmit the infection to several persons. Data of Malaria cases due to the different types of malaria between during 23 to 26 in Thailand indicates that the rate of increase of the incidence of P. vivax is the highest among the four types. This can be seen in figure 1. Issue 3, Volume 3, 29 275
Averaged increasing rate per 1, populations 14 12 1 8 6 4 2 Plasmodium Falciparum Plasmodium Plasmodium Vivax Mixed Type Malariae Type of Malaria suggested [2] that the human population should instead be divided into three states; non-infected, infected but without any acute clinical signs, infected with acute clinical sign, to better reflect the clinical status of the individual. Others believe that the population should be divided into susceptible, infected but not infectious and infected and infectious. In this study, we construct the dynamical model of P.Vivax malaria between villages and houses to see the network distribution of this disease. Fig.1. The average increasing rate per 1, populations of Malaria in Thailand during 23-26[4]. The transmission of malaria is usually described by the Ross-MacDonald (RM) model [2]. This model is only suitable for the transmission of the P. falciparum malaria since it does not contain a role of possible relapses of the illness. In the present state of concern for medical safety, there is no place for human experimentation to see what would happen if new therapies were adopted. Mathematical modeling allows one to simulate what would occur. In this study, the dynamics transmission of P. vivax between houses and villages are considered. We introduce in Section 2, the mathematical model which would make it applicable to the transmission of P. vivax malaria in the villages and district levels. In Section 3, we simulate the consequences of changing the number of villages, number of houses in each village, contact rates and relapse rates. We discuss in Section 4, the implication of the insights obtained from the simulations. Part of the urgency for doing research on P. vivax malaria is due to the fact that P. vivax malaria is becoming an emerging public health problem. It is estimated that about 5% of the malaria cases outside of Africa and 1% in Africa are due to P. vivax and that the percentages are soaring. Fig. 2. The diagram shows the movement of the population between villages. II. TRANSMISSION MODEL In 1911, Ross [5] formulated the mathematical model of the epidemiology of malaria (P. falciparum) and improved on by MacDonald [3]. In the Ross model, an individual in the human population is separated as being in a non-infected or infected state. This gives rise to what is known as a SIS (susceptible-infected-susceptible) model. It has been Issue 3, Volume 3, 29 276
Fig. 3. The diagram shows the movement of the population between houses. Suppose that there are q villages, L persons and u houses in each village. The persons move to any houses in this village by random process. This process is done by random the 1 st person to the L th person (with uniformly distribution) go to each house everyday. Each day, one person can go only one time in one house of this village. The probability for each person to visit each house is equal. There is no person come from the outside of this village. The persons who stay in any house at the beginning time will come back to their houses at the ending time. We assume that at the first day, there is only one infected human in one house and there is no infected human for the other houses in the village. In our model for the transmission of P. vivax in the local level (village), we divide the host (human) population into susceptible, infected, dormant and recovered classes. The last category, the recovered are susceptible to further infections and so they reenter into the susceptible class. The vector population is separated into susceptible and infected classes. The variables in our model are defined as follows: Sh t,i is the number of susceptible persons in house i th after visited at day t, Ih t,i is the number of infected persons in house i th after visited at day t, Dh t,i is the number of dormant persons in house i th after visited at day t, Rh t,i is the number of recovered persons in house i th after visited at day t, Sv t,i is the number of susceptible mosquitoes in house i th after visited at day t, Iv t,i is the number of infected mosquitoes in house i th after visited at day t, L is the total number of persons, q is the number of villages, u is the total number of house, g is the ending time. The dynamics of human and vector populations are given by ΔSh t,i = - γ Iv t,i Sh t,i + s 3 Dh t,i + (1-α) s 1 Ih t,i + s 4 Rh t,i ΔIh t,i = γ Iv t,i Sh t,i s 1 Ih t,i + s 2 Dh t,i - s 5 Ih t,i ΔDh t,i = α s 1 Ih t,i (s 2 + s 3 ) Dh t,i (1) ΔRh t,i = s 5 Ih t,i - s 4 Rh t,i ΔSv t,i = Y t,i β Sv t,i Ih t,i μ v Sv t,i ΔIv t,i = β Sv t,i Ih t,i - μ v Iv t,i Y t,i = μ v (Sv t,i +Iv t,i ) = Basic recruitment number of the mosquitoes in house i th at day t, β is the transmission rate of dengue virus from infectious person to susceptible vector, μ v is the death rate of the mosquitoes, α is the percentage of infected human in whom some hypnozoites remain dormant in the liver, (1-α) is the percentage of infected humans who recover and become susceptible again, s 1 is the rate at which a person leaves the infected class by recovering or by entering into the dormant class, s 2 is the rate at which the dormant human relapses back to the infected human, s 3 is the recovery rate of the dormant human, s 4 is the rate at which the recovered human relapses back to the susceptible human, and s 5 is the rate at which the infected human recovers, γ is the rate at which the P. vivax parasite is transmitted from the mosquito to the human and is given by [6-9] β h (2) γ = b N T + m where b is the specie-dependent biting rate of the mosquitoes; m is the population of other animals that the mosquitoes can feed on; N T is the total human populations and β h is the probability the parasite passed on by the mosquito will continue to thrive in the human. β is the rate at which the mosquitoes become infected with the Plasmodium vivax parasite once the mosquito has bitten an infected human. β is defined by [6-9] βv β = b N T + m where β v is the probability the parasite passed to the mosquito by biting human. III. NUMERICAL SOLUTIONS The simulations of the model are calculated to see the time distributions of this disease. The time distributions of infected and dormant humans for the different situations are shown in the following figures. Issue 3, Volume 3, 29 277
4 4 villages 4 4 villages 6 villages 6 villages 8 villages 8 villages 3 1 villages 3 1 villages 2 Dorm ant hum an 2 1 1 1 41 81 121 161 21 241 281 321 361 5 1 15 2 25 3 35 4 Fig.4 Model outputs display the time distribution of infected human for the different number of villages in each district. The parameters in our model are s 1 = 1/14, s 2 =1/(365*3), s 3 = 1/(25), s 4 =1/(365*1), s 5 =1/3, γ =.55, β =.45, μ v =.4, L = 1,, g = 365. Fig.6. Model outputs display the 2D D h -I h plane for the different number of villages in each district. The parameters are same as fig. 4. 4 3 4 villages 6 villages 8 villages 1 villages 4 3 3 houses 4 houses 5 houses 6 houses 2 2 1 1 1 41 81 121 161 21 241 281 321 361 1 41 81 121 161 21 241 281 321 361 Fig.5. Model outputs display the time distribution of dormant human for the different number of villages in each district. The parameters are same as fig. 4. Fig.7. Model outputs display the time distribution of infected human for the different number of house in each village. The parameters in our model are s 1 = 1/14, s 2 = 1/(365*3), s 3 = 1/(25), s 4 =1/(365*1), s 5 =1/3, γ =.55, β =.45, μ v =.4, L = 1, g = 365. Issue 3, Volume 3, 29 278
4 3 houses 4 houses 5 houses 6 houses 4 gamma =.15, beta =.1 gamma =.25, beta =.15 gamma =.45, beta =.35 3 3 gamma =.55, beta=.45 2 2 1 1 1 41 81 121 161 21 241 281 321 361 1 41 81 121 161 21 241 281 321 361 Fig.8. Model outputs display the time distribution of dormant human for the different number of house in each village. The parameters are same as fig. 7. Fig.1. Model outputs display the time distribution of infected human for the different contact rates in each village. There are 3 houses in each village. The parameters are same as fig. 7 except the contact rates. 4 3 3 houses 4 houses 5 houses 6 houses 4 3 gamma =.15, beta =.1 gamma =.25, beta =.15 gamma =.45, beta =.35 gamma =.55, beta=.45 2 2 1 1 1 2 3 4 1 41 81 121 161 21 241 281 321 361 Fig.9. Model outputs display the 2D D h -I h plane for the different number of houses in each village. The parameters are same as fig. 7. Fig.11. Model outputs display the time distribution of dormant human for the different contact rates in each village. The parameters are same as fig. 1. Issue 3, Volume 3, 29 279
4 gamma=.15,beta=.1 gamma=.25,beta=.15 gamma=.45,beta=.35 gamma=.55,beta=.45 4 s2=1/365 s2=1/(365*3) s2=1/(365*5) s2=1/(365*1) 3 3 2 2 1 1 5 1 15 2 25 3 35 4 1 41 81 121 161 21 241 281 321 361 Fig.12. Model outputs display the 2D D h -I h plane for the different number of houses in each village. The parameters are same as fig. 1. Fig.14. Model outputs display the time distribution of dormant human for the different rate at which the dormant human relapses back to the infected human in each village. The parameters are same as fig. 13. 4 s2=1/365 s2=1/(365*3) 4 s2=1/365 s2=1/(365*3) 3 s2=1/(365*5) s2=1/(365*1) 3 s2=1/(365*5) s2=1/(365*1) n 2 Dormant huma 2 1 1 1 41 81 121 161 21 241 281 321 361 1 2 3 4 Fig.13. Model outputs display the time distribution of infected human for the different rate at which the dormant human relapses back to the infected human in each village. There are 3 houses in each village. The other parameters are same as fig. 11 except s 2. Fig.15. Model outputs display the 2D D h -I h plane for the different rate at which the dormant human relapses back to the infected human in each village. The parameters are same as fig. 13. Issue 3, Volume 3, 29 28
4 3 s4=1/365 s4=1/(365*3) s4=1/(365*5) s4=1/(365*1) 4 3 s4=1/365 s4=1/(365*3) s4=1/(365*5) s4=1/(365*1) 2 2 1 1 5 1 15 2 25 3 35 4 1 41 81 121 161 21 241 281 321 361 Fig.16. Model outputs display the time distribution of infected human for the different rate at which the recovered human relapses back to the susceptible human in each village. There are 3 houses in each village. The other parameters are same as fig. 1 except s 4. Fig.18. Model outputs display the 2D D h -I h plane for the different rate at which the recovered human relapses back to the susceptible human in each village. The parameters are same as fig. 16. 4 s4=1/365 IV. CONCLUSION 3 2 1 1 41 81 121 161 21 241 281 321 361 s4=1/(365*3) s4=1/(365*5) s4=1/(365*1) Fig.17. Model outputs display the time distribution of dormant human for the different rate at which the recovered human relapses back to the susceptible human in each village. The other parameters are same as fig. 16. In this study, we have compared the results of the simulation when different values of several parameters are used. Figures 4, 5 and 6 show the time distributions of infected human, dormant human and the 2D Dormant-Infected human plane, respective, when the number of villages in each district is difference. Figures 7, 8 and 9 show the time distributions of infected human, dormant human and the 2D Dormant- plane, respective; when the number of houses in each village is difference. Figures 1, 11 and 12 show the time distributions of infected, dormant human and the 2D Dormant- plane, respective; when the contact rates in each village are difference. Figures 13, 14 and 15 show the time distributions of infected human, dormant human and the 2D Dormant- plane, respective, when the rate at which the dormant human relapses back to the infected human in each village is difference. Figures 16, 17 and 18 show the time distributions of infected, dormant human and the 2D Dormant- plane, respective, when the rate at which the recovered human relapses back to the susceptible human in each village is difference. We will see that the epidemic sizes are higher when the smaller number of villages, the smaller number of households and the higher contact rates. But when the rate at which the dormant human relapses back to the infected human and the rate at which the recovered human relapses back to the susceptible human are higher, the epidemic sizes are smaller. The important parameters will be the role of traveling of people between villages. Issue 3, Volume 3, 29 281
ACKNOWLEDGMENT This work is supported by Commission on Higher Education and the Thailand Research Fund according to contract number MRG5878. The authors would like to thank Dr.Sopon Iamsirithaworn, Ministry of Public Health, Thailand for supporting the data. REFERENCES [1] J.P.Kreier, Malaria volume 1 epidemiology, chemotherapy, morphology, and metabolism. New York : Academic press; 198. [2] R.M.Anderson, RM.May. Infectious diseases of humans : dynamics and control. Oxford: Oxford University press;1992. [3] G.MacDonald. The epidemiology and control of malaria. London : Oxford University press ;. 1957 [4] Annual Epidemiological Surveillance Report, 23-26, Division of Epidemiology, Ministry of Public Health, Royal Thai Government. [5] R Ross, The preventation of Malaria, 2 th ed. Murray, London, 1911. [6] L.Esteva and C. Vargas, Analysis of a dengue disease trasmission model, Mathematical Bioscience, vol.15, pp.131-151, 1998. [7] P.Pongsumpun, and I. M. Tang, Mathematical model for the transmission of Plasmodium Vivax Malaria, International Journal of mathematical models and methods in applied sciences, vol. 3, pp. 117-121, 27. [8] P.Pongsumpun, and I. M. Tang, Limit Cycle and Chaotic Behaviors for the Transmission Model of Plasmodium Vivax Malaria, International Journal of mathematical models and methods in applied sciences, vol.2, pp.563-57, 28. [9] P. Pongsumpun, and I. M. Tang, The Transmission Model of P.falciparum and P.Vivax Malaria between Thai and Burmese, International Journal of mathematical models and methods in applied sciences, 28, accepted. Issue 3, Volume 3, 29 282