 Research article
 Open Access
 Open Peer Review
 Published:
Populationlevel effects of clinical immunity to malaria
BMC Infectious Diseasesvolume 13, Article number: 428 (2013)
Abstract
Background
Despite a resurgence in control efforts, malaria remains a serious publichealth problem, causing millions of deaths each year. One factor that complicates malariacontrol efforts is clinical immunity, the acquired immune response that protects individuals from symptoms despite the presence of parasites. Clinical immunity protects individuals against disease, but its effects at the population level are complex. It has been previously suggested that under certain circumstances, malaria is bistable: it can persist, if established, in areas where it would not be able to invade. This phenomenon has important implications for control: in areas where malaria is bistable, if malaria could be eliminated until immunity wanes, it would not be able to reinvade.
Methods
Here, we formulate an analytically tractable, dynamical model of malaria transmission to explore the possibility that clinical immunity can lead to bistable malaria dynamics. We summarize what is known and unknown about the parameters underlying this simple model, and solve the model to find a criterion that determines under which conditions we expect bistability to occur.
Results
We show that bistability can only occur when clinically immune individuals are more “effective” at transmitting malaria than naïve individuals are. We show how this “effectiveness” includes susceptibility, ability to transmit, and duration of infectiousness. We also show that the amount of extra effectiveness necessary depends on the ratio between the duration of infectiousness and the time scale at which immunity is lost. Thus, if the duration of immunity is long, even a small amount of extra transmission effectiveness by clinically immune individuals could lead to bistability.
Conclusions
We demonstrate a simple, plausible mechanism by which clinical immunity may be causing bistability in human malaria transmission. We suggest that simple summary parameters – in particular, the relative transmission effectiveness of clinically immune individuals and the time scale at which clinical immunity is lost – are key to determining where and whether bistability is happening. We hope these findings will guide future efforts to measure transmission parameters and to guide malaria control efforts.
Background
Despite extensive efforts to eradicate it, malaria caused by Plasmodium falciparum remains a significant problem resulting in millions of cases and 660,000 deaths in 2010 [1]. A characteristic of falciparum malaria disease that complicates control efforts is clinical immunity – an immune response that develops with exposure to parasites and provides protection against the clinical symptoms of malaria, despite the presence of parasites [2]. Although clinical immunity protects individuals against disease, its effects at the population level are complex.
Malaria is highly variable from region to region, further complicating analysis. Geographic variation in average disease burden (endemicity) leads to variation in acquired immunity [3]. Malaria endemicity ranges from “holoendemic” (defined as having a parasite ratio (PR, the percentage of subjects with parasites found in the blood) consistently greater than 75% of infants [4]) through “hyperendemic” and “mesoendemic” to “hypoendemic”, defined as having a PR of less than 10% of children age 2–9 [4]. This results in variation in the acquisition of clinical immunity. This variation in endemicity and clinical immunity complicates malaria epidemiology and control.
Clinical immunity to malaria develops after exposure to parasites and varies as endemicity varies [5, 6]. In holoendemic regions, exposure to parasitemia is high enough that clinical immunity develops rapidly, and most adults and older children are clinically immune, whereas in hypoendemic regions, most people are not reinfected often enough to develop clinical immunity [6]. Even after it develops, clinical immunity can be lost in 35 years without reexposure [2, 3, 5–7]. When individuals first develop clinical immunity, they are only immune to severe symptoms. If reexposure continues, however, clinical immunity can result in asymptomatic or nearly asymptomatic disease. Full clinical immunity develops slowly and tends to correlate with the onset of puberty [6, 8].
An important aspect of clinical immunity is the possibility that clinically immune individuals are particularly effective at transmitting malaria over the duration of their infections. This phenomenon could arise if clinically immune individuals are more infectious to mosquitoes per unit time, or if they stay infectious for longer (perhaps because they are less likely to seek medical treatment), or both. If clinically immune individuals have a higher reproductive rate, this has potential implications at the population level – in some cases, malaria may be spread more effectively in areas where it is already present, all else being equal.
As transmission of malaria decreases, the proportion of the population protected by clinical immunity decreases as well, since clinical immunity is lost. As a result, decreases in transmission can, under some circumstances, lead to an increase in morbidity and mortality, because fewer people are protected against the symptoms of malaria [9–11].
Águas et al. [12] have shown that under certain circumstances when clinically immune individuals are more infectious over the duration of their infection, than naïve cases, malaria can persist, if established, in areas where it would not be able to invade. In other words, for some sets of parameters, both an endemic equilibrium and a diseasefree equilibrium are stable – a phenomenon known as bistability. Bistability would have important implications for malaria control: in particular, it would imply that there are some areas where, if malaria could be eliminated until clinical immunity wanes, it would not be able to reinvade. Here we use a model of malaria transmission to explore under which conditions we would expect bistability to occur, indicating possible opportunities for malaria elimination.
Malaria elimination has been surprisingly effective in many countries: 75 of the 79 countries that successfully eliminated malaria between 1945 and 2010 remain malaria free, even though many have not sustained control efforts [13]. A recent paper by Chiyaka et al. [14] presents six hypotheses for this phenomenon, and argues that ${\mathcal{R}}_{0}$ may be reduced either by external factors, like demographic and hydrological changes, or by factors driven by malaria elimination itself, for example economic development catalyzed by reduced disease burden, or bistability due to treatment seeking. They argue that, to the extent that malaria elimination reduces ${\mathcal{R}}_{0}$, incentives to aggressively pursue control are increased, since ongoing active control efforts will not be required once malaria is eliminated.
In this paper, we first review what is known about infectiousness and susceptibility to infection of clinically immune individuals. We then build a simple transmission model designed to elucidate what factors make bistability likely, and what measurements could shed light on when and whether bistability is likely to be an important phenomenon in malaria dynamics and control. We derive a simple mathematical criterion for how “effective” transmission by clinically immune individuals must be for bistability to occur.
Effects of clinical immunity on disease transmission
The overall infectiousness of an infected individual is the product of duration of infection, and mean infectiousness. Below, we review what is known about the effects of clinical immunity on these components.
The duration of a malaria infection is highly variable, and treatmentseeking behavior is an important determinant. In a review of populationlevel studies done on malaria treatmentseeking behavior, McCombie [15] found that treatment rates were correlated with severity of symptoms, and that in Africa, 64–95% of individuals who sought treatment received at least one form of treatment; with the majority of studies reviewed reporting over 90% treatment rate.
It is well known that many clinically immune infections are not even recognized by the individual as malaria. Individuals in hyper and holoendemic areas who do not think they have malaria, have been found to test positive at high rates, for example in Ghana [15, 16], Senegal [17], and Kenya [18]. Thus, it seems reasonable to suppose that most clinically immune infections are untreated, and last longer than treated clinical infections because of treatmentseeking behavior. Various studies have been done to estimate the duration of untreated malaria infection. Earle et al. [19] observed the duration of infection in children age 515 years old and found all of them had cleared the infection within a year. However, most of what is known about duration of infection comes from malariatherapy data. These studies infected malarianaïve syphilis patients with P. falciparum strains with low clinical virulence and found the mean duration of infection to be 200 days [20]. A recent study in a highly endemic region found a similar mean duration of infection to that of the malariatherapy data, however, they found a larger variance in duration of infection with many more infections with shorter duration than found in the malariatherapy data [21]. And in an extreme case, an infection was found 8 years after last known exposure to parasites [22]. However, it’s not clear how the duration of untreated clinically immune infections compares to untreated symptomatic cases. Bruce et al. [23] found that “episodes” of parasitemia lasted longer in children than in adults and suggest that this may be due to clinical immunity [23]. Unfortunately, it is not straightforward to relate their measured episodes to infection clearance (partly because infected individuals may be “superinfected” by other strains).
While clinically immune cases are frequently asymptomatic or of low clinical virulence, data from malariatherapy studies suggest naïve cases may also have a wide range of clinical virulence ranging from high virulence to asymptomatic [24]. Other studies suggest that asymptomatic malaria is not limited to areas of high transmission where exposurerelated immunity is expected to develop [25–27]. Thus, not all naïve cases may be terminated with treatment.
Another key component to the populationlevel effects of clinical immunity is the infectiousness per transmission event of infected individuals. Although malaria transmission is much studied, it remains unclear how parasitemia, gametocytemia, and other factors interact to affect malaria transmission. Gametocyte quantity alone is not sufficient to ensure successful transmission; mosquito uptake of gametocytes depends on a wide variety of factors, including transmissionblocking immunity (TBI) [28–30] and cytokine tumor necrosis factor (TNF) [31, 32].
Transmissionblocking immunity (TBI) is a human immune reaction to sexual stages of malaria. TBI develops with exposure to gametocytes and, through a variety of mechanisms, reduces successful transmission of new infections. TBI increases with gametocyte density; consequently, TBI tends to be negatively correlated with clinical immunity [30], but the importance of TBI to populationlevel transmission is not clear. In one study, transmissionblocking immunity was found to reduce transmission by up to 90%, with higher immunity in the younger age groups [30]. Two other studies, which did not include the youngest age group, failed to find correlations between age and TBI [33, 34]; the latter of these found that only 15% of urban and 29% of rural gametocyte carriers had reduced transmission. A model of human infectiousness to mosquitoes found that patterns of EIR across Africa and Papua New Guinea could be explained without invoking TBI [35].
Cytokine TNF is another factor that affects gametocyte success. It is present in the blood serum taken during the crisis of a malaria infection [31, 32]. Cytokine TNF is responsible for the loss of infectiousness during peak parasitemia by killing the gametocytes. Gametocytes present during malaria crisis were found dead before entering the mosquito; as well, gametocytes from crisis serum failed to infect mosquitoes even when washed and resuspended in normal serum [31].
Although there is immunity against gametocytes at peak parasitemia in nonimmune individuals, clinically immune infections often have lower gametocytemia as a result of having antiparasite immunity, conferring protection against highdensity parasitemia [7]. There is also evidence that clinically immune individuals are less infectious per bite, but because clinically immune individuals are less likely to seek treatment, they are consistently infectious at low levels for long periods of time and therefore result in producing a large number of new infections over the course of a single clinically immune infectious period [36, 37].
As we will see below, susceptibility to infection in clinically immune individuals is also important to the populationlevel dynamics of malaria. Although a great deal is known about susceptibility to clinical illness, parasitemia or gametocytemia [38, 39], much less is known about susceptibility to new infection. Individual susceptibility to new infections is complex, and known to be influenced by genotype, parasite virulence, and specific immunity [38], but there is evidence to believe that clinically immune individuals are about as susceptible to disease as nonimmune individuals [40].
Populationlevel effects
We investigate the factors underlying bistability with a simple transmission model that accounts for clinical immunity (see Methods). We assume that individuals infected when they are naïve have a probability of becoming clinically immune when they recover from infection, and that clinically immune susceptible individuals lose immunity at some rate if not infected again, meaning that clinical immunity will be maintained when the force of infection is high, and will often wane if the force of infection is low.
To explore the effects that clinically immune individuals have on the populationlevel disease dynamics, we compare the lifecycle transmission effectiveness of naïve and clinically immune individuals, using subgroup reproductive numbers. These are defined as the average number of secondary infections from a single infectious individual in an otherwise totally susceptible population. We define the reproductive number for naïve cases, ${\mathcal{R}}_{\mathit{\text{NN}}}$, as the average number of secondary infections generated by a single naïvely infectious individual in an otherwise totally naïvely susceptible population. We define the reproductive number for clinically immune cases, ${\mathcal{R}}_{\mathit{\text{CC}}}$, as the average number of secondary infections given by a single clinically immune infectious individual in an otherwise totally clinically immune population.
Because we assume that all individuals are naïve in the absence of infection, the basic reproductive number of our system ${\mathcal{R}}_{0}={\mathcal{R}}_{\mathit{\text{NN}}}$. We stress that this is the basic reproductive number in the presence of baseline control efforts – in particular, we assume that treatment is always available to those who seek it. This ${\mathcal{R}}_{0}$ will typically differ from the ${\mathcal{R}}_{0}$ that would be calculated in the absence of control [14].
The reproductive numbers determine malaria disease dynamics. When ${\mathcal{R}}_{\mathit{\text{NN}}}>1$ the disease will always persist. When ${\mathcal{R}}_{\mathit{\text{NN}}}\le 1$ the disease cannot invade. However, there is evidence that under certain circumstances when ${\mathcal{R}}_{\mathit{\text{CC}}}>1>{\mathcal{R}}_{\mathit{\text{NN}}}$, clinically immune individuals can act as a reservoir and allow malaria to remain endemic, even though ${\mathcal{R}}_{\mathit{\text{NN}}}$ drops below one. In this case, both the diseasefree equilibrium and an endemic equilibrium are stable, this is an example of “bistability”.
Bistability is typically associated with “backwards bifurcations”. In general, as a disease invades, it reduces its reproductive number $\mathcal{R}$, primarily by reducing the number of susceptibles in the population. In such “forward bifurcations”, we expect the disease to go extinct from any starting conditions when ${\mathcal{R}}_{0}\le 1$, and to reach a small equilibrium, when ${\mathcal{R}}_{0}$ is just above 1 [41]. When a disease increases its reproductive number as it invades, backwards bifurcations occur. In a backwards bifurcation, the disease invades to a nonzero level even when ${\mathcal{R}}_{0}=1$, and will be able to persist above a certain threshold when ${\mathcal{R}}_{0}$ is just below 1 [41].
Methods
To explore the dynamics of malaria and determine the conditions in which bistability can occur, we evaluated the following simple transmission model (Figure 1):
S _{ N } represents naïve susceptible individuals: individuals who have never been infected with malaria, those who have been infected but have not developed clinical immunity, or who have lost all immunity. When infected with malaria, they move to the clinically infected class (I _{ N }). Recovered individuals become immediately susceptible again, but with immunity to clinical symptoms (S _{ C }). When nonnaïve susceptibles get infected, they acquire clinically immune infections (I _{ C }). Each class represents a portion of the population. τ _{ N } and τ _{ C } are the transmission rates of naïve and clinically immune cases, respectively. γ _{ NN } is the rate at which naïve individuals recover without clinical immunity, γ _{ NC } is the rate they recover with clinical immunity, and γ _{ C } is the recovery rate for clinically immune individuals. α is the rate at which clinical immunity is lost σ is the relative susceptibility to new infection of clinically immune individuals. $\Lambda =\frac{{\tau}_{N}{I}_{N}+{\tau}_{C}{I}_{C}}{T}$ is the force of infection. Both types of susceptibles (naïve or clinically immune) can be infected by either type of infectious individual (naïve or clinically immune). The total population size is T=S _{ N }+I _{ N }+S _{ C }+I _{ C }. For this model, ${\mathcal{R}}_{0}={\mathcal{R}}_{\mathit{\text{NN}}}=\frac{{\tau}_{N}}{{\gamma}_{\mathit{\text{NN}}}+{\gamma}_{\mathit{\text{NC}}}+\mu}$, and ${\mathcal{R}}_{\mathit{\text{CC}}}=\frac{\sigma {\tau}_{C}}{{\gamma}_{C}+\mu}$. We define ρ to be the ratio of the clinically immune reproductive number to the naïve reproductive number (${\mathcal{R}}_{\mathit{\text{CC}}}/{\mathcal{R}}_{\mathit{\text{NN}}}$). And $\pi =\frac{{\gamma}_{\mathit{\text{NC}}}}{{\gamma}_{\mathit{\text{NN}}}+{\gamma}_{\mathit{\text{NC}}}}$ is the proportion of naïve infections that recover to become clinically immune.
We used the statistical package R [42] to simulate our model. We held the average durations of infectiousness for naïve and immune individuals (1/γ _{ NN } and 1/γ _{ C }) constant at 50 and 200 days respectively. We chose $\pi =\frac{{\gamma}_{\mathit{\text{NC}}}}{{\gamma}_{\mathit{\text{NC}}}+{\gamma}_{\mathit{\text{NN}}}}$ to be 0.5, assuming an equal chance of getting clinical immunity and remaining nonimmune after each naïve infection; results with π=1 were similar. Transmission coefficients τ _{ N } and τ _{ C } were calculated from ${\mathcal{R}}_{\mathit{\text{NN}}}$ and ${\mathcal{R}}_{\mathit{\text{CN}}}$ which varied as described in the Results section. We chose duration of immunity 1/α to be 1282 days (about 3.5 years), and lifespan 1/μ to be 25,550 days (about 70 years). The relative susceptibility of clinically immune individuals was chosen to be σ=0.7 since clinically immune individuals are about as susceptible as naïve cases; we also tested other values.
R code and Maxima code to replicate all of our results will be made available upon publication.
In building a simple transmission model, we have lumped a variety of biological mechanisms into the parameters. For example, transmission parameters τ _{ C } and τ _{ N } include transmission blocking immunity and the reduction of parasites in clinically immune individuals. The duration of infection parameters γ _{ NN }, γ _{ NC }, and γ _{ C } take into account treatment seeking behavior (naïve individuals are likely to seek treatment quickly whereas clinically immune individuals are less likely to seek treatment or they wait longer to seek treatment).
Results
We simulated our malaria transmission model under two scenarios of infectiousness: overall infectiousness of clinically immune cases was either low ($\rho \equiv {\mathcal{R}}_{\mathit{\text{CC}}}/{\mathcal{R}}_{\mathit{\text{NN}}}=0.8$, panels 2a and 2b) or high (ρ=4, panels 2c and 2d). For each scenario, we simulated two different values of ${\mathcal{R}}_{0}\equiv {\mathcal{R}}_{\mathit{\text{NN}}}$. Diseaseinvasion results are shown in Figure 2. When we start near the diseasefree equilibrium, the qualitative behavior is determined by ${\mathcal{R}}_{0}$: when ${\mathcal{R}}_{0}>1$ (panels 2a and 2c), the disease invades and reaches an endemic equilibrium; when ${\mathcal{R}}_{0}<1$, (panels 2b and 2d), the disease does not invade.
In the case where underlying parameters can change over time, however, there are striking differences between the scenarios with low and high relative transmission from immune individuals. Figure 3 shows what would happen in a population with endemic malaria if control efforts moved transmission from the first column of Figure 2 to those shown in the second column. Panel 3a behaves as we would expect: when we change the parameters at day 3000, the system moves to the diseasefree equilibrium. Panel 3b, however, exhibits bistability. Although the parameters in the latter part of the simulation are consistent with disease extinction, the disease does not go extinct from the equilibrium reached under high transmission, but instead finds a lower endemic equilibrium. Whether or not malaria will remain endemic or die out under parameters consistent with disease extinction is dependent both on ${\mathcal{R}}_{0}$ and on initial conditions.
Figure 4 gives a broader perspective on the two scenarios, using “bifurcation diagrams” showing how equilibrium incidence changes as ${\mathcal{R}}_{0}$ increases, while holding the relative infectiousness of clinically immune individuals ($\rho \equiv {\mathcal{R}}_{\mathit{\text{CC}}}/{\mathcal{R}}_{\mathit{\text{NN}}}$) constant. Panel 4a illustrates the scenario where clinically immune individuals are relatively less effective at transmitting disease (ρ<1). Here we see the simple, common, relationship between ${\mathcal{R}}_{0}$ and equilibrium incidence. As ${\mathcal{R}}_{0}$ increases past 1, the system moves smoothly from having a globally stable equilibrium at 0, to having a globally stable endemic equilibrium.
Panel 4b shows shows the scenario where clinically immune individuals are relatively more effective at transmitting disease (ρ>1). In this case, we see a more complicated pattern, where both ${\mathcal{R}}_{0}$ and initial prevalence affect the final outcome. In particular, if we increase ${\mathcal{R}}_{0}$ smoothly past 1 in a population where the disease is absent, the equilibrium jumps abruptly; if the disease invades, and ${\mathcal{R}}_{0}$ is decreased back below 1 the disease does not necessarily go extinct. Similarly, if ${\mathcal{R}}_{\text{crit}}<{\mathcal{R}}_{0}<1$ (the light gray region of the plot), a temporary intervention that sharply reduces disease prevalence could succeed in eliminating the disease even without a longterm reduction in ${\mathcal{R}}_{0}$.
To be specific, we don’t expect the disease to go extinct once established unless ${\mathcal{R}}_{0}$ is reduced beneath the minimum value for which the endemic equilibrium exists. We call this value ${\mathcal{R}}_{\text{crit}}$ and the force of infection Λ that corresponds to it Λ _{ crit }. On the forward bifurcation diagram (Figure 4a) ${\mathcal{R}}_{\text{crit}}=1$ and on the backwards bifurcation diagram (Figure 4b) ${\mathcal{R}}_{\text{crit}}$ is precisely the point where the stability of the endemic equilibrium (in a region of bistability) changes from unstable to stable (the black dot on Figure 4b).
In Additional file 1, we show that a backwards bifurcation will occur when the ratio ${\rho}^{\ast}={\mathcal{R}}_{\mathit{\text{CC}}}/{\mathcal{R}}_{\mathit{\text{NN}}}>\rho $, where:
Here L=1/(α+μ) is the duration of immunity, D=1/(γ _{ NC }+γ _{ NN }) is the duration of naïve infection, and $\pi =\frac{{\gamma}_{\mathit{\text{NC}}}}{{\gamma}_{\mathit{\text{NC}}}+{\gamma}_{\mathit{\text{NN}}}}$ is the proportion of people who become clinically immune after a naïve infection. This criterion determines for what parameters bistability can occur when ${\mathcal{R}}_{0}<1$. The value of ρ ^{∗} is always strictly greater than 1, meaning that bistability only occurs when ${\mathcal{R}}_{\mathit{\text{CC}}}$ exceeds ${\mathcal{R}}_{\mathit{\text{NN}}}$ by a sufficient amount. The amount of excess necessary is determined by how quickly immunity is lost (through death or waning) compared to the duration of infectiousness of the disease: when immunity lasts longer, bistability is more likely.
When ρ>ρ∗, a backwards bifurcation occurs, resulting in a region of bistability where there exists a stable diseasefree equilibrium and a stable endemic equilibrium for the same parameter values. Figure 4b illustrates the backwards bifurcation and the region of bistability (light gray shaded region). Within the region of bistability, if malaria were to be eliminated, it would not be able to reinvade unless ${\mathcal{R}}_{0}$ were increased from below one back above one, making malaria elimination from these regions more sustainable. In order to eliminate malaria from a region of bistability, either the force of infection must be reduced below the unstable endemic equilibrium (dashed curve), or the reproductive number must be reduced below the critical value (${\mathcal{R}}_{\mathit{\text{crit}}}$), or a combination of these. In Figure 4b, this is equivalent to leaving the light gray shaded region (without exceeding ${\mathcal{R}}_{0}=1$, the dashed line).
We explored how different model parameters affected the bifurcation diagram (Figure 5), by varying each parameter individually. To vary ${\mathcal{R}}_{0}$, we changed τ _{ N } and τ _{ C }, while keeping their ratio constant. To vary ρ or to adjust ρ when necessary (e.g., when changing σ), we varied the ratio τ _{ N }:τ _{ C }.
We show that the region of bistability depends strongly on the ratio ρ. When ρ is large, bistability can occur even for large incidence rates, as shown in Figure 5. Consequently, if malaria were to be eliminated then reintroduced, it would jump to being endemic at a higher level in the population than when ρ is small. When ρ is large, Λ _{ crit } (the value of Λ that corresponds to ${\mathcal{R}}_{\mathit{\text{crit}}}$) is also large. A large value of Λ _{ crit } means that the force of infection needs to be reduced by less to move below the unstable endemic equilibrium and be eliminated.
We also found that even when ρ is held constant, the individual components of ρ (σ, π and γ _{ C }) still affect the bifurcation diagram, as shown in Figure 5. When we increased σ (Figure 5) while holding ρ constant, the region of bistabililty decreased. When σ is small (Figure 5) the proportion of infections that occur in clinically immune people increases, thus increasing the accumulation of clinical immunity and magnifying its effect on transmission.
We also explored the dynamics when varying the recovery rates (Figure 5), while holding ρ constant. When individuals recover quickly, the area of bistability is large, since there are more susceptible clinically immune individuals in the population. When γ _{ C } is large, Λ _{ crit } is also large; this means that the force of infection needs to be reduced by less to drop below the unstable endemic equilibrium for malaria to be eliminated. Also, when γ _{ C } is large, the value of Λ when ${\mathcal{R}}_{0}=1$ is larger than when individuals recover slowly. So if malaria is eliminated then reinvades, it would jump to a high endemic level in the population. A study of malaria resurgence found many cases where incidence jumped to a high endemic level from low (or undetectable) levels [43]. We also change the proportion of naïve infections, π. When π is large, so is ${\mathcal{R}}_{\mathit{\text{crit}}}$ and Λ _{ crit } and as π is decreased, so is ${\mathcal{R}}_{\mathit{\text{crit}}}$ and Λ _{ crit }, making the region of bistability larger.
Discussion
The question of whether there are places where malaria is endemic, but where it could remain stably eliminated under current transmission and treatment conditions (ie., places where malaria is bistable) is important to interpreting malaria data and planning control measures. It has been suggested that treatmentseeking behavior can lead to bistability in the dynamics of malaria infections of humans [12].
In areas where malaria is bistable, an aggressive program that held malaria infection at low levels until clinical immunity wanes could result in the disease remaining absent even after the program is terminated. Malaria would not reinvade in this case because treatment seeking by nonimmune infected individuals would lead to shorter duration of infectiousness and less overall transmission. This is a potentially risky strategy, however, because if malaria does reinvade such an area, the fact that clinical immunity has waned could lead to increased morbidity [9–11]. Mass drug administration (MDA) is a possible example of such an aggressive approach. Although MDA has so far proved unsuccessful in permanently interrupting malaria transmission, it is successful at reducing parasitemia and does temporarily reduce transmission [44]. Further investigation of factors underlying bistability could improve understanding of when and where MDA would be likely to lead to longterm elimination.
We analyzed a simple model, and found a simple criterion that determines whether bistability can occur. In particular, we found that the key quantity is the lifecycle “transmission effectiveness” of clinically immune individuals, relative to nonimmune individuals. We encapsulate this relative infectiousness in a ratio, which we call ρ and show that bistability can occur when ρ exceeds 1+D/(π L), where D is the duration of naïve infection, π is the proportion of naïve cases that recover to become clinically immune, and L is the length of immunity. We also show that, in addition to duration of infection and ability to transmit the disease, the relative susceptibility of clinically immune individuals to new infections is a key component of this ratio. Although the relative susceptibility is a key component to understanding bistability in malaria, little is known about the relative susceptibility of clinically immune individuals.
Chiyaka et al. [14] discuss the importance of reducing the reproductive number under control efforts (${\mathcal{R}}_{C}$) below one to eliminate malaria and to gauge the size of an outbreak arising from an imported malaria case. They also point out the importance of investigating the stability of malaria elimination. Bistability of malaria, as we explore here, is one possible mechanism that could make malaria elimination more stable. Estimates of our criterion, ρ∗ could be a valuable component of efforts to assess the stability of malaria control in specific areas.
The details of clinical immunity to malaria are more complex than those in our simple model: in particular, clinical immunity continues to develop, there is not simply one kind of clinical immunity. Nonetheless, we expect our qualitative results to apply to more realistic situations. We expect the possibility of bistability when the relative lifecycle transmission effectiveness of clinically immune individuals is high. Thus, measuring the components of transmission effectiveness, both in clinically immune and nonimmune individuals, is important for evaluating and planning malaria control efforts. Although certain aspects of malaria are well studied, it is surprisingly difficult to find information bearing directly on the components of transmission effectiveness, particularly in clinically immune individuals. Continued investigation of how these components determine transmission effectiveness will be important in understanding the populationlevel patterns of the spread and persistence of malaria.
Although duration of symptomatic infection is well studied [15, 19], little is known about duration of asymptomatic infection [23]. This is a complicated question because malaria infections can be long and variable; failure to detect parasites may not mean an infectious event is over; conversely, parasites that are detected may be due to a new infectious event.
Transmission of infection to mosquitoes is another aspect of malaria biology that is not well understood. Although they are well studied individually, it is not clear how the components of transmission come together to affect the infectiousness per transmission event. These components include: gametocytemia [28, 29]; TBI [30, 33, 34], which increases with gametocytemia, and wanes with age and clinical immunity; and other human and mosquitospecific factors [39]. More information on how these components interact to affect transmission would help to unravel how clinical immunity affects populationlevel transmission.
As we’ve shown, susceptibility of clinically immune individuals to malaria is an important component of the ratio of lifecycle transmission effectiveness. Susceptibility to new infection is known to be complex [38] but not well understood; in the literature “susceptibility” is frequently used to refer to susceptibility to clinical disease [38, 39].
Our model neglects age structure, and oversimplifies the process of clinical immunity. In particular, the population is divided simply into those who are and are not clinically immune. Other omissions include seasonality, biting heterogeneity, and in fact, mosquitoes. For these reasons, the model is not expected to provide quantitatively accurate estimates of malaria dynamics.
The advantage of this simplistic approach, in our opinion, is that the model sheds light on the possible mechanisms and key quantities underlying bistability in malaria dynamics. In particular, we expect the importance of the quantity ρ – the ratio of lifecycle effectiveness of transmission between immune and nonimmune individuals – to be robust to including more model details. Similarly, we expect some analogue of the time scale ratio D/(π L) – that is, the ratio between the time scales of infection and immunity – to be important in a detailed model.
Conclusions
It has been suggested that human malariatransmission dynamics exhibit bistability, which would have important implications for control efforts. We have shown that bistability through treatment seeking by clinically immune individuals is plausible in human malaria transmission dynamics. Using a simple model, we have demonstrated how these dynamics might play out, and determined key parameters underlying when and whether bistability might occur in real populations.
We find that the key quantities underlying whether bistability is expected to occur are: the relative “effectiveness” of clinically immune individuals, compared to naïve individuals, at transmission; and the time scale at which clinical immunity is lost, compared to the time scale of infectiousness. The model also shows that relative susceptibility to malaria infection should be considered part of transmission effectiveness, when making this comparison. We find that bistability can occur for plausible parameters, and suggest that more research into these two ratios may shed light on malaria dynamics, and guide future control efforts.
Abbreviations
 TBI:

Transmissionblocking immunity
 MDA:

Mass drug administration.
References
 1.
Organization WH: World malaria report. 2012, 1195.
 2.
Schofield L, Mueller I: Clinical immunity to malaria. Curr Mol Med. 2006, 6 (2): 205221. 10.2174/156652406776055221.
 3.
Filipe JAN, Riley EM, Drakeley CJ, Sutherland CJ, Ghani AC: Determination of the processes driving the acquisition of immunity to malaria using a mathematical transmission model. PLoS Comput Biol. 2007, 3 (12): e25510.1371/journal.pcbi.0030255.
 4.
Metselaar D, Theil PMV: Classification of malaria. Trop Geogr Malar. 1959, 11: 157161.
 5.
Gatton M, Cheng Q: Modeling the development of acquired clinical immunity to Plasmodium falciparum malaria. Infect Immun. 2004, 72 (11): 653810.1128/IAI.72.11.65386545.2004.
 6.
Gupta S, Snow R, Donnelly C, Marsh K, Newbold C: Immunity to noncerebral severe malaria is acquired after one or two infections. Nature Med. 1999, 5 (3): 340343. 10.1038/6560.
 7.
Doolan DL, Dobano C, Baird JK: Acquired immunity to malaria. Clin Microbiol Rev. 2009, 22: 1336. 10.1128/CMR.0002508.
 8.
Kurtis J, Mtalib R, Onyango F, Duffy P: Human resistance to Plasmodium falciparum increases during puberty and is predicted by dehydroepiandrosterone sulfate levels. Infect Immun. 2001, 69: 12310.1128/IAI.69.1.123128.2001.
 9.
Snow R, Marsh K: Will reducing Plasmodium falciparum transmission alter malaria mortality among African children?. Parasitol Today. 1995, 11 (5): 188190. 10.1016/01694758(95)801561.
 10.
OMeara WP, Mwangi TW, Williams TN, McKenzie FE, Snow RW, Marsh K: Relationship between exposure, clinical malaria, and age in an area of changing transmission intensity. Am J Trop Med Hyg. 2008, 79 (2): 185191.
 11.
O’Meara WP, Bejon P, Mwangi TW, Okiro EA, Peshu N, Snow RW, Newton CR, Marsh K: Effect of a fall in malaria transmission on morbidity and mortality in Kilifi, Kenya. The Lancet. 2008, 372 (9649): 15551562. 10.1016/S01406736(08)616554.
 12.
Águas R, White LJ, Snow RW, Gomes MGM: Prospects for malaria eradication in subSaharan Africa. Plos One. 2008, 3 (3): 16.
 13.
The Global HealthGroup, the Malaria AtlasProject: Atlas of malaria eliminating countries. 2011, San Francisco:The Global Health Group, Global Health Sciences, University of California, San Francisco
 14.
Chiyaka C, Tatem A, Cohen J, Gething P, Johnston G, Gosling R, Laxminarayan R, Hay S, Smith D: The stability of Malaria elimination. Science. 2013, 339 (6122): 909910. 10.1126/science.1229509.
 15.
McCombie S: Treatment seeking for malaria: a review of recent research. Soc Sci Med. 1996, 43 (6): 933945. 10.1016/02779536(95)004467.
 16.
OwusuAgyei S, Koram K, Baird J, Utz G, Binka F, Nkrumah F, Fryauff D, Hoffman S: Incidence of symptomatic and asymptomatic Plasmodium falciparum infection following curative therapy in adult residents of northern Ghana. Am J Trop Med Hygiene. 2001, 65 (3): 197203.
 17.
Port AL, Cot M, Etard JF, Gaye O, MigotNabias F, Garcia A: Relation between Plasmodium falciparum asymptomatic infection and malaria attacks in a cohort of Senegalese children. Malar J. 2008, 7: 19310.1186/147528757193.
 18.
Bousema J, Gouagna L, Drakeley C, Meutstege A, Okech B, Akim I, Beier J, Githure J, Sauerwein R: Plasmodium falciparum gametocyte carriage in asymptomatic children in western Kenya. Malar J. 2004, 3: 1810.1186/14752875318.
 19.
Earle W, Perez M, del Rio J, Arzola C: Observations on the course of naturally acquired malaria in Puerto Rico. Puerto Rico J Public Heath and Trop Med. 1939, 14: 391406.
 20.
Sama W, Killeen G, SMITH T: Estimating the duration of Plasmodium falciparum infection from trials of indoor residual spraying. Am J Trop Med Hygiene. 2004, 70 (6): 625
 21.
Bretscher M, Maire N, Chitnis N, Felger I, OwusuAgyei S, Smith T: The distribution of Plasmodium falciparum infection durations. Epidemics. 2011, 3 (2): 109118. 10.1016/j.epidem.2011.03.002.
 22.
Szmitko PE, Kohn ML, Simor AE: Plasmodium falciparum malaria occurring 8 years after leaving an endemic area. Diagn Microbiol Infect Dis. 2008, 63: 105107.
 23.
Bruce M, Donnelly C, Packer M, Lagog M, Gibson N, Narara A, Walliker D, Alpers M, Day K: Ageand speciesspecific duration of infection in asymptomatic malaria infections in Papua New Guinea. Parasitol. 2000, 121 (03): 247256. 10.1017/S0031182099006344.
 24.
Molineaux L, Diebner H, Eichner M, Collins W, Jeffery G, Dietz K: Plasmodium falciparum parasitaemia described by a new mathematical model. Parasitol. 2001, 122 (04): 379391.
 25.
Bottius E, Guanzirolli A, Trape JF, Rogier C, Konate L, Druilhe P: Malaria: even more chronic in nature than previously thought; evidence for subpatent parasitaemia detectable by the polymerase chain reaction. Trans R Soc Trop Med Hygiene. 1996, 90: 1519. 10.1016/S00359203(96)904630.
 26.
Alves F, Durlacher R, Menezes M, Krieger H, Silva L, Camargo E: High prevalence of asymptomatic Plasmodium vivax and Plasmodium falciparum infections in native Amazonian populations. Am J Trop Med Hygiene. 2002, 66 (6): 641648.
 27.
Laishram DD, Sutton PL, Nanda N, Sharma VL, Sobti RC, Carlton JM, Joshi H: The complexities of malaria disease manifestations with a focus on asymptomatic malaria. Malar J. 2012, 11: 2910.1186/147528751129.
 28.
Schall J: Transmission success of the malaria parasite Plasmodium mexicanum into its vector: role of gametocyte density and sex ratio. Parasitol. 2000, 121 (06): 575580.
 29.
Noden B, Beadle P, Vaughan J, Pumpuni C, Kent M, Beier J: Plasmodium falciparum: the population structure of mature gametocyte cultures has little effect on their innate fertility. Acta Tropica. 1994, 58: 1319. 10.1016/0001706X(94)901171.
 30.
Drakeley CJ, Bousema JT, Akim NIJ, Teelen K, Roeffen W, Lensen AH, Bolmer M, Sauerwein WERW: Transmissionreducing immunity is inversely related to age in Plasmodium falciparum gametocyte carriers. Parasite Immunol. 2006, 28: 185190. 10.1111/j.13653024.2005.00818.x.
 31.
Naotunne T, Karunaweera N, Giudice GD, Kularatne M, Grau G, Carter R, Mendis K: Cytokines kill malaria parasites during infection crisis extracellular complementary factors are essential. J Exp Med. 1991, 173 (3): 523529. 10.1084/jem.173.3.523.
 32.
Eckhoff P: P. falciparum infection durations and infectiousness are shaped by antigenic variation and innate and adaptive host immunity in a mathematical model. PLoS ONE. 2012, 7 (9): e4495010.1371/journal.pone.0044950.
 33.
Mulder B, Tchuinkam T, Dechering K, Verhave J, Carnevale P, Robert V, JHET M: Malaria transmissionblocking activity in experimental infections of Anopheles gambiae from naturally infected Plasmodium falciparum gametocyte carriers. Trans R Soc Trop Med Hygiene. 1994, 88 (1): 121125. 10.1016/00359203(94)905347.
 34.
Boudin C, Kolk MVD, Tchuinkam T, Gouanga C, Bonnet S, Safeukui I, Mulder B, Meunier JY, Verhave JP: Plasmodium Falciparum transmission blocking immunity under conditions Of low and high endemicity in Cameroon. Parasite Immunol. 2004, 26: 105110. 10.1111/j.01419838.2004.00689.x.
 35.
Killeen G, Ross A, Smith T: Infectiousness of malariaendemic human populations to vectors. Am J Trop Med Hygiene. 2006, 75 (2 suppl): 3845.
 36.
Ross A, Killeen G, Smith T: Relationships between host infectivity to mosquitoes and asexual parasite density in Plasmodium falciparum. Am J Trop Med Hygiene. 2006, 75 (2 suppl): 3237.
 37.
Pethleart A, Prajakwong S, Suwonkerd W, Corthong B, Webber R, Curtis C: Infectious reservoir of Plasmodium infection in Mae Hong Son province, northwest Thailand. Malar J. 2004, 3: 3410.1186/14752875334.
 38.
MigotNabias F, Luty A, Ringwald P, Vaillant M, Dubois B, Renaut A, Mayombo R, Minh T, Fievet N, Mbessi J: Immune responses against Plasmodium falciparum asexual bloodstage antigens and disease susceptibility in Gabonese and Cameroonian children. Am J Trop Med Hygiene. 1999, 61 (3): 488494.
 39.
Draper C: Observations on the infectiousness of gametocytes in hyperendemic malaria. Trans R Soc Trop Med Hygiene. 1953, 47 (2): 160165. 10.1016/00359203(53)900728.
 40.
Killeen G, Ross A, Smith T: Infectiousness of malariaendemic human populations to vectors. Am J Trop Med hygiene. 2006, 75 (2 suppl): 3845.
 41.
Dushoff J: Incorporating immunological ideas in epidemiological models. J Theor Biol. 1996, 180: 181187. 10.1006/jtbi.1996.0094.
 42.
R Core Team: R: A Language and Environment for Statistical Computing. 2012, Vienna, Austria: R Foundation for Statistical Computing, [http://www.Rproject.org]
 43.
Cohen JM, Smith DL, Cotter C, Ward A, Yamey G, Sabot OJ, Moonen B: Malaria resurgence: a systematic review and assessment of its causes. Malar J. 2012, 11: 11. 10.1186/14752875111.
 44.
von Seidlein L, Greenwood B: Mass administrations of antimalarial drugs. Trends in Parasitol. 2003, 19 (10): 452460. 10.1016/j.pt.2003.08.003.
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712334/13/428/prepub
Acknowledgements
The authors thank Jake Szamosi, Marta Wayne, Ben Bolker, Mark Loeb and David Earn for useful comments and suggestions. We also thank David Smith for helpful discussions about framing the research project.
Author information
Additional information
Competing interests
Both authors declare that they have no competing interests.
Author’s contributions
LK and JD made and analyzed the model and drafted the manuscript. Both authors have read and approved the manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
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.
About this article
Received
Accepted
Published
DOI
Keywords
 Malaria
 Malaria Infection
 Mass Drug Administration
 Endemic Equilibrium
 Malaria Elimination