Modelling knowlesi malaria transmission in humans: vector preference and host competence

Abstract

Background

Plasmodium knowlesi, a malaria species that normally infects long-tailed macaques, was recently found to be prevalent in humans in Southeast Asia. While human host competency has been demonstrated experimentally, the extent to which the parasite can be transmitted from human back to mosquito vector in nature is unclear.

Methods

Using a mathematical model, the influence of human host competency on disease transmission is assessed. Adapting a standard model for vector-borne disease transmission and using an evolutionary invasion analysis, the paper explores how differential host competency between humans and macaques can facilitate the epidemiological processes of P. knowlesi infection between different hosts.

Results

Following current understanding of the evolutionary route of other human malaria vectors and parasites, an increasing human population in knowlesi malaria endemic regions will select for a more anthropophilic vector as well as a parasite that preferentially transmits between humans. Applying these adaptations, evolutionary invasion analysis yields threshold conditions under which this macaque disease may become a significant public health issue.

Conclusions

These threshold conditions are discussed in the context of malaria vector-parasite co-evolution as a function of anthropogenic effects.

Background

Although principally a disease of long-tailed and pig-tailed macaques, Plasmodium knowlesi has recently been revealed to be a widespread, and potentially life-threatening, malaria infection of humans in Southeast Asia [1]. However, before the public health threat that this disease poses can be assessed, a solid understanding of the fundamental epidemiological processes is needed. Using a simple mathematical model, key components of the disease are explored that are of particular influence in its transmission and yet for which data are lacking, thereby highlighting priority research areas.

Through a series of laboratory experiments, Chin et al demonstrated the ability of this parasite to transmit from humans to both simian and human hosts [2]. However, the extent of human host competency under natural conditions remains unclear, giving rise to controversy over whether P. knowlesi can be classed as the fifth aetiological agent of human malaria [3–5]. Implicit to the epidemiological understanding of this disease is a formalized description of its transmission dynamics. Mathematical methods of Macdonald [6] are, therefore, adapted to allow for multiple host species and used to determine the significance of human host competency in disease transmission.

Traditionally, mosquitoes with higher anthropophily (such as Anopheles gambiae) are considered more efficient vectors of human malaria [7]. This is because blood meals taken from alternative mammalian hosts reduce the transmission intensity of human disease. Therefore, the second component considered in this analysis is the influence of vector host preference on knowlesi malaria transmission dynamics under the condition that humans are either non-competent (dead-end) or competent hosts.

In the case where humans are competent hosts for P. knowlesi, evolutionary theory would suggest that selection favours a parasite that has co-evolved with its natural host, and this has important implications for among-species transmission potential. However, the limited immunological data on this pathogen suggests that serial passage of P. knowlesi through humans can result in increasingly higher parasite densities, until this pathogen reaches life-threatening levels [8]. In this analysis, the consequences of parasite life history that affect infectivity of knowlesi malaria to humans are determined. Using evolutionary invasion analysis, the conditions necessary for a switch of evolutionary strategies in favour of human knowlesi malaria transmission are calculated. Results are discussed in the evolutionary context of all human malarias and the implications to all vector-borne pathogens are described.

Methods

The dynamics are described by a set of ordinary differential equations, adapted from Smith et al (2007) [9] to include multiple hosts (e.g, humans and macaques). The model is partitioned into three equations; one that represents the proportion of infected humans ($\dot{H}$):

$$\frac{d\dot{H}}{dt}=\frac{V}{H}bh{c}_{VH}\dot{V}(1-\dot{H})-r\dot{H}$$

(1)

one that describes the proportion of infected macaques ($\dot{M}$):

$$\frac{d\dot{M}}{dt}=\frac{V}{M}bm{c}_{VM}\dot{V}(1-\dot{M})-r\dot{M}$$

(2)

and finally one that describes the proportion of infected vectors ($\dot{V}$):

$$\frac{d\dot{V}}{dt}=b\left(\dot{H}h{c}_{HV}+\dot{M}m{c}_{MV}\right)({\exp }^{-gT}-\dot{V})-g\dot{V}$$

(3)

where V, H and M are the vector, human and macaque density; b is the bite rate; g is the daily vector mortality rate; T is the extrinsic incubation period; r is the rate of mammalian recovery from infection; c _VH and c _VM are the parasite transmission coefficients from vector to human and from vector to macaque respectively; c _HV and c _MV are the parasite transmission coefficients from human to vector and from macaque to vector respectively; h and m are the proportion of bites that fall on humans and macaques respectively (see below). The vectorial capacity (the capacity of the vector to transmit infection to humans) is determined from:

$$Vc=\frac{\frac{V}{H}{b}^{2}h\left(h{C}_H+m{C}_M\right)\exp \left(-gT\right)}{g}$$

(4)

where, C _H and C _M is the parasite transmission efficacy from human-to-human (the product of c _HV and c _VH ) and macaque-to-human (the product of c _MV and c _VH ), respectively. The basic reproductive number of malaria, R ₀ , is simply the product of Vc with the average rate of human recovery (1/r).

Following Smith et al. (2007) [9], the infected proportions of humans ($\overline{H}$), macaques ($\overline{M}$) and vectors ($\overline{V}$) at equilibrium are:

$$\overline{H}=\frac{R0-1}{R0+bh{c}_{HV}{g}^{-1}}$$

(5)

$$\overline{M}=\frac{R0-1}{R0+bm{c}_{MV}{g}^{-1}}$$

(6)

$$\overline{V}=\frac{b\left(\overline{M}m{c}_{MV}+\overline{H}h{c}_{HV}\right)}{g+b\left(\overline{M}m{c}_{MV}+\overline{H}h{c}_{HV}\right)}\exp \left(-gT\right).$$

(7)

As noted, h and m are the proportion of bites that fall on humans and macaques respectively and are calculated as follows:

$$h=HQ/\left(M(1-Q)+HQ\right)$$

(8)

$$m=M(1-Q)/\left(M(1-Q)+HQ\right).$$

(9)

In this way, the proportion of bites made on alternative hosts accounts for vector preference for human blood (0≤Q≤1) as well as the relative host densities. Using these equations, the relationship between vectorial capacity and vector preference for human blood is examined for systems in which humans are either dead-end (C _H = 0) or competent hosts (C _H > 0).

While the parasite transmission efficacy between vector and human hosts is well established for true human malarias, data for knowlesi parasite transmission are lacking. Using classic parameter values that have been recorded for human malarias (daily mosquito mortality of 0.15 [10], extrinsic incubation period of 10 days, daily bite rate of 0.25 and transmission efficacy from human to vector to human of 0.05 [11]) the infection dynamics are simulated for variable parasite infectivity as a function of alternative transmission routes (via macaque or via human). There are two opposing mechanisms by which infectivity might be expected to differ according to transmission routes: 1) given it has co-evolved with macaques, more effective parasite transmission from macaque to mosquito to human might be expected, but 2) a parasite that has already been passaged through a human immune system might be expected to have greater infectivity to a secondary human. Therefore, the dynamics of both alternative assumptions are presented.

The ultimate aim of this study is to understand what conditions an alternative vector could evolve (invade). Evolutionary invasion analysis explores where the resident vector, human, macaque populations are extant and at equilibrium. Novel strategies invade from rare (V _novel density is very small compared to the resident vector).

Results

In the situation where humans are non-competent (dead-end) hosts, vector preference for human hosts has a non-linear relationship with the vectorial capacity (Figure 1A). Intuitively, at the two extremes in human host preference (Q = 0, Q = 1), a zoonosis cannot persist in humans. In this system, Plasmodium knowlesi transmission in humans is dependent on relative mammalian host densities, with the greatest detriment to human health experienced when there is an even distribution of bites between both humans and macaques.

Figure 1

(A) Vectorial Capacity as a function of human host preference when humans are dead-end hosts (i.e. C _H = 0). The black line corresponds with equal human and macaque density (H:M); blue with H:1/3M; red with H:3M. (B) Vectorial Capacity as a function of human host preference when humans are competent hosts. Colours correspond with (A), broken lines denote vectors that are 'exclusive' in their host choice and solid lines denote 'mixed' preferences. (C) The effect of variable parasite infectivity as a function of its transmission route. Purple corresponds to a parasite that transmits human to mosquito to human more effectively (line, C_H:1/2C_M; cross markers C_H:1/4C_M), and green corresponds to a parasite that transmits macaque to mosquito to human more effectively (line, C_H:2C_M; cross markers C_H:4C_M). Unless stated otherwise we use H:10V; H:M; C _H :C _M (= 0.05); b = 0.25; g = 0.15 and T = 10.

On analysing the dynamics of a system in which humans are competent hosts, the necessity of further scrutinizing the definition of 'host preference' becomes apparent. Under equivalent host preference (Q = 0.5), one extreme scenario is that each individual vector exhibits identically indiscriminate host preference ('mixing vectors') while the alternative extreme entails half of the vectors always choosing one host species over the other ('exclusive vectors'). This latter scenario would mirror the dynamics of a traditional human malaria system because a mosquito that only bites humans cannot be infected by macaques. In general, mosquitoes that are exclusive in their host preference are more effective vectors of malaria (Figure 1B).

The relationship between host preference and vectorial capacity is highly dependent on the parasite response to alternative transmission routes (Figure 1C). Consistent with current understanding of human malaria transmission, results show that an exclusively anthropophilic mosquito (Q = 1) is usually the optimum vector of malaria in humans when the pathogen no longer requires an alternative mammalian host to complete its lifecycle. However, a special case arises whereby the parasite pervades the human population more readily when transmitted by a generalist vector (Figure 1C). For a generalist mosquito to be a superior malaria vector, the following condition must be satisfied:

$$\frac{C_H}{C_M}<\frac{Q}{Q+1}$$

(10)

Parasite transmission from macaque-to-humans must be more than twice as efficient as human-to-human transmission in order to favour a generalist vector. Alternatively, in the case of an uneven density of humans and macaques:

$$\frac{C_H}{C_M}<\frac{HQ}{M(1-Q)+2HQ}$$

(11)

Hence, conditions are even more permissive for enhanced transmission despite incomplete anthropophily when there is a greater difference in infectivity and when human population density exceeds that of macaques.

Using evolutionary invasion analysis, the fitness (per capita rate of change of a strategy) dynamics of a novel malaria vector (V _novel ) are given by:

$$\begin{array}{l}F=\frac{1}{{\dot{V}}_{novel}}\frac{d{\dot{V}}_{novel}}{dt}=\\ \frac{b^{\prime }\left(\overline{M}{m}^{\prime }{c^{\prime }}_{MV}+\overline{H}{h}^{\prime }{c^{\prime }}_{HV}\right)({\exp }^{-g^{\prime }{T}^{\prime }}-{\dot{V}}_{novel})-g^{\prime }{\dot{V}}_{novel}}{{\dot{V}}_{novel}}\end{array}$$

(12)

where the primes denote parameters associated with this novel strategy. This strategy has positive fitness if F > 0, which occurs if:

$$\begin{array}{l}{b}^{\prime }\left(\overline{M}{m}^{\prime }{c^{\prime }}_{MV}+\overline{H}{h}^{\prime }{c^{\prime }}_{HV}\right)({\exp }^{-g^{\prime }{T}^{\prime }}-{\dot{V}}_{novel})\hfill \\ -g^{\prime }{\dot{V}}_{novel}>0\hfill \end{array}.$$

(13)

This leads to:

$$b^{\prime }\left(\overline{M}{m}^{\prime }{c^{\prime }}_{MV}+\overline{H}{h}^{\prime }{c^{\prime }}_{HV}\right)>\frac{g^{\prime }{\dot{V}}_{novel}}{({\exp }^{-g^{\prime }{T}^{\prime }}-{\dot{V}}_{novel})}.$$

(14)

This fitness threshold can be made explicit for any of the novel vector strategy parameters. This study examines the vector biting behaviour that will allow an alternative parasite transmission cycle to evolve. The proportion of bites on humans (h') by the novel vector strategy must exceed losses due to mortality and the proportion of bites on macaques:

$$h^{\prime }>\frac{g^{\prime }{V}_{novel}+{c^{\prime }}_{MV}\overline{M}{m}^{\prime }\left({\dot{V}}_{novel}-\exp (-g^{\prime }{T}^{\prime })\right)}{C_{HV}\overline{H}\left(\exp (-g^{\prime }{T}^{\prime })-{\dot{V}}_{novel}\right)}.$$

(15)

The evolutionary stable strategy (ESS) with respect to h' (the proportion of bites on humans by the novel strategy) is found from setting the partial derivative of equation 12 with respect to h' to equal zero:

$$\frac{\partial F}{\partial h^{\prime }}=\frac{b^{\prime }{c^{\prime }}_{HV}\overline{H}({\exp }^{-g^{\prime }{T}^{\prime }}-{\dot{V}}_{novel})}{{\dot{V}}_{novel}}=0.$$

(16)

This can be simplified to show that at the ESS associated with the proportion of bites on humans, the average mosquito lifespan is:

$$\frac{1}{{g^{\prime }}_{ESS}}=\frac{T^{\prime }}{-\ln ({\dot{V}}_{novel})}.$$

(17)

As noted, the novel vector strategy is rare so, conservatively, vectorial capacity (equation 7) of this malaria vector is:

$$V{c}_{novel}=\frac{\frac{1}{H}b{\text{'}}^{2}h\text{'}\left(h\text{'}C{\text{'}}_H+m\text{'}C{\text{'}}_M\right)\exp \left(-g\text{'}T\text{'}\right)}{g\text{'}}.$$

(18)

Successful spread of this malaria pathogen by an alternative vector occurs if Vc _novel > Vc which is true if:

$$\begin{array}{l}\frac{\frac{1}{H}b{\text{'}}^{2}h\text{'}\left(h\text{'}C{\text{'}}_H+m\text{'}C{\text{'}}_M\right)\exp \left(-g\text{'}T\text{'}\right)}{g\text{'}}>\hfill \\ \frac{\frac{V}{H}{b}^{2}h\left(h{C}_H+m{C}_M\right)\exp \left(-gT\right)}{g}\hfill \end{array}.$$

(19)

Other things equal, spread of the pathogen by the alternative vector occurs if $\frac{1}{H}>\frac{V}{H}$; which is unlikely to be true as the density of resident vectors (V) is always likely to exceed the alternative vector (V'). Greater bite rates (b'), attacks on humans (h') or increased longevity (1/g') (see equation 17) favour greater vectorial capacity.

Further, if it assumed that the density of humans (H) is constant over the evolutionary ecological changes associated with the vector then the invasion potential of the alternative parasite transmission cycle (via humans) can be evaluated. This is illustrated in Figure 2. Novel strategies are more likely to evolve if the proportion of bites on humans exceeds a critical threshold value (Figure 2). Specifically, if resident vectors of the parasite have equal preference for humans and macaques, novel strategies evolve if they have a greater than 50% preference for human hosts. Similarly, even if malaria vectors have no or low preference for humans (h→0) then the proportion of bites on humans by the novel vector (h') needs to be substantial (greater than 40%) but not necessarily requiring that humans are the preferred host in order for this malaria vector to spread (Figure 2).

Figure 2

Fitness boundary of the proportion of bites on humans by a novel vector necessary to exclude the resident vector. Novel vectors need to achieve a threshold level of bites on humans (at least 40%) in order to evolve (invade). Fitness is determined from equation 12 and represents the per capita growth rate of the novel vector strategy in the presence of alternative vectors and hosts when at equilibrium.

Discussion

Human-to-human transmission of P. knowlesi has been shown to occur under experimental conditions [2], and with human populations increasingly encroaching on and supplanting macaque habitat, there is powerful selective pressure for the parasite to switch its natural mammalian host. Plasmodium vivax, the most prominent human malaria species in SE Asia and the closest relative to P. knowlesi, is believed to have a similar evolutionary route [12]. Whether this switch is already occurring in knowlesi malaria is currently unclear [13].

Concurrent with a parasite that is potentially switching host due to increased exposure to a human population, specialization of the associated vector might be expected for the same reason. For example, isomorphic members of the Anopheles gambiae species complex have diverse feeding patterns, several of which have rapidly co-evolved with the relatively recent anthropogenic practices of agriculture and permanent settlement [14, 15]. Speciation, such as that seen in the Anopheles gambiae complex, would be expected to benefit the persistence of the knowlesi malaria parasite as it becomes more dependent on the increasingly prevalent human host population. This analysis highlights the critical threshold conditions for a switch in evolutionary stable strategies.

Here, two distinct mechanisms of measuring host selectivity have been described: 'human host preference' is the proportion of bites that fall on humans, averaged over the population of vectors, whereas 'exclusivity' refers to the individual mosquito level of preference. This is an important distinction because when between-individual differences in host choice are taken into account, exclusive mosquitoes normally constitute better vectors of human infection than mixing vectors. For example, when 20% of the total mosquito blood meals are taken from humans, the parasite pervades the human population more effectively if it is the same 20% of mosquitoes always biting humans than if it is all mosquitoes biting humans 20% of the time. This is analogous to the finding that heterogeneous bite rates across a population enhances the persistence of a vector-borne disease, as found by other studies [16–18]. Whereas the heterogeneity was considered in human "attractiveness" or availability previously, in this analysis, heterogeneity in vector preferences is examined. Also important in malaria transmission are heterogeneities in the susceptibility of hosts to infection [19]. Empirical studies have suggested that these sorts of heterogeneities are important in human infection with knowlesi [20]. Once the extent of heterogeneity of macaque infection is characterised, the significance of this addition can be estimated using the framework developed in this study.

Strong anthropophily is normally considered a prerequisite for an effective malaria vector in humans [7]. In this respect, P. knowlesi might be expected to differ from true human malarias because its spread to the human population necessitates a generalist vector under the circumstances that humans are dead-end hosts. Interestingly, even if the parasite no longer requires the macaque host to persist, this analysis demonstrates how a generalist mosquito may still prove to be a more effective vector of malaria in humans.

Conclusions

Accurately assessing the threat to public health that this pathogen poses clearly necessitates extensive empirical work. The problem is confounded further by the fact that there are numerous vector species with extremely variable anthropophily, and many simian host species [21, 22]. Even in the simplified description of the system presented here, qualitative differences arose from alternative assumptions of human host competency, emphasizing this to be a critical research area. Additionally, given that most vector-borne human infections are zoonoses, clarifying ambiguous definitions of host preference is also a matter of profound public health importance.