Open Access

Optimally timing primaquine treatment to reduce Plasmodium falciparum transmission in low endemicity Thai-Myanmar border populations

  • Saranath Lawpoolsri1, 2Email author,
  • Eili Y Klein3,
  • Pratap Singhasivanon2,
  • Surapon Yimsamran2,
  • Nipon Thanyavanich2,
  • Wanchai Maneeboonyang2,
  • Laura L Hungerford1,
  • James H Maguire1 and
  • David L Smith4, 5
Malaria Journal20098:159

DOI: 10.1186/1475-2875-8-159

Received: 17 May 2009

Accepted: 15 July 2009

Published: 15 July 2009



Effective malaria control has successfully reduced the malaria burden in many countries, but to eliminate malaria, these countries will need to further improve their control efforts. Here, a malaria control programme was critically evaluated in a very low-endemicity Thai-Myanmar border population, where early detection and prompt treatment have substantially reduced, though not ended, Plasmodium falciparum transmission, in part due to carriage of late-maturing gametocytes that remain post-treatment. To counter this effect, the WHO recommends the use of a single oral dose of primaquine along with an effective blood schizonticide. However, while the effectiveness of primaquine as a gametocidal agent is widely documented, the mismatch between primaquine's short half-life, the long-delay for gametocyte maturation and the proper timing of primaquine administration have not been studied.


Mathematical models were constructed to simulate 8-year surveillance data, between 1999 and 2006, of seven villages along the Thai-Myanmar border. A simple model was developed to consider primaquine pharmacokinetics and pharmacodynamics, gametocyte carriage, and infectivity.


In these populations, transmission intensity is very low, so the P. falciparum parasite rate is strongly linked to imported malaria and to the fraction of cases not treated. Given a 3.6-day half-life of gametocyte, the estimated duration of infectiousness would be reduced by 10 days for every 10-fold reduction in initial gametocyte densities. Infectiousness from mature gametocytes would last two to four weeks and sustain some transmission, depending on the initial parasite densities, but the residual mature gametocytes could be eliminated by primaquine. Because of the short half-life of primaquine (approximately eight hours), it was immediately obvious that with early administration (within three days after an acute attack), primaquine would not be present when mature gametocytes emerged eight days after the appearance of asexual blood-stage parasites. A model of optimal timing suggests that primaquine follow-up approximately eight days after a clinical episode could further reduce the duration of infectiousness from two to four weeks down to a few days. The prospects of malaria elimination would be substantially improved by changing the timing of primaquine administration and combining this with effective detection and management of imported malaria cases. The value of using primaquine to reduce residual gametocyte densities and to reduce malaria transmission was considered in the context of a malaria transmission model; the added benefit of the primaquine follow-up treatment would be relatively large only if a high fraction of patients (>95%) are initially treated with schizonticidal agents.


Mathematical models have previously identified the long duration of P. falciparum asexual blood-stage infections as a critical point in maintaining malaria transmission, but infectiousness can persist for two to four weeks because of residual populations of mature gametocytes. Simulations from new models suggest that, in areas where a large fraction of malaria cases are treated, curing the asexual parasitaemia in a primary infection, and curing mature gametocyte infections with an eight-day follow-up treatment with primaquine have approximately the same proportional effects on reducing the infectious period. Changing the timing of primaquine administration would, in all likelihood, interrupt transmission in this area with very good health systems and with very low endemicity.


Plasmodium falciparum is endemic in 87 countries and approximately 2.4 billion persons are at risk [1, 2]. As populations in South Asia and South East Asia regions have grown, so have the number of people at risk [35]. Now, approximately one billion people live in areas at very low risk of P. falciparum [2, 4]. Recently, the WHO Global Malaria Programme has called for countries with low and moderate transmission areas to eliminate malaria transmission from their entire territory [6]. In many of these areas, malaria transmission is suppressed through a combination of insecticide-treated nets, indoor residual spraying, and prompt, effective treatment with anti-malarial drugs [1]. However, despite implementation of all the control measures, in some areas, such as the border region with Myanmar in Thailand, malaria transmission persists. This suggests a need to critically re-examine malaria epidemiology and the parasite life cycle to identify new control points.

One of the components of an elimination programme is to reduce transmission from malaria patients by making the infectious period as short as possible [6]. Prompt treatment of clinical cases is important, but in humans, mature gametocytes are the only infective stage of malaria parasite, and most drugs do not kill mature gametocytes [7]. In addition, in the P. falciparum life-cycle, unlike in other Plasmodium spp., the appearance of infective gametocytes is delayed with respect to the erythrocytic-schizogony cycle, resulting in a delayed appearance of mature P. falciparum gametocytes in the peripheral blood about 7–15 days after the initial acute attack [7, 8]. While, an untreated infection lasts about six months, on average, this duration can be cut short to a few days after the incubation period [7]. Thus, in areas where asexual parasitaemia is cut short by effective treatment with anti-malarial drugs and gametocyte production is, therefore, limited, lingering gametocytes can maintain malaria transmission [911]. This issue is particularly important in areas with good access to medical clinics.

Although artemisinin derivatives have been shown to reduce gametocyte carriage by eliminating asexual parasites and immature gametocytes, only the 8-aminoquinolines, such as primaquine, are lethal to the mature gametocytes [7, 12]. Primaquine, which has been shown to be effective against mature P. falciparum gametocytes is rapidly absorbed with peak plasma concentrations reached about two hours after administration [13, 14]. However, it has a drug elimination half-life of approximately eight hours, so it only remains active against parasites for, at most, a few days [12, 15]. The WHO recommends the use of a single oral dose of primaquine along with an effective blood schizonticide to reduce transmission, particularly in low endemic areas [12]. In areas where an early detection and treatment programme is highly effective, patients generally receive treatment one or two days after the acute attack, or approximately three to five days before gametocyte maturation [16]. While the effectiveness of primaquine as a gametocidal agent is widely documented, the mismatch between primaquine's short half-life, the long-delay for gametocyte maturation, and the proper timing of primaquine administration have not been studied.

Here, a mathematical model was constructed to describe P. falciparum transmission dynamics in seven hamlets along the Thai-Myanmar border, an area with low, seasonal transmission, and a highly effective health system. The transmission model was extended to investigate the timing of follow-up primaquine administration on the duration of gametocyte carriage and the implications for malaria transmission in the area. The models were developed with the goal of critically evaluating malaria transmission and tailoring control measures in areas with extremely good health systems, where elimination of malaria is feasible.


Study site

A malaria transmission model was developed to closely track the observed transmission intensity over an eight-year period of seven hamlets in the Tanaosri subdistrict, Suanphung district in Ratchaburi Province, a mountainous area along the Thai-Myanmar border. The study site covered an area of around 50 km2 with approximately 3,500 inhabitants in about 500 households. Free diagnosis and treatment of malaria has been provided for people in the area since 1997 by the Rajanagarindra Tropical Disease International Center (RTIC), operated by Mahidol University, Thailand, the only malaria clinic in this area. Individuals commonly receive treatment for febrile malaria within two days of an acute attack. After deployment of the early detection and treatment programme, the peak incidence of clinical malaria has gradually decreased from approximately four cases per 100 persons in 1999 to about two cases per 100 persons from 2003 through 2005. A single dose of mefloquine had been used as the first-line drug for treatment of uncomplicated falciparum malaria until the year 2005; the standard regimen was then changed to a two-day artesunate-mefloquine combination therapy, as recommended by the Thai government, because of the increased mefloquine resistance in the area [17]. A single dose of 30 mg primaquine is given to all P. falciparum positive patients on the last day of the treatment course. In Thailand, the haemolysis after primaquine administration among glucose-6-phosphate dehydrogenase (G6PD) deficiency patients is relatively mild [18, 19]. A test for G6PD deficiency is not required before a single dose primaquine administration.

Malaria transmission in this area is markedly seasonal with a peak of transmission during the rainy season (April to July), mainly due to the fluctuation of the mosquito population. However, malaria transmission does continue at a low level during the dry season. The P. falciparum parasite rate (Pf PR) in the area is very low, according to regular active surveillance surveys that have sampled more than half of the study population since 2003 (Table 1). Parasites were not detected in three surveys in 2005. Although some infections with low parasite densities might have gone undetected, the low Pf PR and yearlong continuous transmission suggest that asexual blood-stage infection is rare and thus gametocyte carriage in the absence of asexual parasites must play an important role in maintaining mosquito infection and the persistence of malaria transmission in the area.
Table 1

Plasmodium falciparum parasite rate (Pf PR) according to active surveillance surveys in the study area between 2003 and 2005.



Survey month

Number of blood samples

# P. falciparum positive


Year 2003


























Year 2004


























Year 2005





















Transmission model

A deterministic model of the infection dynamics of human and mosquito populations in the study area was developed. In the models, the human population was divided into five compartments: state variables tracked the proportion that were susceptible (S h ), liver-stage only (L h ), asymptomatic infection with gametocyte carriage (A h ), clinical episode (C h ), and gametocyte carriage only (G h ). In this low malaria transmission setting, malaria immunity and super-infection are rare, so an immune population was not included in the model. In the model for the mosquito population, seasonal mosquito population dynamics was considered, so the model tracked the population density of susceptible (M), infected (Y), and infectious (Z) mosquitoes (Figure 1). The human and mosquito infection dynamics were connected by human blood feeding of uninfected mosquitoes on infectious humans, and of infectious mosquitoes on uninfected humans. The dynamics of malaria infections in the two host populations are, thus, described by the following differential equations.
Figure 1

Schematic illustration of the malaria transmission model. The diagram shows the compartments of human population (Top) and mosquito population (Bottom) and transmission processes. For human population: Sh = susceptible; Lh = liver-stage only; Ah = asymptomatic infection with gametocyte carriage; Ch = clinical episode; and Gh = gametocyte carriage only. For mosquito population: M = susceptible; Y = infected; and Z = infectious.

Equations for the human population:
Equations for the mosquito population:

Susceptible humans become infected when exposed to infectious mosquitoes at the rate baZ, where a is the human feeding rate and b is the probability of transmission [20, 21]. A large proportion of infected patients (P) become symptomatic at the rate ν, an inverse of the latent period [22], these symptomatic individuals are subsequently detected and treated by passive surveillance at the clinic. Because of the delayed development of P. falciparum gametocytes, the rate that patients become infectious (ξ) after an acute attack depends on the duration of the gametocyte maturation process [7, 8]. Individuals normally remain infectious for a period of time even if the infection is treated and cured. Infections with gametocytes lose their infectivity to mosquitoes at the rate σ, which depends in part on the different treatment regimens [13, 23].

Asymptomatic asexual blood-stage infections are occasionally observed in low-endemic malaria areas [24]. Therefore, a small proportion of infected individuals (1-P) was assumed to become asymptomatically infected. In addition, asymptomatic cases imported from neighboring areas occur at the rate μ. These asymptomatic patients are likely to remain untreated and infectious for a long period, until the gametocytes are naturally eliminated at rate r [7]. The imported rate of asymptomatic cases was assumed to be equal to the exported rate of susceptible individuals, so that the population size remains constant.

For the mosquito population, adult mosquitoes emerge from larval habitat, which is modeled with a sinusoidal function for seasonal forcing, λ(1+sin(2πt/365). Mosquito infection occurs at rate ac, where c is the probability of transmission when they feed on infectious humans [21]. These infected mosquitoes subsequently become infectious at the rate q. A constant death rate, g, was applied to all mosquito classes.

All parameter estimates were obtained from published literature and unpublished data from the study area; except for the average mosquito birth rate, λ, which was estimated by fitting the model to the observed malaria occurrence data, and the waiting time to clear gametocytes under different treatment regimens, which was computed using another model (see below). Details and value estimates of parameters in the models are described in Table 2.
Table 2

Description of parameters and parameter estimated used in the transmission model.



Point estimates


Human feeding rate: Number of bites on human, per mosquito, per day, i.e., the product of the number of bites per mosquito per day and the proportion of bites on humans.

0.4 day-1


Parasite transmission probability of mosquitoes to humans: The probability that an infectious mosquito transmit the parasite to a human from a single bite.



Parasite transmission probability of humans to mosquitoes: The probability that a mosquito become infected from a single bite on an infectious human.



The rate that an infected human becomes positive for malaria parasite: An inverse of latent period.

1/18 day-1


Treatment coverage: The proportion that infected humans receive the malaria treatment



The rate that an symptomatic human becomes infectious: An inverse of duration of gametocytogenesis

1/7 day-1


The rate that an asymptomatic human lose the infectivity: An inverse of duration of gametocyte carriage in natural infection

1/188 day-1


The rate that an infectious human lose the infectivity with respect to treatment:


σ1: An inverse of duration of gametocyte carriage after mefloquine treatment

1/10 day-1


σP: An inverse of duration of gametocyte carriage after ACTs treatment

1/6 day-1


σQ: An inverse of duration of gametocyte carriage after primaquine follow-up treatment

1/2 day-1


Import rate of asymptomatic cases: The proportion of imported case per person per day

0.001/365 day-1


Death rate of mosquitoes: An inverse of expected lifespan of a mosquito

1/12 day-1


The rate that an infected mosquito becomes infectious: An inverse of duration of sporogony

1/12 day-1


Average recruitment rate of adult mosquitoes

0.33/12 day-1

Model for gametocyte cycle

The duration of infectivity under different drug regimens was computed using a within-host model of gametocytogenesis (Figure 2). The dynamics of the model are described by the following differential equations.
Figure 2

Schematic illustration of the model for gametocytogenesis process. An asexual parasite produces 16 new parasites at every erythrocytic schizogony cycle. A proportion of these asexual parasites convert to immature gametocytes that appear in the peripheral blood for a day before sequestrating on the blood vessel. The gametocytes become visible again when the maturation process is completed. A combination of schizonticidal and gametocidal treatment affects the gametocyte production process by eliminating the asexual parasites and the mature gametocytes.

In P. falciparum, the erythrocytic stage takes approximately two days, at the end of which each asexual parasite, A g , produces approximately 16 merozoites (θ) [25, 26]. In each cycle of erythrocytic-schizogony, all merozoites were assumed to have an equal chance of undergoing gametocytogenesis; however, only a proportion, p = 0.02, of merozoites commit to gametocytogenesis [27]. The mortality rate (α) is applied to the asexual parasite population due to the schizonticidal treatment. In the model, patients were assumed to receive treatment immediately after asexual parasite density reached 104/μL of blood, the level that normally causes clinical symptoms among people in low-endemic areas [2830]. After artesunate-mefloquine combination treatment, the density of asexual parasites is reduced by a factor of about 1,000 per 2-day schizogony cycle (α), i.e., the asexual parasites density is dropped from 104/μL to 10/μL at the first cycle after treatment [3133]. The surviving merozoites at each cycle then either convert to gametocytes or multiply into newly merozoites that continue to the next schizogony cycle. Therefore, the initial density of early stage gametocytes is a product of the conversion proportion and the net number of merozoites produced with each cycle (). The early stage of immature gametocytes (I g ), which are indistinguishable from asexual parasites, remain in circulation for about one day (δ-1), then sequester on blood vessels while continuing the maturation process (sequestered immature gametocytes; S g ). It takes about 8 days (β-1) for gametocytes to mature and release to the peripheral blood (visible mature gametocytes; G g ). The longevity of a mature gametocyte in the blood stream is approximately 3.5–4 days (ρ-1) [8].

The mortality rate of mature gametocytes due to primaquine treatment (τ) depends on the day of primaquine administration relative to the day of initial schizonticidal treatment (asexual parasites reach 104/μL of blood). Plasma concentration of primaquine was estimated to decrease at elimination rate 0.53 per day (an 8-hour elimination half-life). Since the plasma concentration of primaquine in vivo is difficult to determine, primaquine was assumed to remain effective in killing 90% of mature gametocytes when the plasma concentration was above 10-5 of the maximum concentration, for a net gametocyte reduction ratio of 100 per 2-day cycle, i.e., gametocyte density, which is a function of the mature gametocytes that arise from newly-emerging merozoites at each two-day cycle, eight days earlier, are reduced by 100-fold. The untreated gametocytes remaining in the blood over time were then used to determine the duration of infectiousness. Multiple realizations with different timeframes of primaquine administration were performed to determine the different durations of individual infectiousness after primaquine treatment

The duration of infectiousness is related to mature gametocyte densities and mature gametocyte longevity. The duration of infectiousness was computed based on the density of gametocytes in an individual host, using a log sigmoid relationship between gametocyte density (G g ) and infectivity to mosquitoes [34].
With a 3.6 day half-life, gametocyte densities decline by 90% every eight days. If gametocyte densities were very high initially, then there would be virtually no drop-off in infectivity until the densities reach the S-portion of the logistic curve (Figures 3A and 3B). The time that a person remains infective (probability of infectivity more than zero) differs by about 10 days for every 10-fold difference in the initial gametocyte density (Figure 3C). Thus, for every 10-fold reduction in initial gametocyte densities achieved by curing asexual blood-stage infections alone there will be about 10-day reduction in the duration of human infectivity to mosquito. Therefore, a switch from mefloquine to artemisinin-mefloquine combination therapy (ACT) would reduce infectivity by about 10 days for every 10-fold reduction in mature gametocyte densities. Similarly, a two-day delay in appearing at the clinic could result in a 10-fold increase in gametocyte densities and a 10-day increase in infectiousness. Further reductions can only be achieved by a follow-up treatment with primaquine to reduce the densities of mature gametocytes (Figure 3D).
Figure 3

Relationship among initial gametocyte density, gametocyte longevity, and duration of infectiousness. (A) Changes in gametocyte density over time according to the initial gametocyte density (per μL): 105 (black line), 104 (blue line), 103 (red line). (B) Infectivity to mosquito as a function of gametocyte density (per μL). (C) Probability of infectivity to mosquito over time at different initial gametocyte densities. (D) Duration of infectiousness related to initial gametocyte density without (solid line) and with (dotted line) primaquine follow-up treatment.

Effect of optimally timing primaquine treatment at the population-level

The impact of optimally timing primaquine treatment was assessed by comparing the initial basic reproductive number for malaria derived from the population when no intervention was applied (R 0 ) with the basic reproductive numbers for malaria at different scenarios of the treatment intervention (R c ). R 0 is an estimate of expected number of hosts infected by a single infectious person during his or her entire infectious period [35]. The magnitude of R 0 provides insight into the transmission intensity of the disease and is often used to justify the effect of intervention programmes [36, 37]. The classic formula for R 0 of Ross and Macdonald is shown in equation 5 [35].

Where V denotes vectorial capacity, detailed definitions of other parameters are shown in Table 2.

The basic reproductive number for the control programme was calculated by modifying the Ross and Macdonald formula to consider the proportion treated with ACT (P), or with the follow-up with primaquine (Q):

Where c and r indicate the individual infectivity and duration of gametocyte carriage in natural infections, respectively. Parameters c P and 1/σ P indicate the individual infectivity and duration of gametocyte carriage, respectively, when a schizonticidal drug regimen is applied to a population; the reductions are due to the clearance of asexual parasites. The parameters c Q and 1/σ Q indicate the individual infectivity and duration of gametocyte carriage, respectively, when a primaquine follow-up regimen is applied to the population. A product of c P and 1/σ P , or c Q and 1/σ Q is the cumulative duration of infectiousness for the two different treatment regimens, which is defined as a function of gametocyte density in equation 4.

The impact of the follow-up primaquine regimen was represented by the magnitude of a ratio between R 0 and R c (R 0 /R c ). The vectorial capacity (V) and the transmission probability from mosquitoes to human (b) were assumed to be the same regardless of the intervention, so they cancel out in the ratio. Therefore, R 0 /R C can be calculated by the following equation.

The larger the magnitude of R 0 /R C , the greater the reduction in potential transmission. In addition, the relationship between the percent coverage with the primary treatment (P) and with the follow-up treatment (Q) and the ratio of R 0 and R c was examined.

Finally, the optimal timing of primaquine administration was computed by finding the timing that produced the shortest infectious period in the gametocyte model. By assuming 100% coverage of follow-up primaquine treatment among symptomatic patients, the new infectious period was replaced in the initial transmission model to examine the effect of optimal timing of primaquine administration at the population-level. Also the possibility of malaria elimination in the area was determined when different elimination strategies were implemented.


Transmission pattern

The transmission model provided a good estimate of the dynamics of malaria transmission in the study area. Over the eight-year period, the estimated incidence and observed incidence were comparable in all years except 2000, when the model significantly underestimated malaria incidence (Figure 4A). Malaria transmission in the area was strongly seasonal. Although the model assumed that 99% of the infected population received standard malaria treatment, malaria transmission only decreased gradually in the first three years before settling into a lower, stable orbit, with the annual peak incidence about 1.8 per 100 persons. The switch from mefloquine to ACT in 2005 resulted in a significant reduction in the incidence in year 2006, to less than one per 100 persons. However, estimates from the model suggest that malaria transmission will continue even after the deployment of ACT.
Figure 4

Plasmodium falciparum incidence and dynamic changes in proportion of the human population. (A) Incidence of P.falciparum malaria in seven villages of Ratchaburi province. Solid line represents the observed incidence from 1999 to 2006. The incidence from 1999 to 2009, estimated by the transmission model, is shown by the blue dashed line. (B) Estimated changes in proportion of population in each human compartment over the 10-year period. Solid line indicates population with clinical symptoms, C h ; Red dashed line represents population with gametocytaemia, G h ; Blue dotted line indicates asymptomatic population, A h .

The changes in transmission dynamics are illustrated in Figure 4B. As is typical of low malaria transmission settings, the proportion of asymptomatic gametocyte carriers was low compared with the proportion of symptomatic individuals, which varied seasonally. While the proportion of people with gametocytaemia almost reached zero during the dry season, the model predicts that about three percent of the population remained gametocytaemic when the environment was suitable for mosquito vectors. This small proportion of gametocytaemic people could play an important role in maintaining transmission of the parasite in the area.

Gametocyte carriage regarding primaquine treatment

Different malaria dynamics were observed when primaquine treatment was given at different times. Changes in the density of each parasite developmental stage over time with different primaquine treatment regimens are shown in Figure 5A. The simulated primaquine plasma concentration remained above the killing concentration for up to three days after administration. The duration and density of individual gametocytaemia differed substantially according to the timing of primaquine administration. When primaquine was given on the same day as other anti-malarial drugs, there were no mature gametocytes in the blood to be killed by primaquine, and infectiousness of the patient was not changed in comparison with no primaquine treatment. The duration of individual infectivity to mosquitoes could be as long as 14 days. In contrast, duration of infectiousness would be greatly reduced to two days if the primaquine treatment was delayed until the majority of mature gametocytes were circulating in the bloodstream (day 7 or later). On the other hand, delaying primaquine administration until day 10 or later results in an increased duration of infectiousness as mature gametocytes circulate for several days (Figure 5B).
Figure 5

Gametocyte dynamics and duration of infectiousness with different timings of primaquine administration. (A) Changes in the density of each parasite developmental stage over time with different timings of primaquine administration. Pharmacokinetic curve of primaquine (dotted line) related to waves of asexual parasite, A g (dashed line), immature gametocytes, I g (grey line), and mature gametocytes, M g (solid red line). Initial asexual parasite (A g ) 10/μL (B) Duration of infectivity of humans to mosquitoes over different timings of primaquine administration after acute attack.

Optimal timing of primaquine treatment and malaria transmission

The optimal timing of primaquine treatment showed a substantial effect in reducing malaria transmission in the area. However, the added value of the follow-up primaquine treatment was strongly correlated with the proportion of patients treated with anti-malarial drugs (P). The follow-up primaquine treatment showed a greater effect when a large fraction of clinical malaria was detected and treated. In an area with an effective detection and treatment control programme (P = 99%), the basic reproductive rate (R c ) could be reduced up to 35 times by implementing the optimally timed primaquine treatment programme, compared with the initial basic reproductive rate (R 0 ). More importantly, the proportional reduction in transmission that was achieved through an 8-day follow-up with primaquine in 99% of the patients was approximately equal to the proportional reductions achieved from the initial treatment. In addition, the follow-up primaquine treatment programme appeared to be effective at reducing transmission only when the coverage of the follow-up primaquine treatment was above 90% (Figure 6).
Figure 6

The R 0 /Rc ratio for ACT and follow-up primaquine treatment. The relationship between R 0 /R c and the proportion of P. falciparum patients receiving follow-up primaquine treatment (Q) for different proportions of patients treated with artesunate-mefloquine combination therapy (P).

The transmission model was re-simulated by replacing the duration of infectivity (σ-1) from 10 days to 2 days, according to the gametocyte model. The model indicated that if all individuals received primaquine at day 8 after initial acute attack, the malaria incidence in the area would be reduced substantially (Figure 7). However, changing the primaquine administration regimen alone may not be enough to eliminate P. falciparum malaria in the area; the incidence of malaria still persisted during the wet season, largely because of imported malaria. By including an intervention focused on the effective detection and treatment of imported asymptomatic infections, the incidence of malaria in the area could reach zero within three years after the combination programme was introduced (Figure 7).
Figure 7

Estimated P. falciparum malaria incidence for different control policies. Solid line represents the incidence estimated by the initial model when the standard drug regimen (2-day ACTs and Primaquine at day 2) was used (solid line). The incidence significantly decreased when the timing of primaquine administration was shifted to the eighth day after initial attack (short-dashed line). The elimination of malaria transmission can be reached when the combination of both primaquine follow-up treatment and control of imported asymptomatic cases was implemented (long-dashed line).


Mathematical models were constructed to understand the transmission dynamics of P. falciparum malaria in an area where access to health-care has significantly reduced malaria transmission. Active malariometric surveys suggest that a very high fraction of clinical episodes in the area, perhaps higher than 97% are promptly treated [38]. Prompt and effective treatment may be a very cost-effective strategy for malaria control in low or unstable malaria transmission settings, because most individuals are likely to develop acute febrile illness after P. falciparum infection [3941]. Significant reductions in mortality and morbidity of malaria after deploying a strategy of early detection and prompt treatment has been documented [42]. However, model simulations find that, despite significant reductions in incidence, P. falciparum transmission is likely to continue. This suggests that the elimination of malaria may be difficult even in an area where malaria health systems are highly effective.

Given the fact that multiple surveys of large swaths of the population are unable to find individuals infected asymptomatically, gametocyte carriage in the absence of asexual parasites seems to be vital for maintaining transmission in this heavily controlled area [7, 24]. Findings of our model suggest that there is a strong seasonal fluctuation in the population of residual gametocyte carriers, which is consistent with previous observations in other low transmission areas [7, 9]. An interesting finding is that while the prevalence of gametocytes is relatively small during the dry season, it does not drop to zero. Presumably then these cases are responsible for the source of mosquito infection at the start of wet season [7, 9].

Gametocyte reduction is of great interest for malaria control, particularly in low-endemic malaria areas. One of the strategies for the malaria elimination programme, recommended by the WHO, is to identify and treat all malaria patients as well as to reduce onward transmission caused by gametocytaemia [6]. The policy change from mefloquine treatment to artesunate-mefloquine has been shown to reduce transmission. Artemisinin combination therapies can reduce the asexual parasite burden 100 times faster than mefloquine, which can subsequently inhibit development of more mature gametocytes [3133, 41]. The gametocyte model indicates that there is still an added value of the follow-up primaquine treatment even when the initial gametocyte density is small due to the switch from mefloquine to ACTs. In addition, in areas where artemisinin combination therapy has been used, gametocyte carriage is still common in the 7–21 days following treatment [13, 30, 41, 4345]. In many countries, a single oral dose of primaquine is included in the standard anti-malarial drug regimen with the aim of further reducing gametocyte carriage, even when artemisinin-based therapy is used [6, 12].

However, while in these areas primaquine can be extremely effective at clearing gametocytes that persist after treatment with schizonticidal agents [13, 14, 43], the timing and duration of gametocyte carriage and subsequent infectiousness have not been considered carefully when primaquine is deployed as a transmission-blocking agent. Findings of the gametocyte model indicate that the effectiveness of primaquine in reducing the duration of infectiousness depends critically on timing. Primaquine is most beneficial when the administration is delayed, about eight days following initial treatment, to coincide with the release of a large cohort of mature gametocytes into the blood, which emerged from a large number of merozoites during an acute attack. The effect of primaquine is significantly reduced when the drug is given too early or too late. Although an immediate primaquine treatment can affect a small cohort of mature gametocytes that emerge from the first crop of merozoites that appear in circulation at the time of an acute attack, primaquine will be cleared from the system before the largest cohort of gametocytes mature. If primaquine is given too late, mature gametocytes will be able to circulate and infect mosquitoes until the drug is administered.

The benefits of optimally timed primaquine are greatest in those areas where the early treatment programme to cure asexual blood stage infections is very successful; a high fraction of clinical malaria episodes are expected to receive the standard treatment within one to three days after acute attack. In such areas, optimally timing primaquine administration shows a potential impact on overall malaria transmission at the population-level. The current results show that follow-up primaquine treatment can reduce the duration of infectiousness over the existing strategy of using artesunate-mefloquine alone, with a combined total net reduction in transmission of 98%, a 95-fold reduction in R 0 . An important observation is that the added value of optimally timed primaquine can have relatively large effects on reducing transmission only if a high fraction of patient infections are treated and cured with first-line anti-malarial drugs (i.e. when P is high), suggesting that the first emphasis should be on treating those with clinical malaria. Because primaquine effectively reduces transmission only in those patients who have cleared their asexual parasites, and because the average duration of an asymptomatic infection is approximately six months, the benefit of reducing the duration of gametocyte carriage is of little importance unless at least 90% of clinical malaria episodes are effectively treated. Treatment to clear asexual parasites and prevent asymptomatic infections can only reduce the duration of infectiousness insofar as the gametocytes are also cleared. In such situations, primaquine can be very effective at further reducing the duration of infectiousness, and the added value of good follow-up with primaquine treatment has nearly the same proportional effects on potential transmission as does the primary treatment.

In addition, findings from the model indicate that when all symptomatic P. falciparum patients receive the follow-up primaquine treatment at day eight, the P. falciparum malaria incidence can be reduced nearly to zero. However, the model suggests that a small number of undetected imported cases can pose a big threat for malaria elimination. To reach the elimination goal, vigilance to detect and cure imported asymptomatic cases may also be required. These findings support the WHO recommendations for a malaria elimination programme [6].

Mathematical modeling has been widely used to model transmission dynamics of malaria, and control interventions [46]. However, findings from the models require a careful interpretation. This study intends to construct simple models that provide a valuable insight into the feasibility of malaria elimination in a low malaria transmission area. Results from the models should be considered as approximations that are likely to differ because of the natural variability under field conditions.

In addition, although the model assumed that the level of individual infectiousness follows the log-sigmoid relationship with the gametocyte density among non-immune adults, the duration of infectiousness in the analysis may be underestimated. Infectivity to mosquitoes is observed even when gametocyte densities fall below detection level by microscopy or by a molecular method [24, 47]. However, the trend of duration of infectiousness over different timings of primaquine administration does not change when different infectivity levels are applied. The model also does not take into account the additional gametocyte carriage from recrudescent infections. Though gametocytaemia is estimated to be greater in recrudescent infections than in primary infections, in an area where artesunate-mefloquine combination therapy is used, few recrudescent infections are expected [33, 48]. Lastly, the model assumes that most infected individuals develop clinical symptoms and are treated. In low and unstable transmission areas this assumption is generally correct; however pockets of higher transmission may exist, and the importance of asymptomatic asexual blood-stage infections in these areas in continuing transmission over the dry season can be significant. While surveys suggest that there are basically no asymptomatic carriers in the region in question, studies in other areas have found that sub-patent infections can persist for many months [9, 10, 49]. Thus, while the models suggest that optimally-timed primaquine administration can significantly impact the incidence of malaria in a low-transmission area well served by health centers, asymptomatic individuals in the area and not just imported carriers may also play a significant role in sustaining transmission and should be considered in any elimination plan.


Mathematical models constructed in this study pose an important and testable hypothesis regarding existing control programmes in areas with good health systems where malaria transmission persists. The transmission-blocking effect of primaquine and the timing of its administration should be carefully scrutinized. Given the risks associated with primaquine, it may not be worth giving, except in areas where the early detection and prompt treatment programme is highly effective. In such areas, primaquine should be administered at an appropriate time, or a long acting 8-aminoquinoline should be considered, and combined with surveillance to catch imported malaria it could lead to local malaria elimination. Randomized controlled trials are recommended to determine the most-effective timing of primaquine administration in order to decrease malaria transmission, which is important for planning malaria elimination programmes in low malaria transmission areas.



We would like to thank Dr. Christopher V. Plowe (Department of Medicine, University of Maryland, Baltimore) for helping to review the manuscript. We also thank all officers and staff at Rajanagarindra Tropical Disease International Center, Suan Phung, Ratchaburi, Thailand for their help and support. DLS is supported by a grant from the Bill & Melinda Gates Foundation (#49446) and funding from the RAPIDD program of the Science & Technology Directorate, Department of Homeland Security, and the Fogarty International Center, National Institutes of Health.

Authors’ Affiliations

Department of Epidemiology and Preventive Medicine, University of Maryland School of Medicine
Department of Tropical Hygiene, Faculty of Tropical Medicine, Mahidol University
Department of Ecology and Evolutionary Biology, Princeton University
Department of Biology, University of Florida
Emerging Pathogens Institute, University of Florida


  1. World Health Organization and UNICEF: World malaria report 2005. Geneva. 2005
  2. Guerra CA, Gikandi PW, Tatem AJ, Noor AM, Smith DL, Hay S, Snow RW: The limits and intensity of Plasmodium falciparum transmission: Implications for malaria control and elimination worldwide. PLoS Med. 2008, 5: e38-10.1371/journal.pmed.0050038.PubMed CentralView ArticlePubMed
  3. Hay SI, Guerra CA, Tatem AJ, Noor AM, Snow RW: The global distribution and population at risk of malaria: Past, present, and future. Lancet Infect Dis. 2004, 4: 327-336. 10.1016/S1473-3099(04)01043-6.PubMed CentralView ArticlePubMed
  4. Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI: The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005, 434: 214-217. 10.1038/nature03342.PubMed CentralView ArticlePubMed
  5. Guerra CA, Snow RW, Hay SI: Mapping the global extent of malaria in 2005. Trends Parasitol. 2006, 22: 353-358. 10.1016/ CentralView ArticlePubMed
  6. Global Malaria Programme: Malaria elimination: A field manual for low and moderate endemic countries. 2007, Geneva, World Health Organization
  7. Drakeley C, Sutherland C, Bousema JT, Sauerwein RW, Targett GA: The epidemiology of Plasmodium falciparum gametocytes: Weapons of mass dispersion. Trends Parasitol. 2006, 22: 424-430. 10.1016/ ArticlePubMed
  8. Talman AM, Domarle O, McKenzie FE, Ariey F, Robert V: Gametocytogenesis: The puberty of Plasmodium falciparum . Malar J. 2004, 3: 24-10.1186/1475-2875-3-24.PubMed CentralView ArticlePubMed
  9. Abdel-Wahab A, Abdel-Muhsin AM, Ali E, Suleiman S, Ahmed S, Walliker D, Babiker HA: Dynamics of gametocytes among Plasmodium falciparum clones in natural infections in an area of highly seasonal transmission. J Infect Dis. 2002, 185: 1838-1842. 10.1086/340638.View ArticlePubMed
  10. Nassir E, Abdel-Muhsin AM, Suliaman S, Kenyon F, Kheir A, Geha H, Ferguson HM, Walliker D, Babiker HA: Impact of genetic complexity on longevity and gametocytogenesis of Plasmodium falciparum during the dry and transmission-free season of eastern Sudan. Int J Parasitol. 2005, 35: 49-55. 10.1016/j.ijpara.2004.10.014.View ArticlePubMed
  11. Pethleart A, Prajakwong S, Suwonkerd W, Corthong B, Webber R, Curtis C: Infectious reservoir of Plasmodium infection in Mae Hong Son province, north-west Thailand. Malar J. 2004, 3: 34-10.1186/1475-2875-3-34.PubMed CentralView ArticlePubMed
  12. World Health Organization: Guidelines for the treatment of malaria. Geneva. 2006
  13. Pukrittayakamee S, Chotivanich K, Chantra A, Clemens R, Looareesuwan S, White NJ: Activities of artesunate and primaquine against asexual- and sexual-stage parasites in falciparum malaria. Antimicrob Agents Chemother. 2004, 48: 1329-1334. 10.1128/AAC.48.4.1329-1334.2004.PubMed CentralView ArticlePubMed
  14. Kamtekar KD, Gogtay NJ, Dalvi SS, Karnad DR, Chogle AR, Aigal U, Kshirsagar NA: A prospective study evaluating the efficacy of a single, 45-mg dose of primaquine, as a gametocytocidal agent, in patients with Plasmodium falciparum malaria in Mumbai, India. Ann Trop Med Parasitol. 2004, 98: 453-458. 10.1179/000349804225003550.View ArticlePubMed
  15. Elmes NJ, Bennett SM, Abdalla H, Carthew TL, Edstein MD: Lack of sex effect on the pharmacokinetics of primaquine. Am J Trop Med Hyg. 2006, 74: 951-952.PubMed
  16. Giao PT, Vries PJ, Binh TQ, Nam NV, Kager PA: Early diagnosis and treatment of uncomplicated malaria and patterns of health seeking in Vietnam. Trop Med Int Health. 2005, 10: 919-925. 10.1111/j.1365-3156.2005.01472.x.View ArticlePubMed
  17. Ministry of Public Health: Thailand. Malaria control programme in Thailand. Accessed: 1 November 2008, [http://​eng.​moph.​go.​th/​SpecificHealth/​malaria/​malaria.​htm]
  18. Buchachart K, Krudsood S, Singhasivanon P, Treeprasertsuk S, Phophak N, Srivilairit S, Chalermrut K, Rattanapong Y, Supeeranuntha L, Wilairatana P, Brittenham G, Looareesuwan S: Effect of primaquine standard dose (15 mg/day for 14 days) in the treatment of vivax malaria patients in Thailand. Southeast Asian J Trop Med Public Health. 2001, 32: 720-726.PubMed
  19. Bangchang KN, Songsaeng W, Thanavibul A, Choroenlarp P, Karbwang J: Pharmacokinetics of primaquine in G6PD deficient and G6PD normal patients with vivax malaria. Trans R Soc Trop Med Hyg. 1994, 88: 220-222. 10.1016/0035-9203(94)90306-9.View ArticlePubMed
  20. Garros C, Van Bortel W, Trung HD, Coosemans M, Manguin S: Review of the minimus complex of anopheles, main malaria vector in Southeast Asia: From taxonomic issues to vector control strategies. Trop Med Int Health. 2006, 11: 102-114. 10.1111/j.1365-3156.2005.01536.x.View ArticlePubMed
  21. Sungvornyothin S, Muenvorn V, Garros C, Manguin S, Prabaripai A, Bangs MJ, Chareonviriyaphap T: Trophic behavior and biting activity of the two sibling species of the Anopheles minimus complex in western Thailand. J Vector Ecol. 2006, 31: 252-261. 10.3376/1081-1710(2006)31[252:TBABAO]2.0.CO;2.View ArticlePubMed
  22. Trampuz A, Jereb M, Muzlovic I, Prabhu RM: Clinical review: Severe malaria. Crit Care. 2003, 7: 315-323. 10.1186/cc2183.PubMed CentralView ArticlePubMed
  23. Suputtamongkol Y, Chindarat S, Silpasakorn S, Chaikachonpatd S, Lim K, Chanthapakajee K, Kaewkaukul N, Thamlikitkul V: The efficacy of combined mefloquine-artesunate versus mefloquine-primaquine on subsequent development of Plasmodium falciparum gametocytemia. Am J Trop Med Hyg. 2003, 68: 620-623.PubMed
  24. Shekalaghe SA, Bousema JT, Kunei KK, Lushino P, Masokoto A, Wolters LR, Mwakalinga S, Mosha FW, Sauerwein RW, Drakeley CJ: Submicroscopic Plasmodium falciparum gametocyte carriage is common in an area of low and seasonal transmission in Tanzania. Trop Med Int Health. 2007, 12: 547-553.View ArticlePubMed
  25. Miller LH, Baruch DI, Marsh K, Doumbo OK: The pathogenic basis of malaria. Nature. 2002, 415: 673-679. 10.1038/415673a.View ArticlePubMed
  26. Diebner HH, Eichner M, Molineaux L, Collins WE, Jeffery GM, Dietz K: Modelling the transition of asexual blood stages of Plasmodium falciparum to gametocytes. J Theor Biol. 2002, 202: 113-127. 10.1006/jtbi.1999.1041.View Article
  27. McKenzie FE, Bossert WH: The optimal production of gametocytes by Plasmodium falciparum. J Theor Biol. 1998, 193: 419-428. 10.1006/jtbi.1998.0710.View ArticlePubMed
  28. Hung le Q, de Vries PJ, Binh TQ, Giao PT, Nam NV, Holman R, Kager PA: Artesunate with mefloquine at various intervals for non-severe Plasmodium falciparum malaria. Am J Trop Med Hyg. 2004, 71: 160-166.PubMed
  29. Akim NI, Drakeley C, Kingo T, Simon B, Senkoro K, Sauerwein RW: Dynamics of P. falciparum gametocytemia in symptomatic patients in an area of intense perennial transmission in Tanzania. Am J Trop Med Hyg. 2000, 63: 199-203.PubMed
  30. Tangpukdee N, Krudsood S, Srivilairit S, Phophak N, Chonsawat P, Yanpanich W, Kano S, Wilairatana P: Gametocyte clearance in uncomplicated and severe Plasmodium falciparum malaria after artesunate-mefloquine treatment in Thailand. Korean J Parasitol. 2008, 46: 65-70. 10.3347/kjp.2008.46.2.65.PubMed CentralView ArticlePubMed
  31. White NJ: Assessment of the pharmacodynamic properties of antimalarial drugs in vivo. Antimicrob Agents Chemother. 1997, 41: 1413-1422.PubMed CentralPubMed
  32. White NJ: Qinghaosu (artemisinin): The price of success. Science. 2008, 320: 330-334. 10.1126/science.1155165.View ArticlePubMed
  33. Nosten F, White NJ: Artemisinin-based combination treatment of falciparum malaria. Am J Trop Med Hyg. 2007, 77 (6 Suppl ): 181-192.PubMed
  34. Barnes KI, White NJ: Population biology and antimalarial resistance: The transmission of antimalarial drug resistance in Plasmodium falciparum. Acta Trop. 2005, 94: 230-240.View ArticlePubMed
  35. Macdonald G: The epidemiology and control of malaria. 1957, London, United Kingdom: Oxford University Press
  36. Heffernan JM, Smith RJ, Wahl LM: Perspectives on the basic reproductive ratio. J R Soc Interface. 2005, 2: 281-293. 10.1098/rsif.2005.0042.PubMed CentralView ArticlePubMed
  37. Smith DL, McKenzie FE, Snow RW, Hay SI: Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 2007, 5: e42-10.1371/journal.pbio.0050042.PubMed CentralView ArticlePubMed
  38. Phimpraphi W, Paul RE, Yimsamran S, Puangsa-art S, Thanyavanich N, Maneeboonyang W, Prommongkol S, Sornklom S, Chaimungkun W, Chavez IF, Blanc H, Looareesuwan S, Sakuntabhai A, Singhasivanon P: Longitudinal study of Plasmodium falciparum and plasmodium vivax in a Karen population in Thailand. Malar J. 2008, 7: 99-10.1186/1475-2875-7-99.PubMed CentralView ArticlePubMed
  39. Breman JG: The ears of the hippopotamus: Manifestations, determinants, and estimates of the malaria burden. Am J Trop Med Hyg. 2001, 64 (1-2 Suppl ): 1-11.PubMed
  40. Guerin PJ, Olliaro P, Nosten F, Druilhe P, Laxminarayan R, Binka F, Kilama WL, Ford N, White NJ: Malaria: Current status of control, diagnosis, treatment, and a proposed agenda for research and development. Lancet Infect Dis. 2002, 2: 564-573. 10.1016/S1473-3099(02)00372-9.View ArticlePubMed
  41. Stepniewska K, Price RN, Sutherland CJ, Drakeley CJ, von Seidlein L, Nosten F, White NJ: Plasmodium falciparum gametocyte dynamics in areas of different malaria endemicity. Malar J. 2008, 7: 249-10.1186/1475-2875-7-249.PubMed CentralView ArticlePubMed
  42. Carrara VI, Sirilak S, Thonglairuam J, Rojanawatsirivet C, Proux S, Gilbos V, Brockman A, Ashley EA, McGready R, Krudsood S, Leemingsawat S, Looareesuwan S, Singhasivanon P, White N, Nosten F: Deployment of early diagnosis and mefloquine-artesunate treatment of falciparum malaria in Thailand: The Tak malaria initiative. PLoS Med. 2006, 3: e183-10.1371/journal.pmed.0030183.PubMed CentralView ArticlePubMed
  43. Shekalaghe S, Drakeley C, Gosling R, Ndaro A, van Meegeren M, Enevold A, Alifrangis M, Mosha F, Sauerwein R, Bousema T: Primaquine clears submicroscopic Plasmodium falciparum gametocytes that persist after treatment with sulphadoxine-pyrimethamine and artesunate. PLoS ONE. 2007, 2: e1023-10.1371/journal.pone.0001023.PubMed CentralView ArticlePubMed
  44. Chotivanich K, Sattabongkot J, Udomsangpetch R, Looareesuwan S, Day NP, Coleman RE, White NJ: Transmission-blocking activities of quinine, primaquine, and artesunate. Antimicrob Agents Chemother. 2006, 50: 1927-1930. 10.1128/AAC.01472-05.PubMed CentralView ArticlePubMed
  45. Piyaphanee W, Krudsood S, Tangpukdee N, Thanachartwet W, Silachamroon U, Phophak N, Duangdee C, Haoharn O, Faithong S, Wilairatana P, Leowattana W, Looareesuwan S: Emergence and clearance of gametocytes in uncomplicated Plasmodium falciparum malaria. Am J Trop Med Hyg. 2006, 74: 432-435.PubMed
  46. McKenzie FE: Why model malaria?. Parasitol Today. 2000, 16: 511-516. 10.1016/S0169-4758(00)01789-0.View ArticlePubMed
  47. Schneider P, Bousema JT, Gouagna LC, Otieno S, Vegte-Bolmer van de M, Omar SA, Sauerwein RW: Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infection. Am J Trop Med Hyg. 2007, 76: 470-474.PubMed
  48. Nosten F, van Vugt M, Price R, Luxemburger C, Thway KL, Brockman A, McGready R, ter Kuile F, Looareesuwan S, White NJ: Effects of artesunate-mefloquine combination on incidence of Plasmodium falciparum malaria and mefloquine resistance in western Thailand: A prospective study. Lancet. 2000, 356: 297-302. 10.1016/S0140-6736(00)02505-8.View ArticlePubMed
  49. Babiker HA, Abdel-Muhsin AM, Ranford-Cartwright LC, Satti G, Walliker D: Characteristics of Plasmodium falciparum parasites that survive the lengthy dry season in eastern Sudan where malaria transmission is markedly seasonal. Am J Trop Med Hyg. 1998, 59: 582-590.PubMed


© Lawpoolsri et al; licensee BioMed Central Ltd. 2009

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.